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INTRODUCTION 


1 . 1 BACKGROUND 

The  "energy  crisis"  has  caused  virtually  every 
segment  of  American  society  to  reconsider  energy-use  poli- 
cies. The  immediate  concern  for  running  out  of  fuel  has 
lessened,  but  we  are  left  with  greatly  increased  fuel  costs, 
the  continuing  threat  of  localized  shortfalls,  and  the  reminder 
that  natural  sources  of  hydrocarbon  fuels  eventually  will  run 
dry . 


Throughout  recent  developments,  the  military  has  been 
guaranteed  a war  reserve  of  fuel,  and  the  nation's  defense 
posture  has  not  been  compromised  significantly.  Nevertheless, 
fuel-use  efficiency  is  a central  issue  of  military  prepared- 
ness for  several  reasons.  First,  there  are  obvious  tactical 
and  strategic  reasons:  if  the  average  specific  fuel  consump- 

tion of  a military  aircraft  can  be  improved,  it  can  fly 
farther  and/or  carry  a bigger  payload.  If  the  net  fuel  use 
per  mission  is  reduced,  more  missions  can  be  flown  with  a 
fixed  amount  of  fuel.  From  a budgetary  viewpoint,  the  cost  of 
fuel  per  mission  is  another  important  factor,  with  per-gallon 
costs  averaging  3 to  4 times  what  they  were  a year  ago. 

Impacts  of  fuel  shortages  on  training  schedules,  crew  pro- 
ficiency, reserve-unit  support,  trans-shipment  and  utility 
flying,  and  aircraft  deployment  also  are  major  issues.  Fur- 
thermore, in  times  of  national  emergency,  the  availability  of 
the  war  reserve  of  fuel  itself  could  be  jeopardized. 
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The  seriousness  of  this  situation  suggests  not  only 
that  future  aircraft  be  designed  for  increased  fuel  efficiency, 
but  that  fuel-optimality  of  the  operations  of  current  aircraft 
be  examined  in  detail.  Although  numerous  studies  of  improved 
aircraft  operations  have  been  made,  few  comprehensive  evalua- 
tions of  fuel-optimal  climb-to-descent  flight  profile's  have 
been  conducted  for  existing  aircraft.  The  technology  for 
conducting  such  studies  exists  and  has  been  demonstrated  for 
individual  flight  profiles,  but  it  has  rarely  been  applied 
in  such  a way  that  sensitivities  to  mission  variations  and 
parameter  uncertainties  are  fully  understood. 

In  the  past,  aircraft  fuel  efficiency  often  has 
been  considered  to  be  synonymous  with  cruise  efficiency,  as 
essential  efforts  in  designing  aircraft  have  centered  on 
maximizing  range  (Refs.  1 to  5).  During  maximum-range  flight, 
a transport-  or  bomber-type  aircraft  spends  a ma.jor  portion 
of  its  time  at  the  cruise  condition,  and  theie  is  strong 
justification  for  the  point  of  view  that  minimizing  fuel  us>>d 
during  cruise  minimizes  overall  fuel  use.  When  day-to-day 
operations  of  such  aircraft  are  examined,  however,  it  is 
found  that  average  stage  lengths  ("hops")  are  substantially 
below  maximum  range.  Any  practical  approach  to  saving  fuel 
must  consider  the  reduced  range  of  flight  involved  in  normal 
transport  and  training  operations,  as  well  as  the  maximum 
range  of  the  aircraft. 

In  recent  years,  research  emphasis  has  shifted  from 
range-maximizing  equilibrium  flight  ("Breguet  Cruise")  to 
fuel-minimizing  dynamic  flight  ( "Energy- State"  climb  and 
descent),  and,  in  fact,  the  fuel-optimality  of  cruise  itself 
has  been  questioned  (Ref.  6)  and  defended  (Ref.  7).  Building 
on  the  "energy-maneuverability”  and  "excess  power"  concepts 
introduced  in  the  1950’s,  approximate  techniques  for  analyzing 
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minimum-time-to-climb  and  minimum- fuel -to-c limb  problems 
were  developed  (Refs.  8 to  17).  These  provide  valuable  in- 
sights concerning  the  importance  of  managing  kinetic  and 
potential  energy,  taking  account  of  factors  such  as  specific 
fuel  consumption,  altitude-Mach  number  effects  of  thrust,  and 
Mach-dependent  drag  rise  in  shaping  climb  profiles.  For 
fuel-minimizing  climb,  these  profiles  consist  of  two  types 
of  arcs:  those  which  maximize  total  energy  gained  per  unit 

of  fuel  expended  and  "zoom"  climb  or  dive  segments  which 
provide  kinetic-potential  energy  interchange.  Recent  develop- 
ments have  been  directed  at  computing  and  matching  arcs  using 
asymptotic  expansions  and  singular  perturbation  theory  (Refs. 
18  to  21),  and  practical  approaches  to  specifying  optimal 
flight  path  segments  are  treated  in  Refs.  22  and  23. 

These  approaches  have  identified  the  character  of 
fuel-optimal  flight  paths,  although,  in  most  instances,  the 
comparison  between  approximate  and  exact  results  (obtained  by 
numerical  optimization)  has  not  been  pursued  extensively. 

The  present  study  seeks  to  fill  this  gap  by  investigation  of 
fuel-optimal  flight  paths  which  contain  climb,  cruise,  and 
descent  segments,  which  are  computed  with  a detailed  mathe- 
matical model  of  an  existing  aircraft,  and  which  could  be 
repeated  in  actual  aircraft  operations. 


1.2  SCOPE  AND  TECHNICAL  APPROACH 

This  report  presents  an  analysis  of  the  application 
of  integrated  throttle/flight  path  energy  management  techniques 
to  a military  transport  aircraft.  The  purpose  is  to  establish 
quantitative  data  on  fuel  economies  and  range  extension 
achievable  for  typical  mission  scenarios.  Specific  algorithms 
designed  to  optimize  throttle  setting  and  flight  trajectory 
within  operational  mission  constraints  are  synthesized  and 
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evaluated.  An  on-board  system  to  Implement  the  energy 
management  functions  in  real  time  is  defined,  and  a test  con- 
cept for  system  validation  and  flight  test  is  presented,  con- 
sidering alternative  manual  and  automatic  control  modes. 

Unified  guidance  algorithms  are  developed  using  a 
common  conceptual  and  methodological  approach,  namely  that 
associated  with  modern  control  theory.  The  design  technique 
can  be  summarized  as  follows: 

• Compute  nonlinear,  fuel-optimal  trajectories 
for  a specific  aircraft  type  and  typical 
mission  scenarios. 

• Compute  linear-optimal  perturbation  guidance 
gains  along  these  fuel-optimal  trajectories. 

• Integrate  the  reference  trajectories  and 
perturbation  guidance  gains  to  form  an  on- 
board throttle/energy  management  algorithm. 

Further  details  of  this  technical  approach  can  be  found  in  the 

remainder  of  the  report.  Factors  which  may  constrain  fuel 

saving  measures  (e.g.,  air  traffic  control  and  mission 

scenarios)  are  identified,  and  the  specific  "ground  rules” 

for  system  design  are  established,  using  the  C-141A  transport 

aircraft  as  a basis  for  detailed  analysis. 


1.3  SUMMARY  OF  RESULTS 

The  major  tasks  of  this  project  were  defined  at  the 
outset  as: 

• Aircraft  Selection 

• Constraint  Definition 

• Algorithm  Development 

• Concept  Validation 

• System  Synthesis 

• Test  Concept 
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These  tasks  are  summarized  briefly.  Aircraft  Selection 
required  the  choice  of  a specific  aircraft  for  detailed 
study,  together  with  the  reasons  for  making  the  selection. 
Constraint  Definition  identified  the  operational  limitations 
and  control  objectives  to  be  applied  in  the  investigation. 
Algorithm  Development  provided  throttle/energy  management 
equations  which  were  considered  most  suitable  for  achieving 
the  stated  objectives.  Concept  Validation  was  achieved  by 
a quantitative  comparison  of  fuel-optimal  and  conventional 
flight  policies.  System  Synthesis  was  directed  at  design 
for  real-time  implementation  of  the  energy  management 
algorithms.  Test  Concept  developed  a methodology  for  plan- 
ning and  executing  flight  tests  to  demonstrate  flight  policies 
for  fuel  conservation. 

The  C-141A  jet  transport  was  chosen  as  a baseline 
for  study.  This  aircraft  type  constitutes  the  single  greatest 
user  of  Jet  fuel  in  the  Air  Force,  and  its  flight  operations 
are  flexible  enough  that  substantial  variations  in  vertical 
profile  could  be  flown  without  altering  mission  constraints 
and  objectives. 

Fuel  savings  from  a few  percent  to  more  than  20 
percent  can  be  realized  by  reshaping  the  aircraft’s  vertical 
flight  profile.  The  actual  saving  depends  on  stage  length 
and  the  flight  profile  taken  as  a reference  (the  latter 
ambiguity  is  eliminated  by  choosing  the  optimal  profile  as 
a base  and  evaluating  the  excess  fuel  use  of  various  con- 
ventional flight  paths);  however,  the  comparison  with  hand- 
book profiles  and  cruise  points  used  for  pre-flight  planning 
suggests  that  substantial  savings  can  be  realized  in  day-to- 
day  operations. 
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Fuel-optimal  flight  conditions  for  the  subject  air- 
craft are  not  sharply  peaked,  and  major  fuel  savings  can  be 
accomplished  by  operating  in  the  vicinity  of  the  optimum 
condition.  Aerodynamic  substantiating  data  and  documented 
standard  procedures  provide  valuable  reference  points  for 
fuel  optimization,  although  certain  results  require  considera- 
tion of  the  specific  stage  length  and  vehicle  weight,  par- 
ticularly for  short-range  flights.  For  example,  the  optimum 
climb  schedule  for  short-range  flights  uses  a constant  Mach 
number  and  reduced  throttle  setting,  whereas  the  medium-  to 
long-range  flight  calls  for  constant  indicated  airspeed  and 
maximum  throttle  setting  during  the  climb. 

The  time  penalty  associated  with  fuel-optimal  flight 
paths  is  not  great,  and  the  decreased  fuel  costs  can  lead 
to  decreased  direct  operating  costs.  Except  for  long-range 
cruise  at  an  altitude  substantially  below  the  unconstrained 
optimum  cruise  point  ( in  which  case , the  fuel-optimal  cruise  speed 
is  sharply  reduced),  the  total  flight  times  of  optimal  and 
conventional  profiles  are  comparable.  While  it  is  difficult 
to  account  for  total  operating  costs  <5f  the  military  transport 
(since  several  Air  Force  Commands  share  the  costs),  a com- 
parison with  airline  operating  costs  suggests  that  the  C-141A's 
fuel  costs  now  exceed  time-based  costs  on  most  flights;  hence, 
reduced  fuel  use  would  result  in  reduced  operating  costs. 

The  development  of  fuel-minimizing  guidance  laws  can 
be  separated  into  two  parts:  generating  the  optimal  1 light 

path  and  issuing  commands  which  force  the  aircraft  to  follow 
tho  path  in  the  presence  of  disturbances.  Simplified  refer- 
ence profiles  can  be  established  from  the  precise  optimal 
flight  paths,  and  the  path-following  guidance  laws  are  readily 
formulated  using  linear-optimal  control  theory.  The  latter 
design  process  leads  to  a combined  throttle/energy  management 
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structure  of  very  simple  form.  Such  a system  can  be  synthesized 
for  flight  test  using  a small  digital  computer,  and  it  could 
be  incorporated  within  navigation  and  control  systems  now 
under  development. 

AlteA'natives  f'r  flight  testing  center  around  two 
concepts.  The  first  uses  a dedicated,  well-instrumented  air- 
craft flying  under  carefully  controlled  conditions.  The 
second  approach  is  to  conduct  "piggy-back"  experiments  during 
operational  missions.  Both  types  of  flight  test  should  be 
carried  out  to  maximize  results  with  minimum  cost. 


1.4  ORGANIZATION  OF  THE  REPORT 

This  report  treats  several  aspects  of  the  aircraft 
fuel  conservation  problem.  Chapter  2 reviews  aircraft 
selection,  mission  constraints,  and  prior  developments 
related  to  fuel  conservation.  Chapter  3 describes  the 
mathematical  basis  for  the  investigation.  Optimal  perfor- 
mance results  for  the  C-141A  are  given  in  Chapter  4.  and 
these  results  are  compared  with  conventional  flight  paths 
in  Chapter  5.  Fuel-optimal  guidance  algorithms  are  developed 
in  Chapter  6,  and  test  planning  is  discussed  in  Chapter  7. 

The  report  is  concluded  by  Chapter  8.  Mathematical  tools 
are  described  in  Appendices  A and  B. 
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AIRCRAFT  USE  FACTORS  AND  PROSPECTS 


1 1 . 

FOR  FUEL  CONSERVATION 


2.1  OVERVIEW 

With  an  annual  fuel  usage  in  the  vicinity  of  more 
than  4 billion  gallons  and  an  associated  bill  of  about  $2 
billion  (Ref.  24),  the  U.S.  Air  Force  has  strong  motivation 
to  consider  alternatives  for  conserving  fuel.  Projected 
cost  increases  will  make  fuel  conservation  even  more  impor- 
tant in  years  to  come.  In  the  long-term,  improved  aircraft 
configurations  and  propulsion  systems  can  be  expected  to 
reduce  fuel  use,  but  for  the  short-  and  mid-term,  it  is 
necessary  to  consider  improvements  in  the  maintenance, 
operations,  and  subsystems  of  existing  aircraft. 

Although  this  investigation  is  directed  at  improving 
the  technology  base  for  fuel-efficient  aircraft  operations, 
particular  attention  is  devoted  to  obtaining  results  which 
can  assist  the  Air  Force  in  reducing  aircraft  fuel  use  in 
the  near-term.  To  achieve  this  objective,  aircraft  use 
factors  must  be  considered  in  the  context  of  mission 
requirements,  airspace  limitations,  crew  workload,  and  the 
probability  that  fuel-saving  procedures  can  be  incorporated 
in  day-to-day  operations.  This  chapter  discusses  these  fac- 
tors 


There  are  three  topics  covered  here.  Section  2.2 
treats  the  selection  of  a subject  aircraft,  including  the 
basic  characteristics  of  four  aircraft  types  and  the  criteria 
for  selection.  Section  2.3  presents  mission  profiles  and 
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constraints  for  the  chosen  aircraft,  the  C-141A  transport. 
Section  2.4  reviews  energy  managenent  concepts  for  fuel- 
efficient  flight,  and  the  chapter  is  sumnarized  in  Section  2.5. 


2.2  AIRCRAFT  SELECTION 


To  assure  that  the  energy  management  development 
provides  usable  results  at  an  early  stage,  a single  aircraft 
type  must  be  defined  for  detailed  study.  Breakdowns  of  air- 
craft fuel  use  by  aircraft  class  and  type  assist  the  selec- 
tion of  this  study  aircraft.  Nearly  90  percent  of  aircraft 
fuel  used  by  the  Air  Force  is  used  by  four  aircraft  classes: 


• 

Cargo 

32.7% 

• 

Fighter /Reconnaissance 

25.4% 

• 

Bomber /Reconnaissance 

18.1% 

• 

Tanker 

11.9% 

Six  aircraft 

types  use  almost  70  percent  of 

• 

C-141 

14.9% 

• 

B-52 

14.2% 

• 

F-4 

14.1% 

KC-135 

11.2% 

# 

C-5 

7.2% 

• 

C-130 

6.9% 

1973) 


These  lists  indicate  that  over  half  the  fuel  is  used 
by  large  aircraft  (cargo,  bomber,  and  tanker),  and  four  air- 
craft types  consume  the  bulk  of  this  (C-141,  B-52,  KC-135,  and 
C-5).  These  aircraft  were  candidates  for  study  in  this  inves- 
tigation (Ref.  25). 

A useful  analysis  could  be  made  for  any  one  of  these 
aircraft,  and  the  four  types  are  sufficiently  similar  that 
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the  results  obtained  for  one  type  are  generally  applicable 
to  the  other  three.  In  addition  to  the  total  fuel  consumed 
by  each  fleet,  other  factors  affect  the  selection,  e.g.  , 
availability  of  aerodynamic  and  engine  data,  mission  sce- 
narios, operational  constraints,  and  Air  Force  Command  fly- 
ing policies.  These  aspects  are  summarized  below  for  each 
aircraft  type. 

2.2.1  The  Candidate  Aircraft  Types 

C-141A  - Over  275  of  these  aircraft  have  been 
delivered  to  the  Air  Force  for  use  in  long-range  strategic 
transport.  Maximum  range  with  31, 000- lb  payload  is  about 
5,600  nm,  and  normal  cruise  Mach  number  is  0.72  to  0.767. 
Maximum  loaded  weight  is  approximately  300,000  lb.  This 
aircraft  has  been  the  subject  of  past  throttle/energy 
management  studies  (Refs.  1 and  2);  hence,  there  is  a body 
of  performance  data  available  for  ready  comparison.  As  a 
major  element  of  the  Military  Aircraft  Command  (MAC),  poten- 
tial fuel  savings  with  this  aircraft  would  be  significant 
if  the  throttle/energy  management  techniques  developed  in 
this  program  were  implemented  at  the  squadron  level.  The 
aircraft  was  first  test  flown  in  1963,  and  its  low-bypass- 
ratio  turbofan  engine  shares  its  core  (TF33)  with  the  earlier 
vintage  B-52.  The  C-141A  is  in  the  same  weight  class  as  the 
KC-135,  both  of  which  are  lighter  than  the  C-5A  or  B-52; 
hence,  fuel  costs  for  testing  throttle/energy  management 
techniques  would  be  relatively  low. 

C-5A  - This  aircraft  is  more  than  twice  as  heavy 
as  the  C-141A,  and  it  carries  three  times  the  payload.  The 
number  of  C-5A's  in  MAC  inventory  is  one-third  the  number  of 
C-141A's,  and  the  potential  C-5A  fleet  fuel  saving  is 
approximately  one  half  that  of  the  C-141A  (subject  to  actual 
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use  (actors).  Cruise  characteristics  of  the  aircraft  are 
similar  to  the  C-141A,  although  policies  directed  at  increas- 
ing aircraft  fatigue  life  by  reducing  average  flight  loads 
provide  additional  constraints  on  this  aircraft.  The  C-5A 
possesses  the  most  modern  standard  avionics  suite  in  this 
list  of  aircraft,  and  a fuel  management  computer  is  included 
as  part  of  the  standard  equipment.  The  actual  costs  of  test- 
ing are  likely  to  be  higher  than  those  of  any  other  aircraft 
in  the  group.  The  moderate-bypass-ratio  turbofan  engine 
(TF39)  is  more  representative  of  current  and  future  technology 
than  the  engines  in  the  other  three  aircraft  types. 

KC-135  - Since  1956,  the  Air  Force  has  ordered  over 
400  of  this  aircraft  type.  It  was  initially  equipped  with 
the  J57  turbojet  engine,  although  a few  aircraft  are  equipped 
with  the  TF33  turbofan  engine.  The  original  version  has  a 
"normal”  cruising  speed  of  510  kt  , somewhat  faster  than  the 
C-141A  and  C-5A , and  it  is  the  lightest  aircraft  in  the  group 
(270,000-lb  maximum  loaded  weight).  The  potential  Air  Force- 
wide  fuel  savings  with  this  Strategic  Air  Command  (SAC)  air 
craft  (some  also  are  assigned  to  MAC)  must  be  evaluated  in 
1 ight  of  refueling  mission  constraints,  which  could  preclude 
flving  for  best  fuel  consumption  for  extended  periods  of  time. 

B-52  - The  Stratofortress  was  designed  over  20  years 
ago,  first  flying  in  1952.  Production  of  the  Air  Force's  711 
aircraft  of  this  type  was  completed  by  1962.  Early  versions 
have  a maximum  loaded  weight  of  400,000  lb  and  are  equipped 
with  eight  J57  turbojet  engines;  later  models  are  100,000-1) 
heavier  and  use  TF33  turbofan  engines.  Normal  cruise  speed 
of  the  aircraft  is  490  kt . The  initial  analysis  and 
development  for  fuel-minimizing  cruise  which  preceded  Refs. 

1 and  2 was  directed  at  the  turbojet-powered  B-52  (Ref.  26); 
hence,  data  for  comparison  are  available.  This  aircraft  is 
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SAC's  major  user  of  fuel,  and  It  is  likely  to  retain  this 
status  for  the  next  several  years;  thus,  squadron-level 
implementation  of  fuel-saving  flight  policies  could  have  sub- 
stantial effect.  Mission  constraints  could  be  more  restrictive 
for  this  aircraft  than  for  the  other  three,  although  fuel 
optimality  is  a consideration  in  achieving  strategic  flexibility. 

2.2.2  The  Selected  Aircraft 

There  are  several  reasons  for  choosing  the  C-141A. 

This  aircraft  type  constitutes  the  single  greatest  user  of 
jet  fuel  in  the  Air  Force,  and  virtually  all  of  the  C-141A 
aircraft  delivered  to  the  Air  Force  have  identical  aero- 
dynamic, thrust,  and  inertial  characteristics  (allowing  for 
nominal  aircraft-to-aircraft  variation);  hence,  the  results 
of  the  study  are  applicable  to  all  aircraft  in  the  C-141A 
fleet.  (By  comparison,  the  B-52  and  KC-135  f Leets  contain 
both  turbojet-  and  turbo fan-powered  aircraft,  and  results 
of  the  study  depend  on  engine  type;  therefore  direct  results 
would  be  less  widespread  than  for  the  C-141A).  The  C-141A 
has  been  the  subject  of  previous  fuel-use  studies,  and  all 
data  required  to  formulate  the  mathematical  model  are  avail- 
able. Most  of  the  aircraft  in  the  C-141A  fleet  are  approach- 
ing the  mid-point  of  their  service  life  (now  projected  at 
40,000  flight  hours);  therefore,  the  specific  energy  manage- 
ment techniques  developed  here  could  be  applied  over  an 
extended  period  of  time.  The  proposal  to  modify  the  entire 
C-141A  fleet  to  the  "stretched"  C-141B  configuration  would 
have  little  effect  on  this  result,  as  the  range  of  param- 
eter variations  studied  includes  the  C-141B  performance 
characteristics . 

The  use  of  a transport  aircraft  rather  than  a bomber 
or  tanker  type  provides  fewer  constraints  in  determining 


12 


fuel-optimal  flight  policies,  as  point-to-point  navigation 
rather  than  strategic  (or  tactical)  navigation  is  involved. 
This  further  increases  the  possibility  that  throttle/energy 
management  techniques  developed  in  this  contract  ultimately 
could  find  their  way  into  early  operational  use.  In  addition, 
fuel-efficient  transport  flight  sets  the  stage  for  later 
studies  incorporating  strategic  or  tactical  constraints. 


2.3  MISSION  DEFINITIONS  AND  CONSTRAINTS 

For  fuel-u&e  studies  and  flight  profile  planning, 
C-141A  missions  can  be  separated  into  three  classes  with 
associated  utilization  percentages  (Ref.  25): 


Training 

16% 

Exercise 

38% 

Airlift 

46% 

Excluding  aborted  or  diverted  missions,  touch-and-go  train- 
ing, air  drop  exercises,  and  other  low-altitude  missions, 
over  80  percent  of  the  flight  hours  can  be  expected  to  be 
logged  on  missions  which  reach  normal  cruisirg  altitudes 
(greater  than  15,000  ft).  These  missions  are  considered 
candidates  for  detailed  study.  The  average  duration  of 
missions  in  this  category  is  5 hours,  with  minimum  and 
maximum  values  of  1 and  10  hours,  respectively.  Corresponding 
ranges  are  2,000,  300,  and  4,200  nm.  Payloads  range  from 
negligible  to  24  tons,  with  an  average  value  of  10.4  tons. 

2.3.1  Air  Traffic  Control 

While  air  traffic  control  (ATC)  presents  one  of  the 
most  restrictive  constraints  on  aircraft  energy  management, 
it  should  be  recognized  as  a factor  which  can  be  altered 
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and  influenced  to  achieve  desirable  objectives.  The  need 
for  an  ATC  system  arises  from  possibly  conflicting  demands 
for  the  use  of  airspace,  as  shown  in  Fig.  1;  hence,  if  these 
demands  can  be  coordinated,  or  if  air  traffic  is  light,  ATC 
becomes  less  of  a constraint.  Conversely,  if  airspace  sur- 
veillance is  inadequate,  if  ATC  policies  are  particularly 
rigid  and  unyielding,  or  if  air  traffic  is  heavy,  then  the 
constraint  is  increased.  Designs  for  throttle/energy  manage- 
ment in  a variable  ATC  environment  must  be  sufficiently  flexi- 
ble as  to  take  best  advantage  of  the  prevailing  constraints 
(or  lack  there, J)  to  save  fuel. 


Figure  1 Overview  for  Constraint  Definition 
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The  ATC  system  imposes  "hard"  constraints  on 
position  and  velocity  in  a time-based  framework.  The  ATC 
system  effectively  establishes  a network  of  corridors  and 
time-slots  within  which  an  aircraft  can  travel.  Explicitly 
or  implicitly,  the  aircraft  is  assumed  to  be  constrained  to 
a "4_D"  volume  in  time  and  space,  surrounded  by  a region 
which  the  aircraft  must  not  enter.  Current  civil  ATC  rules 
make  this  assumption  implicitly,  restricting  aircraft  to 
predetermined  altitudes  and  air lanes;  future  systems  pre- 
dicated on  all  aircraft  having  area  navigation  ("R-Nav") 
capabilities  could  assign  allowable  "4-D"  volumes  explicitly 
in  a very  flexible  way. 

For  the  time  being,  a military  transport  operating 
without  special  priorities  in  the  civil  ATC  system  must 
fly  according  to  fairly  restrictive  ri’les.  For  example, 
over  the  continental  United  States,  cruise  altitude  must 
be  held  constant,  altitude  is  assigned  with  2,000-ft  spacing 
below  29,000  ft  and  with  4,000-ft  spacing  above.  Assigned 
altitude  (a  function  of  other  traffic)  may  not  bo  optimal 
(Refs.  27  to  29).  Indicated  airspeed  must  be  held  below 
250  kt  at  altitudes  under  10,000  ft.  The  throttle/energy 
management  technique  developed  here  handles  such  constraints 
in  an  efficient  and  direct  way. 

2.3.2  Atmospheric  Conditions 

Winds,  atmospheric  profiles  of  pressure  and  tem- 
perature, and  severe  weather  (Fig.  1)  all  must  be  considered 
in  establishing  optimal  flight  paths,  although  only  the  latter 
falls  into  the  category  of  a "hard"  constraint.  These  fac- 
tors regularly  play  a part  in  military  and  commercial  route 
and  altitude  preflight  planning  (Refs.  30  and  31),  which 
attempts  to  minimize  fuel  or  direct  operating  costs. 
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Avoiding  severe  weather  is  similar  to  collision  or 
obstacle  avoidance,  and  a similar  approach  can  be  applied 
in  both  cases.  The  more  difficult  problem  is  identifying 
the  potential  conflict,  which  is  beyond  the  scope  of  the 
study.  Wind  and  atmospheric  profile  variations  can  be 
accounted  for  in  the  basic  guidance  algorithms,  while 
guidance  to  avoid  severe  weather  can  be  treated  as  an 
altitude/route  determination  problem,  thus  fitting  into 
the  same  framework  as  ATC  constraints. 

2.3.3  Vehicle  Constraints 

The  following  chapters  illustrate  that  fuel-optimal 
flight  of  a jet  transport  does  not  call  for  abrupt  maneuvers 
and  takes  place  at  moderate  dynamic  pressure  (and,  there- 
fore, indicated  airspeed)  and  Mach  number.  Consequently, 
the  vehicle's  structural  limitations  do  not  have  any  bear- 
ing on  shaping  the  vertical  flight  profile  to  minimize  fuel 
use , and  the  flight  condition  does  not  exceed  low  dynamic 
pressure  buffet  limits.  Although  the  performance  ceiling 
of  a lightly  loaded  C-141A  is  in  the  vicinity  of  50,000 
ft  (Ref.  32),  cabin  pressurization  cannot  be  maintained 
at  adequate  levels  for  passengers  above  41,000  ft.  As  for 
most  jet -powered  aircraft,  the  ability  to  make  idle-thrust 
descents  also  is  limited  by  cabin  pressurization  require- 
ments. (In  order  to  maintain  adequate  cabin  pressure 
levels,  one  or  more  powerplants  must  have  sufficient  speed 
to  drive  the  appropriate  subsystems.)4  Cabin  deck  angle 
(pitch  Euler  angle)  during  fuel-optimal  climb  may  be  greater 
than  conventional  climb  angle,  but  the  vehicle  contains  no 


♦Pressurization  constraints  have  not  been  applied  in  the 
following  analyses,  as  they  could  be  revised  by  subsystem 
modification,  or,  for  cargo  flights,  by  the  air  crew  using 
breathing  equipment.  During  descent,  the  C-141A  maintains 
adequate  pressurization  at  idle  thrust. 
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inherent  limitation  which  would  prevent  the  higher  angle. 
Flight  paths  investigated  here  pertain  to  the  "clean"  air- 
craft configuration  (flaps,  spoilers,  and  landing  gear  are 
retracted  and  the  rear  cargo  door  is  closed),  so  restrictions 
on  the  deployment  of  these  devices  are  not  pertinent. 


2.4  ENERGY  MANAGEMENT  FOR  FUEL  CONSERVATION 

Fuel-conservative  flight  has  a. vays  been  a major 
goal  in  aeronautical  development,  so  the  recent  concern 
for  saving  fuel  is  an  intensification  of  prior  interests 
rather  than  a new  subject  for  study.  Numerous  aspects  of 
aircraft  operations  can  be  examined  for  possible  fuel 
savings,  but  the  area  which  is  addressed  here  is  energy 
management . i.e. , the  regulation  of  an  individual  aircraft's 
kinetic  and  potential  energy  so  that  a minimum  amount  of 
powerplant  fuel  (chemical  energy)  is  expended  to  achieve 
mission  objectives.  As  mentioned  in  the  Introduction  (Sec- 
tion 1 . 1 ),  much  prior  analytical  work  in  this  area  has  been 
directed  at  the  equilibrium  cruise  condition,  during  which 
specific  total  energy  is  approximately  constant,  and  at 
climbs  and  descents,  which  result  in  transition  to  higher 
or  lower  energy. 

Figure  2 illustrates  the  relationships  between 
specific  energy  and  basic  flight  variables,  expressed  in 
conventional  units.  The  plot  presents  contours  of  specific 
total  energy,  E,  and  indicated  airspeed  (in  knots),  KIAS, 


♦It  will  be  recalled  that  the  aircraft's  kinetic  energy  is 
mV2/2,  where  m is  aircraft  mass  and  V is  velocity.  Poten- 
tial energy  is  mgh,  where  g is  gravitational  acceleration 
and  h is  height.  The  aircraft's  total  energy  is  the  sum 
oi  the  two,  and  its  specific  total  energy  (or  energy  per 
unit  weight),  E,  is  V*/2g  + i. 
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Figure  2 Contours  of  Constant  Specific  Total 
Energy  (E)  and  Indicated  Airspeed, 
in  Knots  (KIAS).  Initial  and  Final 
Point  for  Climb-Cruise-Descent-Studies 

Indicated  by  "Or 

as  functions  of  altitude  and  true  Mach  number,  M.  (M  and 
KIAS  are  computed  using  the  1962  Standard  Atmosphere, 

Ref.  33).*  Note  that  E has  the  dimensions  of  altitude,  and, 
in  fact,  altitude  is  equivalent  to  specific  potential 
energy.  On  a typical  flight,  the  aircraft  climbs  from  the 


♦Indicated  airspeed  (IAS)  and  calibrated  airspeed  (CAS)  are 
assumed  to  be  identical.  IAS  is  computed  as  1.689(2q/pgL)- , 
where  q is  the  dynamic  pressure  (lb/ft2)  and  pgL  is  the  sea 
level  air  density  (slugs/ft^). 
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lower  left  to  the  upper  right  of  this  figure  and  back, 
with  commensurate  variations  in  specific  energy.  Although 
M is  held  constant  during  fuel-optimal  cruise  of  the  C-141A, 
the  figure  shows  that  specific  energy  increases  (due  to 
altitude  increase)  as  fuel  is  burned  and  the  aircraft  loses 
weight.  This  leads  to  a subtle  effect  on  specific  range 
(nautical  miles  travelled  per  pound  of  fuel  expended),  which 
is  described  in  Chapter  4. 


Figure  3 illustrates  the  general  pattern  for  climb 
to  and  descent  from  the  cruise  point.  Contours  of  constant 
specific  energy  change  per  unit  of  fuel  burned  (3E/3m)  can 
be  defined;  if  these  contours  are  considered  to  portray  a 
three-dimensional  convex  figure,  the  "ridge  line"  of  this 
contour  defines  a fuel-minimizing  climb  path.  Similarly, 
contours  of  constant  drag  force  define  a minimum-drag  (or 
maximum- range)  glide  path  for  descent  from  the  cruise  condi- 
tion. 

n-ixwr 


MINIMUM-FUEL  ENERGY  CLIMB 


MINIMUM-DRAG  ENERGY  DESCENT 


CONTOURS  OF 


CONTOURS  OF  CONSTANT 


Figure  3 Fuel-Optimal  Energy  Transitions 
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As  a result,  fuel-optimal  flight  profiles  tend  to 
look  like  Fig.  4,  which  shows  climb,  descent,  and  two  fuel- 
minlmizing  cruise  alternatives  — a Breguet  cruise-climb  and 
a decelerating,  constant-altitude  cruise.  The  latter  cruise 
is  less  efficient  than  the  cruise-climb,  but  it  is  more  effi- 
cient than  a const ant -Mach,  constant-altitude  cruise.  In 
later  chapters  it  will  be  seen  that  all  fuel-optimal  flight 
paths  have  the  general  shape  shown  in  Fig.  4,  including  short- 
range  trajectories  without  well-defined  cruise  segments.  Spe- 
cifically, climb  speeds  are  faster  than  descent  speeds,  and 
the  flight  profile  approaches  operational  boundaries  (defined 
by  buffet  onset  or  by  maximum  operating  velocity,  V^,  or 
Mach  number,  MMO)  only  in  the  vicinity  of  cruise. 


Ml  30M 


TRUE  MACH  NUMBER 

Figure  4 Fuel-Optimal  Climb-Cruise-Descent  Profile 
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Although  these  trends  in  the  vertical  flight  path 
are  known,  actual  flight  operations  may  depart  from  this 
type  of  profile.  Standard  operating  procedures  frequently 
call  for  the  descent-speed  profile  to  be  the  reverse  of  the 
climb-speed  profile,  and  many  operators  erroneously  use 
higher-speed  "barber-pole  descents"  (flying  along  the  MMq- 
Vjjo  boundary)  to  save  fuel.  Because  the  duration  of  the 
barber-pole  descent  is  short,  the  apparent  use  of  fuel  during 
descent  is  low,  and  the  aircraft  can  stay  at  cruising  altitude 
(where  specific  range  is  high)  longer.  Nevertheless,  the 
total  fuel  used  during  the  cruise  segment  and  descent  is 
greater  than  that  used  on  a fuel-optimal  trajectory, 
especially  if  spoilers  are  used  to  increase  descent  rate  by 
increasing  drag. 

Cruise  segments  often  are  flown  one  or  two  per- 
cent faster  than  the  speed  for  best  specific  range  for  two 
reasons.  The  first  reason  is  to  reduce  time-based  direct 
operating  costs,  although,  at  present  fuel  rates,  the  net 
cost  saving  may  be  negligible.  The  second  reason  is  to 
provide  speed/thrust  stability,  i.e.,  to  assure  that  speed 
excursions  due  to  atmospheric  disturbances  are  nulled 
inherently  without  throttle  adjustment  (Refs. 34  and  35).  The 
additional  fuel  used  in  manual  throttle  adjustment  can  more 
than  offset  the  penalty  associated  with  the  higher  speed,  and 
conventional  autothrottles  may  have  a similar  effect  in 
strong  turbulence.  This  suggests  that  improved  automatic 
systems  could  help  by  augmenting  speed  stability  without 
throttle  "hunting"  and  by  automatically  following  a fuel- 
optimal  altitude/speed  schedule. 

A number  of  general  suggestions  related  to  energy 
management  for  fuel  conservation  can  be  found  in  Ref.  36. 

In  addition  to  recommendations  regarding  routing,  fuel 


21 


reserves,  ATC  clearances,  and  center-of-gravity  location, 
there  are  the  following: 

• Consider  the  use  of  area  navigation 
equipment 

• Adjust  en  route  speed  control  to 
optimum  values 

• Cruise  at  optimum  altitude 

• Use  optimum  climb  schedules 

• Use  idle-thrust  descent 

• Use  lowest  practical  flap  settings 
during  landing 

• Consider  advanced  onboard  avionics 
for  real-time  selection  of  optimum 
flight  profiles. 

These  suggestions  indicate  directions  for  fuel  saving  which 
are  considered  practical  by  industry  representatives,  and 
they  make  note  of  the  special  importance  attached  to 
augmenting  navigation,  guidance,  and  control  systems  to 
assist  the  pilot  in  saving  fuel. 

Practical  implementation  of  fuel-saving  flight 
policies  for  an  individual  aircraft  requires  detailed 
information  regarding  the  actual  sta';e  of  the  aircraft's 
engines,  aerodynamic  efficiency,  and  loading.  There  are 
a number  of  reasons  why  handbook  data  on  a particular  Air- 
craft type  may  be  inadequate  for  realizing  fuel  savings  on 
an  individual  aircraft  of  that  type.  Aircraft  of  a par- 
ticular design  series  often  have  modifications  which  occur 
during  production.  Rigging  of  aircraft  control  surfaces 
and  of  engine  actuator  settings  changes  with  each  overhaul , 
while  aircraft  use  patterns  and  aging  processes  affect 
aerodynamic  efficiency  and  the  calibration  of  engine  instru- 
ments and  fuel  gauges.  Consequently,  fuel-use  efficiency 
for  a particular  aircraft  is  dependent  on  a real-time  method 
for  sensing  these  changes  and  estimating  their  effects. 


22 


This  study  addresses  crucial  elements  of  a throttle/ 
energy  management  system  which  can  optimize  the  flight  paths 
of  individual  aircraft,  namely  the  determination  of  fuel- 
optimal  flight  paths  for  the  C-141A  and  the  synthesis  of 
guidance  algorithms  for  following  these  paths.  This  develop- 
ment assumes  that  the  aircraft's  characteristics  are  known, 
and  it  demonstrates  sensitivities  to  variations  in  the 
important  parameters.  The  complete  system  should  include 
fuel-use  estimation  logic  to  provide  a continuing  calibration 
of  the  relationship  between  throttle  setting  and  the  thrust/ 
fuel-consumption  characteristics  of  each  engine,  and  it  would 
be  aided  by  improved  thrust  measuring  devices,  as  reported, 
for  example,  in  Ref.  37. 


2.5  CHAPTER  SUMMARY 

This  chapter  has  presented  the  process  by  which 
the  C-141A  was  chosen  as  the  study  aircraft,  the  mission 
characteristics  and  limitations  for  this  aircraft,  and 
basic  principles  of  energy  management  for  aircraft  fuel 
conservation.  It  was  decided  that  the  potential  fuel- 
savings  payoff  of  the  study  would  be  greatest  by  examining 
the  C-141A.  This  single  jet  transport  type  consumes 
15  percent  of  the  Air  Force's  aircraft  fuel,  and  over  80 
percent  of  its  missions  could  be  flown  with  fuel -minimizing 
procedures.  Prospects  for  flying  fuel-optimal  profiles 
on  a day-to-day  basis  are  enhanced  by  an  automatic  throttle' 
energy  management  system  which  follows  optimal  altitude/ 
speed  schedules,  which  augments  cruising  speed  stability 
without  undue  throttle  "hunting",  and  which  adapts  to 
variations  in  aircraft  characteristics. 
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Ill . 


MATHEMATICAL  FORMULATIONS  FOR 


OPTIMIZATION  AND  ANALYSIS 


3.1  OVERVIEW 

Trajectory  optimization  cons' sts  of  determining  a 
"best"  trajectory  to  accomplish  a desired  mission.  The  mis- 
sion is  defined  by  the  aircraft's  initial  conditions,  the  ter- 
minal conditions,  and  the  in-flight  constraints  imposed  upon 
the  trajectory.  To  be  "best",  or  optimal,  implies  that  a 
trajectory  satisfies  the  mission  requirements  and  also 
minimizes  some  cost  function  used  to  evaluate  trajectory 
performance.  For  the  present  purpose  of  minimizing  fuel 
use,  the  obvious  cost  function  is  the  amount  of  fuel  con- 
sumed in  flying  a given  mission. 

An  aircraft  trajectory  is  characterised  by  the 
time  histories  of  the  vehicle  states  (altitude,  velocity, 
flight  path  angle,  etc.)  and  the  associated  control  his- 
tories (e.g.,  angle-ol -attack  and  throttle  seiting). 

The  dynamic  relationships  between  the  control  and  state 
histories  are  defined  by  the  aircraft's  equations  of  motion. 

As  these  equations  must  be  satisfied  by  any  physically 
ralizable  trajectory,  they  are  included  in  the  constraints 
which  the  optimal  trajectory  must  satisfy.  The  aircraft 
equations  of  motion  which  are  used  for  the  current  fuel 
minimization  analysis  are  presented  in  Section  3.2. 

All  types  of  aircraft  obey  the  same  set  of  motion 
equations.  What  differs  from  one  aircraft  type  to  another 
is  not  the  functional  form  of  the  equations,  but  rather  the 
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magnitudes  of  the  effects  of  the  control  histories  upon  the 
trajectory  state  response.  For  example,  a particular  throttle 
setting  results  in  different  values  of  thrust  and  fuel  flow 
rate  from  orte  aircraft  to  another.  Similarly,  a given  value 
of  angle  of  attack  yields  values  of  lift  and  drag  which  are 
aircraft-dependent.  The  mathematical  model  which  is  used 
to  describe  the  C-141A  aircraft  is  reviewed  in  Section  3.3 
of  this  report  and  presented  in  detail  in  Appendix  B.  Com- 
bining these  models  with  the  aircraft  equations  of  motion 
fully  specifies  the  interrelationships  of  the  vehicle's 
control  and  state  histories. 

In  order  to  examine  trajectory  performance  and  to 
exercise  optimization  techniques,  a means  of  trajectory 
generation  must  be  available.  Given  the  aircraft's  equations 
of  motion  and  mathematical  model,  all  that  is  needed  to  gen- 
erate a trajectory  by  integration  techniques  is  some  means 
of  specifying  the  control  history.  Reference  trajectories 
are  computed  using  control  histories  which  represent  standard 
aircraft  operations,  as  described  in  Section  3.4.  The  present 
study  embodies  three  alternatives  for  guidance  command  gener- 
ation. These  three  techniques  (Open-loop,  ad  hoc,  and  linear- 
optimal  perturbation  guidance)  are  discussed  in  Section  3.4. 
Trajectory  generation  by  any  of  these  means  does  insure  a 
physically  realizable  trajectory.  As  long  as  the  guidance 
commands  are  maintained  within  an  allowable  range,  the 
vehicle  model  and  motion  equations  accurately  reflect  the 
vehicle's  state  history. 

Numerous  methods  exist  for  performing  trajectory 
optimization  (Refs.  38  to  40).  Parametric  optimization, 
variational  optimization  (using  first-  or  higher-order 
algorithms),  dynamic  programming,  and  quasi-linearization 
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algorithms  are  but  a few  of  the  commonly  used  approaches  for 
optimization.  Each  technique  has  its  inherent  advantages 
and  disadvantages.  For  this  reason,  two  optimization  tools 
have  been  applied  to  the  fuel-minimizaticn  problem  under 
consideration . 

Parametric  optimization  is  a straightforward  method 
for  finding  the  maximum  or  minimum  of  a function  which  is 
determined  by  a finite  set  of  parameters.  Partial  deriva- 
tives are  not  required,  and  the  integration  of  sensitivity 
equations  is  unnecessary,  thus  making  parametric  optimization 
the  simplest  of  the  optimization  methods.  A number  of  prob- 
lems concerned  with  fuel-minimization  can  be  formulated  as 
parameter  optimization  problems,  and  this  analysis  tool  has 
proved  to  be  useful  in  this  study . 

Variational  optimization,  while  more  complex  than 
the  parametric  method,  provides  the  capability  to  perform 
optimization  in  a more  general  framework.  A first-order 
gradient  algorithm  is  used  here,  as  it  has  demonstrated 
its  effectiveness  in  numerous  aerospace  trajectory  optimi- 
zation applications  (Refs.  41  to  45).  Higher-order  algo- 
rithms have  somewhat  better  convergence  characteristics  in 
some  test  cases  (Ref.  39),  but  their  performance  is  not 
uniformly  superior,  and  the  complexities  involved  in  imple- 
menting these  higher-order  techniques  are  not  justified  for 
the  current  problem. 

Descriptions  of  the  parametric  and  steepest -descent 
optimization  algorithms  are  presented  in  Sections  3.5  and  3.6, 
respectively.  These  descriptions  provide  only  the  details  of 
the  optimization  techniques  without  the  burden  of  theoretical 
derivation,  which  can  be  found  in  Ref.  39.  The  chapter  is 
summarized  in  Section  3.7. 
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3.2 


AIRCRAFT  EQUATIONS  OF  MOTION 


In  this  section,  the  "point-mass"  equations  of  motion 
applicable  to  the  aircraft  fuel  minimization  are  derived.* 
First,  the  three-degree-of -freedom  (3-DOF)  equations  are  pre- 
sented. The  two-degree-of-f reedom  (2-DOF)  equations  describ- 
ing motion  in  a vertical  plane  then  are  derived  as  a subset  of 
the  3-DOF  model.  These  2-DOF  equations  provide  the  basis  for 
all  fuel  minimization  analysis  performed  in  this  report.  A 
change  of  independent  variable  from  time  to  range  is  shown  o 
further  simplify  trajectory  optimization  analysis.  Finally, 
the  usefulness  of  a new  technique  for  reducing  the  computer 
time  necessary  for  trajectory  integration  is  demonstrated. 


The  3-DOF  aircraft  equations  of  motion,  assuming 
zero  sideslip  angle  and  neglecting  round,  rotating  earth 
effects,  can  be  written  as 


V 

Y 


* T(^nP')  COS  (ot+otT)  " f(V+Vw)2  CD  m ~g  sin  ' (3.2-1) 


^)Sin(a-aT)  * ift<V+Vw>2  CL  S 


cos 


cos  y 


(3.2-2) 


h = V sin  y 


(3.2-3) 


R = V cos  Y cos  5 (3.2-4) 

CR  = V cos  y sin  £ (3.2-5) 

l (a+otT^  + (V+Vw)2  CL  §]sin<J>/cos  Y (3.2-6) 

m = uj(T , NRT,  M,  h)  (3.2-7) 

where  the  dependent  variables  are 

♦"Point-mass"  equations  of  motion  describe  the  dynamics  of 
the  aircraft's  center  of  gravity,  and  they  assume  that  angu- 
lar maneuvers  can  be  accomplished  in  negligible  time. 
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V ■ Earth-relative  velocity,  fps 

Y * Flight  path  angle,  deg 
h * Altitude,  ft 

R * Range,  ft 

CR  ■ Crossrange,  ft 

£ * Horizontal  heading  angle,  deg 

m **  Aircraft  mass,  slugs 

to  * Fuel  mass  flow  rate,  slugs/sec 

The  model  parameters  include  engine  normal-rated 
thrust  (NRT),  throttle  setting  (T),  vehicle  lift  and  drag 
coefficients  (CL  and  Cp),  angle  of  attack  (a),  thrust-vector 
incidence  angle  (aT),  vehicle  reference  area  (S),  atmospheric 
density  (p),  gravitational  acceleration  (g),  and  vehicle  Mach 
number  (M).  The  term  Vw  represents  the  atmospheric  wind 
velocity  and  is  assumed  to  be  aligned  with  the  vehicle’s 
velocity  vector.  This  assumption  is  appropriate  for  the 
determination  of  horizontal-wind  effects  during  cruising 
flight  (y  - 0),  which  is  discussed  in  Section  4.2.4.  In 
the  preceding  equations,  total  engine  thrust  is  modeled  as 
the  product  of  normal-rated  thrust  (summed  for  each  opera- 
ting engine),  and  throttle  setting. 

The  2-D0F  equations  are  obtained  from  the  above  set 
by  constraining  roll  angle,  crossrange,  and  heading  to  be 
zero.  The  resulting  set  of  equations  can  be  expressed  as 

V - cos  (a+aT)  - X (V+Vw)2  C„  f -g  sin  i (3.2-8) 

T ’ T(^r)  <“+“T>  - % <V+Vw>2  CL  I -V  cos  * (3-2'9> 
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h ■ V sin  y 


(3.2-10) 


R - V cos  y (3.2-11) 

m = uj(T,NRT,M,h)  (3.2-12) 

Since  the  above  dynamic  equations  have  no  explicit 
dependence  on  time,  their  number  can  be  reduced  by  redefining 
the  independent  variable  to  be  one  or  a combination  of  the 
dependent  variables.  The  new  independent  variable  must  be 
monotonic  along  the  trajectory  to  avoid  singular  points  and 
multivalued  control  histories.  For  certain  applications,  spe- 
cific total  energy,  E,  is  an  appropriate  choice  of  independent 
variable  (Refs.  42  and  45).  During  aircraft  cruising  flight, 
however,  E,  is  only  slightly  changing  and  may  not  be  monotonic. 

For  the  present  analysis,  the  trajectory  path  length, 
R,  proves  to  be  a good  choice  of  independent  variable.  Not 
only  is  path  length  a monotonic  trajectory  variable,  but  the 
initial  and  final  points  of  an  aircraft  mission  actually 
define  total  path  length  for  planar  trajectories.  Use  of 
path  length  as  the  independent  variable  also  enables  the 
optimization  problem  to  be  formulated  in  the  "fixed-"  (rather 
than  "free-")  endpoint  framework.  Typically,  fixed-endpoint 
optimization  analysis  is  simpler  to  perform  because  the 
sensitivities  of  the  terminal  state  conditions  with  respect 
to  terminal  point  variations  need  not  he  considered  (Refs. 

38  and  39 ) . 

Conversion  of  the  dynamic  equations  to  utilize  path 
length  rather  than  time  as  the  independent  variable  is  accom- 
plished by  invoking  the  chain  rule  for  derivatives.  For 
example,  the  rate  of  change  of  velocity  with  respect  to  path 
length  is 
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dV  , /dy\/dt\  . dv/dR 

dR  \dt /\dR/  dt/dt  ^ 13) 

Equation  (3.2-11)  can  be  rearranged  to  yield  the  relationship 
between  time  and  path  length  as 

» l/( V cos  y)  (3.2-14) 

The  use  of  the  chain  rule  and  the  t ime/path-length 
relationship  of  Eq.  (3.2-14)  results  in  the  following  set  of 
2-DOF  motion  equations: 

V'  = [t(^)cos  (*+aT)  - § (V+Vw)2  CD  | - g sin  yJ/(V  cos  y ) 

(3.2-15) 

Y'-  sin  (a+aT)  + § <VfVw)2  CL  I " g COS  ^]/(v2  cosy) 


(3.2-16) 

h'=  tan  y 

(3.2-17) 

m'  = u)(T,NRT,M,h)/(V  cos  y) 

(3.2-18) 

where  ()'  is  used  to  represent  derivatives  with  respect  to 
path  length  rather  than  time. 

Although  angle  of  attack,  a,  appears  as  an  indepen- 
dent control  parameter  in  the  equations  of  motion,  the  use  of 
pitch  angle,  0 , as  a control  variable  is  desirable.  Refer- 
ence 46  demonstrates  that  the  use  of  0 in  place  of  a for  con- 
trol tends  to  damp  the  undesirable  phugoid  oscillations  which 
arise  with  constant  angle-of-attack  commands.  In  the  case  of 
the  planar  trajectory  equations,  the  relationship  between 
pitch  angle  and  angle  of  attack  is  given  simply  as 
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0 = a + y 


(3.2-19) 


The  equations  of  motion  derived  in  this  section  are 
driven  by  many  vehicle-dependent  effects  on  lift,  drag,  and 
thrust.  A description  of  these  vehicle  parameters  is  pre- 
sented in  the  following  section. 


3.3  C-141A  PARAMETERS 

The  aircraft  equations  of  motion  presented  in  Section 
3.2  are  valid  for  all  types  of  aircraft.  The  study  of  a 
particular  aircraft  necessitates  the  specification  of  those 
parameters  driving  the  equations  of  motion  which  are  vehicle- 
dependent.  The  parameters  which  must  be  defined  are  the  air- 
craft's lift  and  drag  coefficients,  CL  and  Cp,  reference 
area,  S,  normal-rated  thrust  from  the  engines,  NRT,  and  fuel 
consumption  rate,  u>. 

A complete  description  of  the  mathematical  model  of 
the  C-141A  is  presented  in  Appendix  B.  This  model  includes 
all  the  numerical  values  of  the  vehicle-dependent  parameters. 
Only  the  functional  relationships  of  these  parameters  are  given 
here,  in  order  to  review  those  quantities  which  affect  tra- 
jectory performance. 

The  complex  functional  interrelationships  of  the 
C- 141A  vehicle  and  atmospheric  parameters  are  summarized  by 
the  following  equations: 

CL  - f (a , h,  M,  q,  NRT,  T)  (3.3-1) 

CD  ■ f(CL,  M)  (3.3-2) 

NRT  = f (h,  M,  TEMP)  (3.3-3) 
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U) 


= f ( NRT , T,  6,  TEMP) 


(3.3-4) 


M = f(V, 

TEMP) 

(3.3-5) 

q =f(V, 

P) 

(3.3-6) 

TEMP  = f ( h) 

(3.3-7) 

P =f(h) 

(3.3-8) 

6 * f(h) 

(3.3-9) 

The  atmospheric  model  parameters  include  total  temperature 
(TEMP),  density  (p),  and  pressure  ratio  (6).  The  vehicle's 
dynamic  pressure  (q)  and  Mach  number  (M)  are  computed  from 
both  the  dependent  variables  of  the  equations  of  motion  and 
the  environmentally  determined  (atmospheric)  quantities. 

All  of  the  above  expressions  are  highly  nonlinear 
functions  of  the  model  parameters.  The  environmental  models 
employed  in  the  current  study  are  those  of  the  1962  Standard 
Atmosphere  (Ref.  33).  A number  of  sources  (Refs.  1,  2,  32, 
47,  and  48)  are  used  to  define  the  C-141A  model.  These 
sources  are  chosen  to  provide  a complete  and  valid  mathe- 
matical vehicle  representation.  The  fuel-minimization  study 
accounts  for  all  the  major  factors  which  influence  vehicle 
performance. 


3.4  REFERENCE  TRAJECTORY  COMPUTATION 

Generation  of  aircraft  trajectories  is  an  important 
part  of  fuel-minimization  analysis.  The  vehicle  equations  of 
motion  and  the  mathematical  model  of  the  C-141A  are  two  of  the 
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three  components  necessary  for  trajectory  computation.  Addi- 
tion of  a guidance  algorithm  enables  the  generation  of  com- 
plete trajectory  histories  by  means  of  numerical  integration. 
The  resulting  profiles  are  also  physically  realizable  as  long 
as  the  guidance  commands  are  maintained  within  allowable 
bounds.  (In  the  current  context,  the  terms  "guidance"  and 
"control"  are  used  interchangeably;  in  a complete  on-board 
system,  the  "control  system"  would  produce  the  forces  and 
moments  required  to  follow  the  "guidance  commands."). 

Both  open-loop  and  closed-loop  guidance  techniques 
may  be  used  to  compute  aircraft  controls.  Open-loop  guidance 
systems  specify  control  histories  as  functions  of  an  inde- 
pendent dynamic  variable,  while  closed-loop  commands  are  based 
upon  values  of  the  dependent  variables  of  the  equations  of 
motion.  In  the  present  study,  control  commands  (pitch  angle 
and  throttle  setting)  which  are  specified  as  functions  of 
path  length  (the  independent  variable  in  the  motion  equations) 
provide  open-loop  guidance.  Controls  computed  from  dependent 
trajectory  variables,  such  as  altitude  and  velocity,  form  a 
closed-loop  guidance  law. 

Open-loop  guidance  algorithms  have  only  limited  use- 
fulness for  actual  vehicle  implementation.  The  effects  of 
winds  and  other  uncertain  conditions  preclude  specifying  guid- 
ance commands  in  a totally  open-loop  fashion.  Mathematical 
trajectory  optimization  techniques  do,  however,  employ  this 
form  of  control  specification.  The  variational  optimization 
analysis  described  in  Section  3.6  demonstrates  the  use  of 
open-loop  guidance  systems  for  analytical  trajectory  shaping. 

Closed-loop  guidance  algorithms  can  be  employed  for 
both  mathematical  vehicle  studies  and  direct  aircraft  con- 
trol. A perturbation  guidance  system,  which  generates  control 


commands  based  upon  the  difference  between  desired  and 
actual  flight  conditions,  is  a form  of  closed-loop  guidance. 
The  use  of  this  guidance  technique  for  aircraft  implementa- 
tion (linear-optimal  perturbation  guidance)  is  examined  in 
Section  6.2.  Simple  ad  hoc  perturbation  techniques  serve 
the  purposes  of  the  parametric  optimization  efforts  under- 
taken in  this  study. 

The  ad  hoc  guidance  command  generation  used  for 
parametric  optimization  in  the  current  analysis  is  described 
by  the  following  equations: 

8 - % <vv> + Ke2  <V*> 

T = Kt  (Ejj-E) 

where  0 and  T are  the  time  derivatives  of  the  commanded 
pitch  angle  and  throttle  setting,  V and  E are  the  aircraft's 
velocity  and  specific  total  energy,  respectively,  and  the  sub- 
script D denotes  desired  values  of  these  quantities.  The 
derivative  of  desired  velocity,  Vp , is  computed  by  dividing 
the  change  in  VD  over  each  computation  interval  by  the 
interval's  time  increment. 

The  specific  energy  term  in  Eq . (3.4-2)  is  a func- 
tion of  the  vehicle's  velocity  and  altitude, 

E - h + — ■ (3.4-3) 

where  g is  the  earth's  gravitational  acceleration.  The  use 
of  energy,  rather  than  altitude,  in  the  throttle  command 
equation  reduces  the  effects  of  phugoid  oscillation  (the  in- 
terchange of  kinetic  and  potential  energy)  on  the  trajectory 
response. 


(3.4-1) 

(3.4-2) 


The  desired  values  in  Eqs.  (3.4-1)  and  (3.4-2)  may 
themselves  be  derived  quantities.  For  example,  it  may  be 
desired  to  fly  a constant-altitude  (hjj)/constant-Mach  number 
(Mp)  trajectory.  In  this  case,  desired  values  of  velocity 
and  energy  would  be  computed  as 


V 


D 


E 


D 


(3.4-4) 

(3.4-5) 


where  a is  the  speed  of  sound. 


The  vehicle  equations  of  motion  derived  in  Section 
3.2  belong  to  a class  of  differential  equations  commonly 
called  "stiff"  (Refs.  49  and  50).  This  class  of  differential 
equations  is  characterized  by  the  fact  that  large  differences 
exist  between  the  response  times  of  the  individual  equations. 
In  the  current  study,  the  flight  path  angle  equation  is  much 
faster  than  the  velocity,  altitude,  and  mass  equations. 
Typically,  this  necessitates  using  extremely  sma.l  step 
sizes  in  the  numerical  integration  procedure  to  accurately 
determine  the  trajectory  characteristics. 


In  the  current  efforts,  a simple  time-scaling  pro- 
cedure has  been  employed  to  enable  increased  integration 
intervals.  The  time  constant  of  the  flight  path  angle  differ- 
ential equation  is  modified  to  be  more  compatible  with  that 
of  the  velocity  equation.  This  is  accomplished  by  using  a 
suppression  factor,  e,  to  lower  the  value  of  y'  obtained 
from  Eq.  (3.2-16).  Note  that  in  the  1 ‘.mit , e = 1 , and  the 
flight  path  angle  equation  is  exact;  as  c approaches  0,  the 
dynamic  effect  of  flight  path  angle  change  becomes  negligible, 
and  the  equation  becomes  an  algebraic  constraint.  This  con- 
straint frequently  is  employed  in  approximate  trajectory 
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optimization  methods,  e.g.,  the  energy-state  method  mentioned 
in  Section  1.1.  To  insure  that  this  suppression  factor  does 
not  substantially  alter  the  resultant  trajectories,  flight 
paths  generated  both  with  and  without  the  suppression  factor 
are  compared.  Altitude  profiles  for  suppressed  (r.  ■ 0.05)  and 
unsuppressed  (e  = 1.0)  sets  of  motion  equations  are  pre- 
sented in  Fig.  5.  The  profiles  are  generated  employing  the 
ad  hoc  guidance  algorithm  demonstrating  a climb  in  altitude 
from  34,000  to  35,000  ft.  The  small  differences  in  overall 
profiles  ensure  that  the  suppressed  equation  set  provides  an 
excellent  representation  of  the  vehicle  dynamics  while  allow- 
ing the  computer  time  needed  for  trajectory  integration  to  be 
reduced  by  a factor  of  nearly  five.  The  unsuppressed  tra- 
jectory settles  at  the  desired  altitude  more  quickly  and 


Figure  5 Altitude  Profiles  for  Nominal  and 
Suppressed  Equations  of  Motion 
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with  less  overshoot  than  the  suppressed  version  does  because 
the  unsuppressed  motion  equations  have  a faster  flight  path 
angle  response. 


3.5  PARAMETRIC  OPTIMIZATION  OF  FLIGHT  PATHS 


The  basic  concepts  of  parametric  optimization  and 
its  applications  in  determining  fuel-optimal  aircraft  tra- 
jectories are  reviewed  in  this  section.  Parametric  optimi- 
zation's major  advantage  over  other  optimization  methods 
lies  in  the  fact  that  it  requires  only  funct ion  evaluations, 
not  partial  derivative  or  gradient  evaluations.  It  tends  to 
be  a computationally  efficient  technique  (Ref.  51).  Of 
course,  the  use  of  this  optimization  tool  is  limited  to  that 
class  of  problems  in  which  the  cost  function  can  be  minimized 
by  a finite  set  of  parameters.  Furthermore,  if  the  initial 
choice  of  parameters  is  inappropriate  or  inadequate,  further 
improvements  over  the  parametric  optimum  are  possible. 

Relative  comparisons  of  trajectory  performance  are 
made  possible  by  assigning  a scalar  cost  value  to  each  tra- 
jectory. Trajectory  A is  considered  to  be  superior  to  tra- 
jectory B if  the  cost  associated  with  A is  lower  than  that 
of  B.  An  optimal  profile  is  one  which  provides  the  lowest 
possible  cost . 

For  the  fuel  minimization  study  at  hand,  the  scalar  cost 
function,  J,  used  to  characterize  trajectory  cost  has  the  form 


J 


where 


*(x(Rf))  + 


L(x,  u,  R)  dR 


(3.5-1) 
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total  scalar  cost  function 


x(Rf ) 


R0’Rf 


L 

* 


final  values  of  dependent  trajectory 
state  variables:  velocity,  (V), 

altitude  (h),  flight  path  angle  ( y ) , 
mass  (m)  and  time  (t) 

vector  of  control  variables:  pitch 

angle  (0)  and  throttle  setting  (T) 

initial  and  final  values  of  the 
independent  variable:  range 

scalar  in-flight  cost  function 

scalar  cost  function  of  final  stat o variables 


The  terminal  cost  function,  ♦,  is  used  to  penalize  the  tra- 
jectory for  not  meeting  desired  terminal  conditions,  such  as 
final  altitude  and  velocity,  and  to  provide  a direct  penalty 
for  the  amount  of  fuel  consumed.  The  in-flight  cost  function. 

L,  penalizes  the  trajectory  for  violating  desired  flight  path 
quantities,  such  as  exceeding  a maximum  altitude  or  Mach  number. 
A complete  list  of  penalty  functions  available  for  use  in  the 
C-141A  fuel  minimization  analysis  is  presented  in  Appendix  A. 2. 


As  previously  noted,  parametric  optimization  deter- 
mines the  lowest  cost,  J*,  of  a function  which  can  be  de- 
scribed by  a finite  set  of  parameters,  i.e., 

J*  = min  J(p1,  P2  ...  PN)  (3.5-2) 

where  the  p^  are  the  parameters  which  can  be  varied.  The  re- 
sults of  the  optimization  procedure  are  J* , the  optimal  cost, 
and  the  p^*(i=l  to  N).  which  are  the  parameter  values  that 
yield  the  lowest  cost. 

The  ad  hoc  closed- loop  guidance  algorithm  described 
in  Section  3.4  provides  the  means  by  which  the  fuel-minimiza- 
tion analysis  can  be  cast  in  parametric  form.  The  parameters 
which  can  be  varied  are  those  quantities  which  define  the 
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desired  trajectory  characteristics  (velocity.  V^,  and  energy 
Kp)  and  are  inputs  of  the  ad  hoc  guidance  equations.  The 
parametric  optimization  provides  the  optimum  values  of  thest 
parameters  for  particular  trajectory  segments. 

Two  examples  of  parametric  optimization  are  illus- 
trated in  Fig.  6.  In  the  first  case  (Fig.  6a),  a two- 
segment  profile  which ‘ minimizes  fuel  used  during  a climb  is 
generated,  with  the  profile  terminating  at  an  optimal  cruise 
condition  (specified  by  weight  over  atmospheric  pressure 
ratio,  W/S.  and  Mach  number,  M).  The  two  parameters  which 
are  optimized  are  1)  the  value  of  indicated  airspeed,  IAS. 
during  the  initial  climb  phase,  and  2)  tin*  Mach  number 
during  the  second  portion  of  the  climb.  The  results  of  the 
optimization  are  the  best  values  of  IAS  and  M to  fly  if  this 
t "po  of  profile  is  to  be  implemented  during  the  climb. 


• ? SEGMENT  CLIMU  PHOFII.E  AT  MAXIMUM 
THROTTLE  SETTING 

• MINIMUM  FUEL  USE  USING  CONVENTIONAL 
PROCI  OURES 


• CONVENTIONAL  PHOCtOURI  S THROUGHOUT 
ENTIRF  PROFILE 

• CONSTRAINTS  to  9 , conttjm  crime  -Hitmlcl  WINDS. 
AND  VEHICLE  VARIATIONS  INCOHI'OR  ATE  0 


a)  CLIMB  OPTIMIZATION 


bl  DESCENT  POINT  OPTIMIZATION 


Figure  6 


Examples  of  Parametric  Optimization 


The  second  example  (Fig.  6b)  demonstrates  the  use 
of  single-parameter  optimization.  In  this  case,  the  point 
for  initiating  descent  is  chosen  to  match  the  desired  values 
of  terminal  altitude  and  velocity.  Thus,  complete*  near- 
optimal  trajectory  profiles  can  be  generated  by  appropriate 
use  of  parametric  optimization  methods. 


The  parametric  optimization  procedure  determines 
the  best  parameter  values  by  successive  function  evaluations. 
A single  parameter  case,  Fig.  7a,  serves  as  a simple  example 
of  the  optimization  concepts.  As  shown  in  the  figure,  the 
cost  function  of  parameter  p,  J(p),  is  evaluated  three  times: 
for  values  of  p equal  to  x,  (x-x),  and  (x  + x).  A parabola 
is  then  fitted  to  the  three  values  of  cost  associated  with 
these  parameter  values,  and  a final  function  evaluation  is 
porformea  at  the  predicted  minimum  of  the  parabola.  If  the 
cost  function  were  truly  parabolic  in  the  parameter,  p,  only 
one  such  iteration  would  serve  to  locate  the  true  minimum. 

For  most  situations,  however,  the  cost  function  is  not 
exactly  parabolic,  and  the  sequence  of  function  evaluations 
must  be  repeated  several  times  to  approach  the  optimum  point 
as  closely  as  desired. 


o)  SIMPLE  parameter  minimization 
Figure  7 Parametric  Minimization  Concepts 
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A two-dimensional  parametric  optimization  is  de- 
picted in  Fig.  7b.  The  basic  concepts  are  the  same  as  for 
tho  single  parameter  case,  but  more  function  evaluations  are 
required.  The  first  set  of  function  evaluations,  denoted  by 
points  1 to  4 in  the  figure,  represents  the  three  trial  points 
(1-3)  and  the  predicted  minimum  (4)  in  the  p^  - Pg  parameter 
space.  Function  evaluations  5 and  6 are  again  trial  steps; 
however,  they  are  performed  along  a direction  which  is 
orthogonal  to  the  initial  search  direction.  Function  evalua- 
tion 7 is  performed  at  the  predicted  minimum  of  this  new 
search  direction. 

As  in  the  single  parameter  optimization,  successive 
iterations  enable  the  optimum  point  to  be  approached  as 
closely  as  desired.  The  program  which  is  used  for  the  para- 
metric optimization  procedure,  EXTREM,  is  a modified  version 
of  a program  presented  in  Ref.  40.  This  program  also  modifies 
the  magnitude  of  the  search  steps  to  improve  convergence,  and 
it  is  capable  of  handling  parametric  optimizations  of  any 
dimension . 

Parametric  techniques  have  been  shown  to  have  the 
capability  to  define  optimal  segments  and  complete  trajec- 
tories. These  trajectories  are  optimal  in  that  they  are 
the  best  trajectories  which  follow  the  profiles  determined 
by  the  specified  parameters.  Variational  optimization 
methods  determine  optimal  trajectories  which  are  not  confined 
to  pre-specif ied  profiles  (e.g.,  constant  IAS  and  Mach 
number  climb  modes).  The  parametric-optimal  flight  paths, 
however,  facilitate  the  use  of  variational  optimization 
methods  by  providing  a good  starting  trajectory  for  the  exact 
optimal  solution. 
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3.6 


VARIATIONAL  OPTIMIZATION  OF  FLIGHT  PATHS 


Variational  optimization  methods  determine  the  "best'' 
trajectory  profiles  subject  to  any  number  of  desired  terminal 
conditions  and  in-flight  constraints.  Although  variational 
techniques  are  more  difficult  to  implement  than  parametric 
algorithms,  they  are  applicable  in  those  cases  where  the 
proper  choice  of  parameters  is  not  obvious.  Thus,  the  re- 
sultant optimal  profiles  are  not  constrained  to  follow  a 
parametrically  defined  (e.g.,  constant  Mach  number ) prof i le  . 

The  particular  variational  technique  employed  here 
combines  a steepest-descent  algorithm  with  the  direct  search 
used  for  parametric  optimization.  It  is  a first-order  gra- 
dient algorithm  which  iteratively  modifies  the  aircraft  con- 
trol history  to  improve  overall  trajectory  performance.  The 
control  perturbations  are  chosen  to  maximize  the  reduction 
in  trajectory  cost  per  unit  change  in  control  history,  i.e., 
the  controls  are  applied  along  projections  of  the  gradient. 
Thus,  the  technique  used  here  is  a modified  steepest-descent 
approach  which  retains  the  essential  features  of  many 
accelerated  gradient  methods. 


The  steepest-descent  gradient  algorithm  is  derived 
by  specifying  a cost  function  similar  to  that  of  Eq.  (3.5-1): 


J = <t>(x(Rf)) 


L(x,u,R) 


(f (x  ,u,R) 


(3.6-1) 


where  X is  the  costate  (or  adjoint)  vector  of  the  state,  x. 
Because  the  state  derivatives  with  respect  to  path  length 
are  given  by 
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(3.6-2) 


x ' = f(x,u,R) 

the  magnitude  of  the  second  term  under  the  integral  in 
Eq . (3.6-1)  is  zero.  Adjoining  the  dynamic  equations  to  the 
cost  function,  however,  enables  consideration  of  the  effects 
of  state  perturbations  at  R^,  caused  by  control  perturbations 
applied  at  previous  values  of  the  independent  variable  (Rq> 
for  example),  on  the  trajectory  cost.  The  relationship  be- 
tween the  state  ard  costate  variables  is  defined  by  the 
following  differential  equation: 

X’(R)  = - f^(R)  A(R)  - L^(R)  (3.6-3) 

T 

Here,  the  matrix  f is  the  transpose  of  the  partial  deriva- 
tive  of  the  vector,  with  respect  to  x.  The  vector,  Ll  , is 

tiie  partial  derivative  of  L (a  scalar)  witli  respect  to  x. 

This  differential  equation  is  integrated  backwards  from  Rf 
to  Rqi  as  only  the  terminal  values  A_  are  immediately  avail- 
able. The  terminal  values  are  defined  by  the  sensitivity 
of  the  terminal  cost  to  state  variations: 

A(Rf)  = 9 4>(x(Rj. ) )/9x  (3.6-4) 

The  Hamiltonian,  H,  is  a scalar  function  defined 
as 

H(x,u,X^R)  = L(x,u,R)  + xTf(x,u,R)  (3.6-5) 

It  is  shown  in  Ref.  39  that  the  Hamiltonian  has  the  prop- 
erty that  its  partial  derivative  with  respc?et  to  control 
perturbations,  H , is  zero  at  locally  stationary  points. 

This  is  a necessary  condition  for  a trajectory  to  be  optimum. 
When  H f 0,  it  indicates  which  way  to  adjust  the  controls 
to  reduce  total  cost. 
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The  gradient,  Hu,  at  a non-optimum  point  indicates 
the  direction  of  control  perturbation  which  effects  a reduc- 
tion in  total  cost.  It  can  be  computed,  using  the  costate 
vector,  as 

HU(R)=  -[lu(x.H.R)  + (x,u,R)  XT(R)J  (3.6-6) 

where  the  vector  Lu  is  the  partial  derivative  of  L with  re- 
spect to  the  controls  u,  and  is  the  partial  derivative 
matrix  of  f with  respect  to  u. 


Having  determined  the  direction  of  control  perturba- 
tion, actual  perturbations  can  be  generated  as 

6u(R)  - C HU(R)  (3.6-7) 

where  the  matrix  C contains  unspecified  gain  values.  In 
order  to  accelerate  the  convergence  of  the  gradient  algorithm 
to  the  optimal  trajectory,  the  present  study  makes  use  of  the 
parametric  optimization  tools  described  in  Section  3.5.  Th  s 
is  accomplished  by  specifying  a diagonal  form  for  C,  with  the 
diagonal  elements  being  independent  parameters.  For  a two- 
dimensional  case,  the  gain  matrix  would  be  defined  as 


C 


(3.6-8) 


and  the  resultant  control  perturbations  are  given  by 


f6u1(R)" 

m 

'Cl  HUl(B)' 

i 

os 

CJ 

3 

«o 

1 

c2  HU2CR)- 

(3.6-9) 


where  the  terms  5u^,  fiUg,  Hu^,  are  the  components  of 
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the  control  perturbation  and  gradient  vectors  respectively. 

The  new  control  history,  Hi+i*  *s  computed  as 

ui+i<R)  - ui(R)  + 6u(R)  (3.6-10) 

where  u^  is  the  previous  control  history  and  the  control 
perturbations  are  those  specified  in  Eq.  (3.6-9).  The 
iteration  is  complete  when  no  further  significant  reduction 
of  the  cost  function  can  be  obtained. 

The  gain  values,  c^  and  Cg.  can  be  chosen  through 
the  parametric  optimization  process.  Once  the  gradient,  Hu> 
is  determined,  the  changes  in  the  cost  function  are  only  ~ 
dependent  upon  the  choice  of  gain  parameters.  This  use  of 
parametric  optimization  in  the  variational  technique  enhances 
the  convergence  characteristics  of  the  first-order  algorithm, 
thereby  permitting  the  optimal  trajectory  to  be  determined 
with  fewer  steepest-descent  iterations.  As  the  computation 
of  trajectory  partials  (L.*,  L , f^,  f^)  can  be  burdensome, 
the  reduction  of-  iterations  results  in  a substantial  reduc- 
tion of  computer  time. 

The  execution  of  a steepest-descent  optimization 
iteration  consists  of  performing  a series  of  mathematical 
operations.  In  order  to  clarify  the  technique,  this  section 
is  concluded  with  an  ordered  list  of  these  operations: 

• Generate  an  initial  trajectory,  xq(R), 
given  an  arbitrary  initial  control 
history,  Uq(R),  using  open-loop  guidance 
commands . 

• Compute  trajectory  and  cost  function 

partial  derivatives:  L • 

• Integrate  (backwards)  the  costate  differ- 
ential equation  to  determine  A^R). 
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• Evaluate  the  gradient,  Hn(R),  in  order 
to  determine  control  perturbations, 

6u(R) . 

• Parametrically  optimize  the  control  per- 
turbation gains,  C,  to  obtain  the  lowest 
trajectory  cost,  J. 

• Either  terminate  execution  or  repeat,  starting 
at  the  second  step  and  using  the  best  trajec- 
tory as  the  new  initial  point. 


This  section  has  presented  the  basic  concepts  of 
variational  optimization.  The  results  of  applying  these  con- 
cepts to  the  fuel-minimization  problem  are  shown  in  the 
following  chapters  of  this  report. 


3 . 7 CHAPTER  SUMMARY 

The  optimization  and  trajectory  analysis  techniques 
brought  to  bear  on  the  C-141A  fuel-minimization  problem  have 
b<en  presented  in  this  chapter.  The  equations  of  motion 
applicable  to  the  current  study  and  the  mathematical  aircraft 
model  serve  to  define  trajectory  state  response  as  a function 
of  the  input  controls.  By  working  within  the  vehicle  model 
framework,  the  optimal  trajectories  generated  are  physically 
realizable,  i.e.,  they  can  be  flown  by  actual  aircraft. 

The  key  points  and  significant  advantages  of  the 
current  trajectory  analysis  can  be  summarized  as  follows: 

• Two-degree-of -freedom  equations  are 
employed  to  simplify  fuel-minimization 
analysis  in  the  vertical  plane. 

• The  C -141A  vehicle  model  includes  all 
the  major  factors  which  influence  air- 
craft trajectory  characteristics. 
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• A scalar  cost  is  assigned  to  each  tra- 
jectory. This  enables  the  relative 
comparisons  of  trajectory  performance. 

• Parametric  techniques  provide  for  com- 
putationally efficient  trajectory 
segment  optimization. 

• Steepest-descent  optimization  methods  are 
used  to  generate  fully  optimal  flight  pro- 
files and  vehicle  control  histories.  Para- 
metric methods  are  employed  to  enhance  the 
convergence  characteristics  of  the  steepest- 
descent  algorithm. 


The  analysis  tools  described  in  this  section  provide 
flexible  and  accurate  means  for  trajectory  shaping  and 
analysis.  The  results  of  applying  these  tools  to  the  fuel- 
minimization  problem  are  presented  in  the  following  chapters. 
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IV. 


AIRCRAFT  PERFORMANCE  ON  FLIGHT  PATH  SEGMENTS 


4.1  OVERVIEW 

In  this  chapter,  C-141A  performance  on  individual 
flight  path  segments  is  considered.  The  flight  profile  is 
divided  into  the  three  natural  flight  modes:  climb,  cruise, 
and  descent.  The  climb  segment  typically  is  flown  with 
maximum  throttle  setting  (Refs.  32  and  47),  and  during  this 
phase,  the  aircraft  increases  altitude  and  velocity  until 
cruising  flight  conditions  are  achieved.  Cruising  flight 
is  characterized  by  a constant,  or  slowly  varying, altitude 
and  velocity  profile.  The  descent  mode  is  initiated 
following  cruise,  and  it  serves  the  purpose  of  reducing 
altitude  from  the  cruise  to  landing  phases  of  flight. 

By  considering  these  flight  modes  individually, 
simple  analysis  tools  can  be  employed  to  optimize  fuel 
consumption.  The  parametric  methods  described  in  Section 
3.5  are  applicable  for  optimizing  the  climb  and  descent 
flight  segments,  while  direct  search  techniques  are  useful 
in  determining  the  conditions  which  define  fuel-optimal 
cruising  flight. 

Minimizing  fuel  consumption  on  the  individual  flight 
segments  provides  results  which  either  can  be  implemented 
directly  or  can  be  used  to  assist  in  complete  trajectory  op- 
timization. If,  for  example,  the  climb  and  descent  flight 
modes  are  constrained  (due  to  ATC  considerations),  at  least 
the  cruise  phase  may  be  flown  in  a fuel-minimal  fashion.  If, 
on  the  other  hand,  the  cruising  flight  segment  is  constrained, 
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some  fuel  savings  may  still  be  achieved  through  optimization 
of  the  climb  and  descent  schedules. 

By  piecing  together  the  individually  optimized 
segments,  fuel-efficient  complete  trajectory  profiles  are 
defined.  These  complete  trajectory  profiles,  however,  may 
not  be  fully  optimal  for  two  reasons.  First,  the  parametric 
optimization  methods  constrain  the  aircraft  to  fly  a specific 
profile  shape,  which  is  defined  by  the  parameters  to  be  opti- 
mized. Second,  by  considering  each  flight  mode  separately, 
the  optimality  of  the  transitions  from  one  mode  to  the  next 
is  ignored.  Comparisons  of  the  parametric-optimal  trajec- 
tories (formed  by  piecing  together  fuel-minimal  segments), 
and  fully  optimal  profiles  (determined  by  means  of  variational 
optimization),  are  presented  in  Chapter  5.  There  it  is 
demonstrated  that  parametrically  optimized  flight  segments 
provide  for  nearly  optimal  fuel  minimization. 

The  cruise,  climb,  and  descent  flight  modes  are 
discussed  in  Sections  4.2,  4.3,  and  4.4,  respectively.  The 
cruise  segment  is  examined  first,  as  the  terminal  conditions 
for  the  climb  mode  and  the  initiation  point  for  the  descent 
phase  are  determined  by  the  cruising  flight  conditions.  A 
summary  of  the  chapter  is  then  provided  in  Section  4.5. 


4.2  CRUISING  FLIGHT  PARAMETERS 

A straightforward  method  of  examining  the  fuel 
efficiency  of  a trajectory  is  to  divide  the  trajectory  path 
length  by  the  amount  of  fuel  consumed.  The  resultant  value, 
commonly  called  specific  range,  SR,  indicates  the  nautical 
miles  flown  per  pound  of  fuel.  The  instantaneous  value  of 
specific  range  at  any  point  on  the  trajectory  is  given 
approximately  as 


49 


(4.2-1) 


where  V is  the  aircraft's  velocity  in  knots  and  w is  the 
fuel  flow  rate  in  pounds  per  hour. 

The  expression  in  Eq.  (4.2-1)  is  approximate,  in  that 
it  assumes  that  the  velocity  vector  is  parallel  to  the  ground 
i.e.,  that  the  flight  path  angle  is  zero.  This  is  the  same 
approximation  that  is  made  in  deriving  the  maximum  range  flight 
defined  by  the  Breguet  equations  (Ref.  52).  As  the  flight 
path  angle  during  optimal  cruise  is  on  the  order  of  a few  hun- 
dredths of  a degree,  the  approximation  is  a good  one. 

Fuel-optimal  cruising  flight  can  be  achieved  by 
flying  at  conditions  such  that  the  aircraft's  instantaneous 
value  of  specific  range,  as  given  by  Eq.  (4.2-1),  is  maxi- 
mized. The  maximum  value  of  specific  range  is  determined  by 
the  point  at  which  V/w  attains  a maximum,  as  shown  in  Fig.  8. 
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It  is  shown  in  Refs.  1 and  2 that  maximum  specific 
range  for  a turbojet  aircraft  can  be  achieved  by  flying  at 
constant  values  of  W/6  (weight  over  atmospheric  pressure 
ratio)  and  Mach  number.  A plot  of  pressure  ratio  (6)  vs. 
altitude  is  presented  in  Fig.  9a  (Ref.  33).  At  sea  level, 

6 equals  one.  As  altitude  increases  to  50,000  ft,  atmos- 
pheric pressure  decreases,  and  6 is  approximately  0.09.  As 
the  aircraft  burns  fuel,  weight  decreases,  and  altitude  must 
be  increased  to  keep  W/6  constant.  Profiles  of  altitude 
versus  weight  for  three  constant  levels  of  W/6  are  shown  in 
Fig.  9b. 


•I  PRESSURE  RATIO  VERSUS  AITITUOS  M W/l  VERSUS  AlTlTUOE  AN0  WEIGHT 


Figure  9 Optimal  Cruise  Parameters 

These  values  of  W/6  are  typical  of  the  optimal 
values  for  the  C-141A  aircraft.  At  a vehicle  weight  of 
320,000  lb  (the  maximum  C-141A  weight),  the  W/6  profiles 
indicate  that  the  fuel-minimal  ^ruise  altitude  is  on  the 
order  of  33,000  ft.  At  a weight  of  140,000  lb,  the  optimal 
cruise  altitude  is  near  50,000  ft.  It  is  for  this  reason 
that  typical  aircraft  trajectories  include  climb  and  descent 
phases.  The  climb  segment  brings  the  vehicle  near  the 
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minimum- fuel  cruise  altitude,  and  the  descent  segment  returns 
the  aircraft  to  low  altitude  prior  to  landing. 

Sections  4.2.1  through  4.2.2  present  the  optimal 
cruise  parameters,  W/6  and  Mach  number,  for  the  C-141A. 

Nominal  parameter  values  and  the  effects  of  non-standard 
atmospheric  conditions  and  vehicle  configurations  on  these 
parameters  are  discussed. 

4.2.1  Nominal  Configuration 

Specific  range  contours  for  the  C-141A  aircraft 
model  employed  in  this  study  are  shown  in  Fig.  10a.  Those 
contours  are  based  upon  a nominal  aircraft  center  of  gravity 
(<  .g.  = 0.241  MAC),  a vehicle  weight  of  180,000  lb,  standard 
aimospheric  temperatures,  and  zero  wind  velocity.  The  specif- 
i<  range  values  presented  in  this  figure  neglect  trim  drag, 
as  outlined  in  Appendix  B. 

The  maximum  value  of  specific  range  occurs  for  the 
following  cruise  conditions: 

• Mach  number  ■ 0.727 

• W/6  - 1.24  x 106  lb 

At  a vehicle  weight  of  180,000  lb,  the  maximum  specific 
range  is  seen  to  be  0.0565  nm/lb.  For  other  vehicle  weights, 
the  maximum  specific  ranges  are  presented  in  Fig.  10b.  Also 
plotted  here  are  the  altitudes  which  correspond  to  maximum 
specific  range  for  W/6  equal  to  the  optimum  value, 

1.24*  106  lb. 

The  maximum  specific-range  numbers  quoted  in  this 
report  are  based  upon  the  horizontal  flight  path  assumption 
described  in  Section  4.2;  however,  in  order  to  follow  a 
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<1  SPECIFIC  NANCE  CONTOURS 


■•1*1*1 


Figure  10  C-141A  Nominal  Specific  Range  Parameters 

cruise-climb  ( constant-W/ 6 ) profile,  typical  values  of 
flight  path  angle  for  the  C-141A  are  on  the  order  of  0.02 
deg.  At  a vehicle  weight  of  229,000  lb,  the  horizontal- 
cruise  specific  range  is  0.04432  nm/lb,  while  the  climbing- 
cruise  is  0.04401  nm/lb.  Thus,  the  horizontal  flight 
instantaneously  has  a value  of  specific  range  which  is 
0.7  percent  greater  than  that  of  the  cruise-climb.  This 
percentage  is  equal  to  the  percentage  increase  in  energy 
which  must  be  expended  to  gain  altitude  on  the  climbing 
trajectory.  If  the  horizontal  flight  path  is  maintained, 
however,  the  aircraft  will  no  longer  be  flying  at  t he 
optimal  cruise  conditions,  and  specific  range  will  fall  far 
below  that  obtained  on  the  climbing  profile. 


53 


AITITUOC  (»*> 


4.2.2  Center-of-Gravity  Shift  Effects 


As  the  vehicle  center-of -grav ity  (c.g.)  moves  aft, 
the  pitching  moment  supplied  by  the  aircraft  tail  for  trim 
is  decreased.  Thus,  the  magnitude  of  negative  lift  supplied 
by  the  tail  is  decreased,  and  the  wing  lift  required  for 
level  flight  is  reduced  by  an  equal  amount.  As  the  lift  of 
the  wing  decreases,  so  does  the  wing's  contribution  to  total 
aircraft  drag.  In  order  to  examine  the  effects  of  specific 
range  caused  by  c.g.  location  shifts,  the  vehicle  model  which 
incorporates  tail-lift  effects  (see  Appendix  B)  is  employed 
for  the  present  and  all  ensuing  off-nominal  specific  range 
eva luat ions . 

The  optimal  value  of  W/6  to  maximize  specific  range 
is  shown  in  Fig.  11  as  a function  of  aircraft  c.g.  location. 
The  values  of  (w/5)0pt  have  been  computed  at  forward  (18 
percent  MAC),  mid  (24  percent  MAC),  and  aft  (34  percent 
MAC)  center-of-gravit y locations.  As  the  specific  range 
computations  are  performed  at  discrete  altitudes  (2,500  ft 
increments)  and  vehicle  weights  (10,000-lb  increments),  there 
is  some  uncertainty  in  the  resultant  mean  value  of  (W/6) 

opt 

For  this  reason,  the  standard  deviation  associated  with  the 
data  points  is  indicated.  The  variation  in  (w/5)opt  small 

(W/6opt  = 1 ’ 18  x 106  lb  at  forward  c.g.,  122  x 106  lb  at  aft 
c.g.)  and. is  within  the  uncertainty  of  the  estimates.  Thus, 

for  most  practical  purposes,  the  variation  in  ( W / 6 ) t with 
c.g.  location  shifts  can  be  neglected.  Also,  the  optimal 
cruising  Mach  number  remains  constant  at  a value  of  0.73  (a 
slight  increase  from  the  0.725  computed  without  trim  drag 
e f foots  ) . 


A plot  of  maximum  specific  range  vs.  aircraft 
weight  is  presented  in  Fig.  12.  These  specific  ranges  are 
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derived  from  the  aerodynamic  model  which  incorporates  the 
trim  lift/drag  effects.  Changes  in  these  values  of  spe- 
cific range  with  variation  in  c.g.  location  are  shown  in 
Fig.  13.  These  percentage  deviations  from  the  nominal  values 
are  independent  of  aircraft  weight.  Also,  the  percentage 
variations  in  specific  range  are  seen  to  be  1 i near  with  c.g. 
location  shifts.  The  maximum  deviation  in  specific  range, 
as  the  center  of  gravity  moves  from  the  forward  to  aft  loca- 
tion, is  about  5.7  percent. 


Figure  13  Effects  of  c.g.  Location  on  Specific  Range 
*1.2.3  Aircraft  Drag  Coefficient  Effects 

The  major  emphasis  of  this  study  is  concentrated  on 
trajectory  profile  shaping  to  minimize  fuel  consumption.  Air- 
craft con  1 igurat ion  alterations  also  could  reduce  fuel  use. 

In  particular,  reduced  aircraft  drag  serves  to  increase 
the  vehicle's  specific  range.  Each  aircraft  in  the  C-141A 
fleet  has  values  of  Cp  which,  in  all  probability,  are  not 
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exactly  equal  to  the  model  values  of  this  study.  The  effects 
of  those  Individual  drag  variations  on  fuel  performance  can 
bo  determined  from  the  following  results. 

As  part  of  the  proposed  C-141A  fleet  "stretching" 
operations,  aircraft  modifications  also  would  be  undertaken 
to  reduce  total  drag.  Current  projections  suggest  a drag 
reduction  of  up  to  6 percent  for  the  stretched  version,  the 
C-141B  (Ref.  53).  The  effects  of  ± 10  percent  variations  in 
CD  on  the  cruise-phase  specific  range  determined  in  this  sec- 
tion can  be  used,  therefore,  to  estimate  the  improvement  in 
performance  associated  with  the  stretching  modifications. 

As  vehicle  drag  is  reduced,  the  engine  thrust  re- 
quired for  cruising  flight  is  lowered  by  an  equivalent 
amount.  Because  fuel  consumption  is  monotonically  increas- 
ing with  thrust  (see  Fig.  B-l),  this  in  turn  lowers  the  air- 
craft’s fuel  consumption  rate.  The  cruising  flight  param- 
eters, W/6  and  Mach  number,  are  presented  in  Table  1 as 
functions  of  Cp  variations.*  The  optimum  Maeh  number  change 
for  + 10  percent  drag  variations  is  small  (presumably  because 
the  transonic  drag-rise  Mach  number  is  not  changed  signifi- 
cantly by  uniform  changes  in  Cp).  However,  the  values  of 
W/6  for  optimum  specific  range  indicate  a definite  decrease 
with  increasing  Cp.  Thus,  as  drag  increases,  the  cruising 
altitudes  associated  with  maximum  specific  range  are  lowered. 
Conversely,  as  Cp  is  improved,  maximum  specific  range  is 
achieved  by  flying  at  higher  altitudes. 


♦Percentage  variations  in  Cd  correspond  directly  to  per- 
centage variations  in  thrust  (at  a fixed  fuel  flow  rate) 
of  opposite  sign. 
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TABLE  1 


OPTIMUM  CRUISE  PARAMETERS  VS. 
AIRCRAFT  DRAG  COEFFICIENT  (CD) 


CD  Variations 

■mm 

(Mach  Number) 

opt 

-10  percent  Cq 

1 . 28x10®  lb 

0.74 

Nominal 

1.19xl06  lb 

0.73 

+10  percent  Cq 

l.lOxlO6  lb 

0.72 

The  changes  in  fuel  performance  caused  by  drag 
variations  are  presented  in  terms  of  the  aircraft's 
maximum  specific  range,  Fig.  14.  A reduction  in  drag  of 
10  percent  produces  an  increase  in  specific  range  of  approxi- 
mately 10  percent  for  all  vehicle  weights.  Increases  in  CD 
demonstrate  corresponding  decreases  in  fuel  economy.  The 
predicted  6 percent  reduction  in  the  value  of  CD  for  the 
stretched  C-141B  would,  therefore,  significantly  improve  the 
vehicle'^  fuel  efficiency. 


Figure  14  Specific  Range  as  a Function 
of  Cd  Variation 
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4.2.4  Atmospheric  Effects 


In  this  section,  the  effects  of  variations  in  atmo- 
spheric conditions  on  fuel-optimal  cruising  flight  are 
described.  The  atmospheric  variations  which  are  of  interest 
here  are  the  deviations  in  temperature  (from  the  1962 
Standard  Atmosphere  Model)  and  non-zero  wind  velocities.  It 
is  reasonable  for  these  parameters  to  change  over  the  course 
of  a single  mission;  by  accounting  for  known  variations, 
fuel-optimality  can  be  assured  along  the  cruise  profile. 

The  effects  of  temperature  variations  on  cruise  are  found  to 
be  small  compared  to  the  vehicle  configuration  effects  dis- 
cussed in  the  previous  sections.  Table  2 presents  the  opti- 
mal cruise  parameters  associated  with  atmospheric  temperature 
variations  of  ±20°F.  The  effect  of  these  variations  on  fuel 
performance  (specific  range)  also  is  shown  here.  The  optimal 
cruise  altitudes,  as  defined  by  the  parameter,  W/S,  change 
insignificantly  with  temperature  variations,  while  the  "best” 
Mach  number  decreases  with  increasing  temperature.  The  re- 
sultant variation  in  maximum  specific  range  caused  by  the 
temperature  variations  is  on  the  order  of  only  1 percent. 

TABLE  2 

EFFECTS  OF  TEMPERATURE  DEVIATIONS  ON  CRUISE  PARAMETERS 


Temperature 

Variations 

(°F) 

(wy6)opt  (lb) 

(Mach  Number )q 

Specific  Range 
Variat ion 
( percent ) 

! -20 

1.18xl06 

0.74 

+0.81 

Nominal 

1.19xl06 

0.73 

— 

+20 

1.19xl06 

0.72 

CM 

O 

1 

Atmospheric  winds,  on  the  other  hand,  have  a large 
effect  on  optimal  cruise  conditions.  In  this  study,  wind 


59 


Figure  15  Optimal  Cruise  Mach  Number 
vs.  Wind  Velocity 

values  of  -300  fps  (head  wind)  to  +300  fps  (tail  wind)  are 
examined.  Optimal  cruise-phase  Mach  numbers  are  presented 
as  a function  of  wind  velocity  in  Fig.  15.  For  a head  wind, 
the  best  cruise  Mach  number  increases  from  its  nominal  value 
of  0.73.  With  a tail  wind,  the  desired  Mach  number  is  re- 
duced, but  the  optimal  value  of  (W/6)  for  all  wind  veloci- 
ties  remains  constant  at  a value  of  about  1.19  x 10  lb. 

The  observed  optimal  Mach  number  trends  with  wind 
variation  can  be  explained  as  follows.  In  a tail  wind  envir- 
onment, vehicle  ground-relative  speed  increases  if  Mach  num- 
ber is  not  changed.  Therefore,  a reduction  in  Mach  number 
to  reduce  fuel  flow  rate  can  be  accomplished  while  still  main- 
taining good  ground  speed.  In  a head  wind  case,  ground  speed 
is  reduced.  Mach  number  is  increased  in  this  situation,  in 
order  that  the  total  fuel  consumption  for  the  flight  is 
minimized. 
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Figure  16  presents  the  maximum  values  of  specific 
range  for  the  C-141A  aircraft  for  several  wind  environments. 
Specific  range  curves  for  wind  speeds  between  -300  and  +300 
fps  are  plotted  for  wind  speed  increments  of  100  fps.  In  the 
maximum  wind  speed  variation  cases,  the  specific  range  varia- 
tions are  as  large  as  40  percent.  Because  the  variations 
are  this  large,  it  would  be  reasonable  (in  some  situations)  to 
increase  trajectory  path  length  in  order  to  maximize  tail 
winds  and  minimize  head  winds  along  the  flight  path.  Trade- 
offs would  have  to  be  performed  to  insure  that  the  fuel 
reductions  afforded  by  increasing  specific  range  are  not  off- 
set by  the  increase  in  fuel  needed  to  fly  the  longer  path 
length.  In  all  cases,  however,  the  cruise  profile  should  be 
flown  at  optimal  values  of  the  cruise  parameters  to  provide 
the  maximum  specific  range  for  the  wind  conditions  encountered. 
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Figure  16  Specific  Range  as  a Function 
of  Wind  Velocity 
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4.3 


CLIMB  PROFILES 


During  the  climbing  segment  of  flight,  the  aircraft's 
altitude  and  Mach  number  are  changed  from  their  takeoff  values 
to  the  values  which  define  the  initiation  of  the  cruise  phase. 
Thus,  the  vehicle's  total  energy  (kinetic  plus  potential)  is 
greatly  increased  in  this  flight  mode.  In  this  section, 
flight  profiles  which  minimize  the  fuel  used  to  accomplish 
these  energy  changes  are  determined. 

The  conditions  defining  the  end  of  the  climb  seg- 
ment are  obtained  from  the  optimal  cruise  parameters  described 
in  Section  4.2.1.  When  these  cruise  conditions  (W/6  and 
Mach  number)  are  met,  the  aircraft  profile  switches  from 
the  climb  to  the  cruise  phase.  In  this  study,  the  aircraft 
altitude  at  the  beginning  of  climb  is  chosen  to  be  2,000  ft, 
and  the  initial  velocity  is  340  fps  true  airspeed.  At  the 
2,000-ft  initial  altitude,  this  corresponds  to  a value  of 
195  KIAS.  These  conditions  are  similar  to  those  of  Ref.  32, 
for  a vehicle  which  has  completed  takeoff  and  has  fully  re- 
tracted the  flaps.  Thus,  the  climb  profiles  simulated  in 
this  study  assume  a clean  vehicle  configuration  and  well- 
defined  initial  and  terminal  cruise  conditions. 

A conventional  C-141A  climb  profile  is  described  in 
Section  4.3.1.  Optimal  (minimum-fuel ) climb  segments  for  a 
range  of  initial  vehicle  weights  are  presented  in  Section 
4.3.2.  The  utility  of  flying  minimum-fuel  climb  modes  also 
is  demonstrated  in  Chapter  5,  where  complete  trajectory 
climb,  cruise,  and  descent  fuel  consumption  is  discussed. 
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4.3.1  Conventional  Profiles 


The  conventional  C-141A  climb  profiles  simulated  in 
this  study  closely  approximate  the  handbook  schedules 
documented  in  Ref.  32.  The  climb  profiles  consist  of  two 
constant- IAS  modes  followed  by  a climb  at  constant  Mach 
number.  The  values  of  climb  IAS  and  Mach  number,  as  well 
as  the  mode  switching  points,  are  summarized  in  Table  3. 

All  the  climb  modes  are  flown  with  throttle  set  for  normal 
rated  thrust  (NRT).  This  corresponds  to  full  throttle 
setting  in  this  study. 


TABLE  3 

CLIMB  SCHEDULE  FOR  CONVENTIONAL  PROFILES 


Mode  Condition 

Mode  Termination 
Point 

Throttle 

Setting 

n 

IAS  - 250  kt 

h - 10,000  ft 

1.0  (NRT) 

IAS  * 280  kt 

Mach  Number  = 0.720 

1.0 

Mach  Number  ■ 0.720 

W/6  - 1.24  x 106  lb 

1.0 

A plot  of  altitude  vs.  range  for  a conventional 
climb  is  presented  in  Fig.  17.  The  profile  is  generated  for 
a C-141A  vehicle  with  an  initial  weight  (at  2,000-ft  altitude) 
of  260,000  lb.  The  climb  trajectory  has  a path  length  of 
228  nm  and  consumes  10,691  lb  of  fuel.  During  the  constant- 
I AS  climb  modes,  the  aircraft's  altitude  rate  decreases  from 
3,000  to  750  ft/min.  Altitude  rate  during  the  constant- 
Mach  number  phase  is  on  the  order  of  1,000  ft/min. 

The  climb  profile  is  shown  to  change  slope  at  the 
switch  point  between  the  constant-IAS  and  constant-Mach 
number  modes.  This  change  in  slope  at  the  IAS/Mach  number 
twitch  point  also  is  found  to  occur  in  the  fuel-optimal  climb 
profiles  described  in  the  following  section. 
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Figure  17  Conventional  Climb  Profile 
4.3.2  Optimal  Profiles 

Conventional  climb  profiles  are  defined  in  terms  of 
constant-IAS  and  constant-Mach  number  modes,  in  order  to 
ease  the  task  of  piloting  the  aircraft.  In  the  present 
ptudy  of  fuel-minimal  climbs,  these  same  flight  modes  are 
preserved,  so  that  flying  a fuel-efficient  climb  will  be 
no  more  difficult  than  following  a conventionally  defined 
profile. 


The  parameteric  optimization  techniques  described 
in  Section  3.5  have  been  employed  for  improving  fuel  consump- 
tion during  climb.  Rather  than  using  two  constant-IAS  climb 
modes,  only  one  constant-IAS  mode  followed  by  a constant-Mach 
number  mode  is  presumed  for  the  generation  of  optimal  climb 
profiles.  This  enables  the  use  of  two-parameter  parametric 
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optimization  methods  to  determine  fuel-minimal  climbs.  The 
parameters  which  are  to  be  optimized  are  the  constant  values 
ol  IAS  and  Mach  number  to  be  flown  during  climb. 

Parametric-optimal  climb  profiles  have  been  gener- 
ated for  C-141A  aircraft  with  initial  weights  of  200,000, 
200,000,  and  320,000  lb.  The  "best"  values  of  airspeed  and 
Mach  number,  as  a function  of  these  initial  aircraft  weights, 
are  presented  in  Fig.  18a.  A Mach  number  of  0.721  provides 
fuel-minimum  climbs  for  all  the  trajectories  simulated.  Best 
indicated  airspeed  increases,  however,  with  increasing  ini- 
tial weight.  For  an  initial  weight  of  200,000  lb,  the  "best" 
constant  value  of  airspeed  is  191  kt . This  value  increases 
to  245  kt  for  a 320,000-lb  vehicle.  Also  shown  on  the  fig- 
ure are  the  altitudes  at  which  the  switch  from  constant-IAS 
to  constant-Mach  number  is  performed. 

Although  the  heavy  aircraft  terminates  the  climb 
phase  at  a lower  altitude  than  does  a lighter  weight  vehicle*, 
it  consumes  more  fuel  in  the  climb  segment.  This  is  demon- 
strated by  the  plot  of  climb-phase  fuel  consumption  vs. 
initial  vehicle  weight  of  Fig.  18b. 

Altitude-vs. -range  profiles  for  the  three  optimal 
climb  trajectories  are  presented  in  Fig.  19.  It  is  seen 
here  that  as  initial  aircraft  weight  increases,  both  the 
total  climb  range  and  the  portion  of  flight  spent  in  the 
constant-Mach  number  mode  increase.  However,  the  general 
profile  shapes  are  similar  to  the  conventional  profile  pic- 
tured in  Fig.  17. 
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Optimal  Climb  Conditions 
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Figure  19  Parametric-Optimum  Climb  Profiles 
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The  optimal  climb  profile  for  the  260,000-lb  air- 
craft traverses  a range  of  211  nm  during  climb  and  consumes 
10,277  lb  of  fuel.  This  is  414  lb  less  than  the  fuel  needed 
for  the  conventional  climb  schedule.  However,  the  conventional 
profile  travels  17  nm  further  in  the  climb  segment.  Thus, 
the  specific  range  (nm  traveled/lb  fuel  consumed)  is  nearly 
identical  for  the  fuel-optimal  and  conventional  trajectories. 
The  major  advantage  of  fuel-optimal  climbs  is  demonstrated 
by  the  short-range  trajectories  where  the  fuel  minimal  pro- 
files do  not  contain  a cruise  phase.  The  fuel  savings 
ai forded  by  optimal-climb  profiles  for  these  trajectories 
are  discussed  in  detail  in  Section  5.2. 


4.4  DESCENT  PROFILES 

The  results  of  the  fuel-minimum  cruise  phase 
analyses,  Section  4.2,  demonstrate  that  the  best  C-141A 
cruising  altitudes  are  typically  greater  than  30,000  ft. 

The  descent  flight  segment  serves  to  reduce  altitude  from 
the  cruise  levels  to  a value  at  which  the  landing  phase 
cun  be  initiated.  In  this  study,  the  altitude  at  which  the 
descent  phase  terminates  is  chosen  to  be  2,000  ft  above  sea 
level.  This  is  the  same  value  of  altitude  that  has  been 
used  to  initiate  the  climb  phase  analyses.  T e three  fuel- 
optimal  trajectory  segments--climb,  cruise,  and  descent-- 
cin  be  pieced  together,  therefore,  to  provide  complete  tra- 
jectories which  begin  and  end  at  a 2,000-ft  altitude. 

As  in  the  case  of  the  C-141A  climb  profiles,  both 
conventional  and  fuel-optimal  descent  profiles  are  examined 
in  this  section.  The  fuel-optimal  descent  profiles  are 
shown  to  have  longer  flight  path  ranges  than  their  con- 
ventional counterparts.  This  extra  descent  range  reduces 
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the  duration  of  the  cruise  segment  for  a given  mission,  and 
results  in  a lower  value  of  net  fuel  consumption.  A descrip- 
tion of  the  conventional  descent  profiles  is  given  in  Section 

1.4.1,  and  the  fuel-optimal  descents  are  reviewed  in  Section 

1.4.2. 

4.4.1  Conventional  Profiles 

The  conventional  C-141A  descent  profiles  simulated 
in  this  study  are  representative  of  the  "Enroute  Descent" 
schedules  of  Ref.  32;  however,  the  descents  simulated  here 
employ  zero  thrust  (idle)  throttle  settings.  Actual  C-141A 
descents  are  typically  flown  with  posi five  throttle  settings, 
and  pilots  frequently  increase  Mach  number  on  descent  in 
order  to  begin  the  scheduled  constant-IAS  mode  more  quickly. 
Thus,  the  simulated  conventional  descents,  by  which  the 
I'ue  L-opt  imal  profiles  of  this  study  are  judged,  are  themselves 
optimistic. 

The  descent  profiles  each  consist  of  a constant -Mach 
number  mode  followed  by  two  constant-IAS  modes.  Thus,  the 
descent  schedule  contains  the  same  modes  as  the  climb  schedule, 
but  they  are  flown  in  reverse  order  for  descent.  The  values 
of  IAS  and  Mach  number  flown  on  the  conventional  descent 
schedule  are  presented  in  Table  4.4-1.  Also  included  in  the 
(able  are  the  mode  switching  points  and  the  descent  phase 
i hrot  1 1 e set  t i ngs . 


TABLE  4 

DESCENT  SCHEDULE  FOR  CONVENTIONAL  PROFILES 


Mode  Condition 

Mode  Termination 
Point 

Thrott le 
Set  t i nn 

1 . 

Mach  Number  = 

0.725 

Mach  Number  < 0.725 

0 (Idle) 

2. 

IAS  = 300  kt 

h = 10,000  ft 

0 

3 . 

IAS  = 250  kt 

h = 2,000  ft 

0 
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During  descent,  the  throttle  is  set  to  the  idle 
position,  providing  zero  thrust.  This  does  not  imply,  how- 
ever, that  fuel  flow  is  also  zero,  us  a certain  amount  of  fuel 
is  required  just  to  keep  the  powerplants  operating.  The  fact 
that  fuel  flow  rate  is  not  zero  at  idle  throttle  setting  is 
demonstrated  by  the  fuel  flow-vs . -thrust  curves  of  Fig.  B-l  . 

At  the  zero-thrust  level,  the  fuel  flow  curves  demonstrate 
positive  values  of  flow  rate  which  increase  with  increasing 
Mach  number. 

A plot  of  altitude  vs.  range  for  a conventional 
descent  is  presented  in  Fig.  20.  The  profile  represents  a 
simulated  descent  for  a C-141A  vehicle  cruising  at  Mach 
0.740  at  an  altitude  of  35,000  ft.  At  the  initiation  of  the 
descent,  the  aircraft  weight  is  232,537  lb,  and  the  descent 
path  length  is  approximately  95  nm.  Descent  rate  during 
the  constant-Mach  number  phase  reaches  a maximum  value  of 
about  4,000  ft/min  and  reduces  to  approximately  1,500  ft/min 
during  the  constant-IAS  modes. 
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Figure  20  Conventional  Descent  Profile 
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Time  of  flight  and  fuel  consumption  values  for 
this  descent  are  shown  in  Fig.  21.  These  variables  are  plot- 
ted as  functions  of  initial  altitude  on  the  descent  profile. 

For  example,  at  an  altitude  of  20,000  ft,  approximately  9.f»min 
of  flight  time  and  800  lb  of  fuel  »"ll  be  expended  along  the 
remainder  of  the  descent  profile.  .he  larger  increase  in 
flight  time  for  altitudes  above  33,000  ft  occurs  on  this  pro- 
lile  beeauso  of  the  transition  between  the  cruise  and  descent 
segments  of  flight.  This  transition  period  is  also  discern  - 
able  on  the  altitude  profile  of  Fig.  20.  The  transition 
between  cruise  and  descent  spans  the  first  10  nm  of  the 
descent  profile  and  contributes  an  altitude  decrease  of  about 
2,000  ft. 
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Figure  21  Conventional-Descent  Fuel 

Consumption  and  Flight  Time 
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Although  most  of  the  performance  results  concerning 
i omplete  trajectory  profiles  are  discussed  in  Chapter  5, 

1 he  effects  of  using  non-zero  throttle  setting  during  descent 
on  total  fuel  consumption  are  presented  here.  Total  fuel 
consumption  for  four  different  1,000-nm  missions  are  shown 
in  Table  5.  (Initial  weight  of  the  vehicle  is  260,000  lb.) 

The  missions  are  classified  according  to  the  constant  altitude 
and  constant  Mach  number  flown  during  the  cruise  flight  seg- 
ment. The  two  columns  of  fuel  consumption  data  correspond 
to  flying  each  mission  with  either  idle  throttle  or  20-per- 
cent throttle  setting  during  descent.  For  both  types  of 
descent,  the  Mach  number  and  indicated  airspeed  profiles 
flown  are  those  of  Table  3. 


TABLE  5 

MISSION  FUEL  CONSUMPTION  FOR  1,000-NM  FLIGHT  PATHS 
WITH  IDLE  AND  20-PERCENT  THROTTLE  SFTTING  DESCENTS  * 


Cruise  Conditions 

— 
Fuel  Consumption  (lb) 

Altitude  (ft) 

Mach  Number 

Idle  Throttle 

20-percont 
Thro  1. 1 1 e 

31,000 

0.740 

30,032 

30,956 

35,000 

0.740 

28,452 

29,520 

31,000 

0.767 

31 , 199 

32 . 022 

35,000 

0.767 

29,233 

30,237 

Use  of  20-percent  throttle  setting  on  the  descent 
profiles  causes  an  increase  in  fuel  consumption  of  about  H00 
to  1,100  lb  for  these  1,000-nm  missions.  Thus,  scheduling 
of  idle  throttle  on  descent  does  tend  to  reduce  overall  fuel 


♦The  climb  and  cruise  profiles  are  identical,  however,  the 
idle  descents  begin  later  in  the  flight  than  do  the  20  percent 
throttle  descents. 
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consumption  on  conventional  trajectories.  In  the  following 
section,  it  is  demonstrated  that  idle  throttle  setting  also 
is  a characteristic  of  the  fuel-optimal  descent  profiles. 
These  results  also  suggest  that  the  "conventional”  profiles 
using  idle-throttle  descents,  which  are  compared  with  opti- 
mal profiles  in  the  next  chapter,  actually  are  optimistic. 

4.4.2  Optimal  Profiles 

If  fuel  flow  rate  at  idle  throttle  setting  were 
zero,  an  optimal  descent  profile  would  be  one  which  maximizes 
gliding  path  length.  By  providing  maximum  range  during 
descent,  the  duration  of  the  cruise  phase  of  flight  would 
be  reduced,  resulting  in  a minimum  value  of  fuel  consump- 
tion for  the  complete  mission.  Although  fuel  use  is  not  zero 
at  idle  throttle,  it  is  substantially  lower  than  the  values 
experienced  for  typical  cruising  flight.  For  this  reason, 
the  investigation  of  near-optimal  descent  profiles  described 
in  this  section  concentrates  on  determining  maximum  path- 
length,  idle-throttle  descents. 

It  is  well  known  that  maximum  path  length  of  a 
gliding  aircraft  is  achieved  by  flying  at  the  vehicle’s 
maximum  value  of  lift-to-drag  ratio  (L/D).  At  idle  throttle, 
the  net  thrust  of  the  C-141A  is  zero,  so  the  aircraft  is,  in 
effect,  flying  as  a glider.  The  problem  of  determining 
near-optimal  descent  profiles  is  to  find  the  maximum  L/D 
gl iding  fl ight  path. 

The  conventional  descent  profiles  described  in 
Section  4.4.1  are  not  difficult  for  a pilot  to  fly  because 
they  are  defined  in  terms  of  constant-Mach  number  and  con- 
stant-IAS  modes.  Near-maximum  L/D  descent  profiles  for 
the  C-141A  can  be  achieved  by  maintaining  a constant  aircraft 
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pitch  angle  (0).  For  the  C-141A  aircraft  modeled  in  this 
study,  a pitch  angle  of  0.935  degrees  provides  a near-maximum 
average  value  of  L/D,  equal  to  about  19.5  along  the  descent 
profiles.  The  validity  of  the  assumptions  (small  fuel  flow 
rate  at  idle  throttle  setting  and  near-maximum  L/D  at  constant 
pitch  angle)  is  demonstrated  in  Chapter  5 of  this  report.  It 
is  shown  that  the  descent  portion  of  the  fully  optimal  tra- 
joctory  profiles,  as  determined  by  complex  variational  opti- 
mization methods,  is  nearly  identical  to  the  simple  near- 
optimal  profiles  described  here. 

The  increased  glide  range  afforded  by  the  constant-0  , 
maximum-L/D  descent  schedule  over  the  conventional  descent 
schedule  is  demonstrated  by  Fig.  22.  Both  descent  profiles 
are  initiated  at  an  altitude  of  about  40,000  ft  with  an 
initial  vehicle  weight  of  approximately  234,000  lb.  The 
optimal  descent  traverses  a path  length  of  143  nm,  while  the 


Figure  22  Comparison  of  Optimal  and  Conventional 
Descent  Profiles 
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conventional  descent  flies  only  108  nm.  Thus,  the  trajectory 
which  incorporates  a conventional  descent  must  fly  35  nm 
longer  in  the  cruise  phase. 

The  fuel  consumption  for  the  150-nm  profiles  plotted 
is  1,552  lb  for  the  optimal  descent  and  2,161  lb  for  the  con- 
vontional  descent.  For  a 1,000-nm  mission  with  a typical 
total  fuel  requirement  of  30,000  lb,  this  descent  profile's 
savings  of  609  lb  represents  a net  reduction  in  total  fuel 
use  of  about  2 percent.  Thus,  flying  optimal  descent  pro- 
files is  an  excellent  means  of  reducing  fuel  consumption  for 
the  C-141A  fleet. 

4 . 5 CHAPTER  SUMMARY 

In  this  chapter,  fuel  consumption  along  the  three 
major  flight  path  segments  — climb,  cruise , and  descent  — 
has  been  examined.  By  applying  the  optimization  techniques 
described  in  Chapter  3,  minimum-fuel  profile  shapes  are 
determined  for  these  flight  segments.  These  fuel-minimal 
profiles  are  shown  to  provide  significant  fuel-use  reduc- 
tions compared  to  currently  flown,  conventional  flight  paths. 

The  fuel-optimal  cruise  profile  is  found  to  be 
characterized  by  constant  values  of  Mach  number  and  W/6 
(weight  over  atmospheric  pressure  ratio).  As  fuel  is  burned, 
aircraft  weight  decreases,  and  altitude  must  be  increased 
to  main1 Jin  constant  W/6.  For  this  reason,  the  optimal 
cruise  trajectories  slowly  climb  along  the  mission  profiles. 

The  effects  of  off-nominal  vehicle  and  atmospheric 
conditions  on  the  optimal  cruise  profiles  also  are  presented 
in  this  chapter.  The  changes  in  fuel  efficiency  are 
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indicated  by  the  specific  range  (nm/lb  fuel)  associated 
with  the  cruise  profiles.  As  the  center  of  gravity 
of  the  aircraft  is  moved  aft,  the  cruise-phase  specific 
range  increases;  however,  the  optimal  values  of  Mach  number 
and  W/5  are  insensitive  to  these  c.g.  location  shifts. 

Head  winds  are  found  to  decrease  specific  range 
while  increasing  optimal  cruise  Mach  number,  and  the  opposite 
effects  are  true  of  tail  winds;  however,  wind  variations  have 
no  effect  on  the  best  value  of  W/6.  Atmospheric  temperature 
variations  have  only  a small  effect  on  cruising  specific 
range.  Temperature  variations  of  ±20°F  cause  a 2 percent 
change  in  the  aircraft's  specific  range. 

Decreases  in  aircraft  drag  coefficient.  Cp,  are 
shown  to  increase  both  the  cruise  specific  range  and  the 
optimal  value  of  W/6.  The  opposite  effects  result  for 
increasing  values  of  drag  coefficient. 

The  examination  of  fuel-minimal  climb  profiles  con- 
sists of  determining  the  best  values  of  indicated  air  speed 
(IAS)  and  Mach  number  to  fly  during  climb.  The  results  of 
this  analysis  indicate  that  the  optimal  climb  Mach  number. 
0.721,  is  the  same  for  all  vehicle  weights.  The  best  IAS 
varies,  however,  as  a function  of  vehicle  weight  from  191 
kt  for  a 200,000-lb  aircraft  to  245  kt  for  a 320.000-lb 
aircraft.  The  fuel  consumed  during  the  climb  also  increases 
substantially  with  increasing  vehicle  weight. 

The  optimal  descent  profiles  generated  in  this 
chapter  are  based  upon  flying  near-maximum  L/D  glides.  It 
is  shown  that  these  fuel-efficient  gliding  descents  can  be 
achieved  with  idle  throttle  setting  and  constant  aircraft 
pitch  angle.  For  a typical  descent  from  40,000  ft  to  2,000  ft  . 
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the  optimal  descent  profile  provides  a fuel  consumption 
value  609  lb  less  than  that  obtained  by  flying  a conventional 
enroute  descent . 

The  optimal  flight  segments  defined  in  this  chapter 
can  be  utilized  individually  to  provide  some  reduction  in 
overall  fuel  consumption  or  they  can  be  pieced  together  to 
provide  near-optimal  integrated  trajectory  profiles.  The 
efficiency  of  these  near-optimal  profiles  as  compared  to 
fully  optimal  complete  flight  paths  is  determined  in  the 
following  chapter. 


V. 


COMPARISON  OF  OPTIMAL  AND  CONVENTIONAL 


INTEGRATED  FLIGHT  PATHS 


5.1  OVERVIEW 

The  results  presented  in  Chapter  4 demonstrate 
the  potential  for  reducing  aircraft  fuel  consumption  on 
individual  flight  segments.  By  simply  piecing  optimized 
flight  segments  (climb,  cruise,  and  descent)  together, 
trajectory  profiles  which  are  more  fuel-efficient  than 
current  handbook  profiles  can  be  generated.  There  is 
no  guarantee,  however,  that  the  resulting  piecewise- 
optimal  trajectories  will  be  fully  optimal. 

Fuel  consumption  comparisons  for  conventional 
(handbook  defined),  piecewise-optimal , and  f ul ly-opt imal 
C-141A  flight  profiles  are  presented  in  this  section. 

For  convenience,  the  trajectories  are  categorized  according 
to  path  length:  short  range  (0-500  nm),  medium  range  (500- 

2,000  nm)  and  long  range  (2,000-5,600  nm).  These  trajec- 
tories are  presented  as  Sections  5.2,  5.3,  and  5.4,  respec- 
tively. The  comparisons  in  all  categories  demonstrate  that 
the  piecewise-optimal  trajectories  provide  fuel  consumption 
values  which  are  very  nearly  fuel-optimal.  Piecewise-optimal 
trajectories  tend  to  use  one  to  two  percent  more  fuel  than  do 
the  fully  optimal  flight  paths.  The  conventional  trajectory 
fuel  usage  is  approximately  10  to  20  percent  greater  than  the 
optimum  values.  These  percentages  vary  as  a function  of 
total  trajectory  length  and  handbook  cruise  conditions  (alti- 
tude and  Mach  number). 
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In  addition  to  the  fuel  consumed,  the  duration, 
in  hours,  is  recorded  for  all  trajectories.  The  current 
optimization  efforts  are  directed  toward  reducing  total 
fuel  usage,  not  total  cost  in  dollars.  The  actual  costs 
associated  with  each  profile,  are  tabulated  assuming  a fuel 
cost  of  45^/gallon  and  a time-related  cost  of  $400/hr. 

In  addition,  a range  of  fuel  and  time  costs  is  considered. 
With  this  choice  of  numbers,  some  of  the  conventional 
trajectories  which  appear  to  consume  excessive  amounts  of 
fuel  provide  for  total  dollar  costs  which  are  not  signifi- 
cantly different  from  the  fuel-minimum  profiles.  As  fuel 
costs  appear  to  be  increasing  more  rapidly  than  personnel 
and  other  time-related  expenses,  trajectory  fuel -opt imal i ty 
is  becoming  increasingly  significant.  The  computer  tools  and 
optimization  techniques  employed  for  fuel  minimization 
(see  Chapter  3 and  Appendix  A)  can  easily  handle  the 
"dol lar-optimizat ion"  problem.  The  means  by  which  this 
can  be  done  and  the  present  fuel-time  tradeoffs  are  dis- 
cussed in  Section  5,5.  A summary  of  the  chapter  is  pre- 
sented in  Section  5.6. 

5.2  SHORT-RANGE  FLIGHT  PATHS 

In  the  current  study,  flight  profiles  with  ranges 
of  500  nm  or  less  are  considered  to  be  "short-range."  These 
short-range  profiles  are  character lzeri  by  the  fact  that  the 
optimal  trajectories  contain  only  little  or  no  cruise  seg- 
ment. Almost  ail  fuel -use  reduction  for  these  trajectories 
is  achieved  during  the  climb  and  descent  flight  modes. 

The  parametric-optimum  a 1 1 i t ude- vs  -range  pro- 
file for  a 200-nm  trajectory  is  presented  in  Fig.  23a. 
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t-wto* 


a)  ALTITUDE  VS.  RANGE  PROFILE 


b>  FUEL  USE  VS.  CLIMB  IAS 


Figure  23  Minimum-Fuel  200-nm  Trajectories 

The  initial  aircraft  weight  is  260,000  lb,  and  the  trajec- 
tory begins  and  ends  at  an  altitude  of  2,000  ft.  Not  only 
does  the  profile  contain  no  cruise  segment,  but  the  constant- 
Mach  number  phase  of  the  parametric -opt imum  climb  segment 
(see  Section  4.3.2)  also  is  absent.  Because  the  trajectory 
shape*  is  different,  the  constant-IAS  climb  profiles  for  the 
short-range  case  may  be  different  from  the  values  obtained 
by  optimizing  the  climb  segment  alone. 

In  order  to  determine  the  parametric-optimum  200-nm 
range  profile,  fuel  use  vs.  indicated  air  speed  (IAS)  during 
climb  is  determined.  Fuel  use  as  a function  of  climb  IAS 
for  this  short-range  case  is  plotted  in  Fig.  23b.  The 
minimum  total  fuel  use,  7,515  lb,  is  obtained  with  a climb 
speed  of  279  kt . 

The  utility  of  parametric  optimization  techniques 
for  determining  minimum-fuel  trajectories  is  demonstrated 
by  the  comparison  of  fuel  consumption  for  conventional, 
paramet ric-opt ima 1 . and  fully  optimal  ( steepest-descent ) 
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technique)  profiles.  The  short-range  conventional  profile 
is  characterized  by  handbook  climb  and  descent  schedules  and 
a one-third  path-length  constant-altitude  cruise  (h = 23,000  It) 
Fuel  consumptions  for  the  three  trajectories  are  presented  in 
Table  6.  The  parametric  results  are  within  1.3  percent  of  the 
steepest-descent  defined  optimal,  while  the  conventional  pro- 
file uses  12.5  percent  more  fuel. 


TABLE  6 

FUEL  CONSUMPTION  FOR  200-NM 
TRAJECTORIES 


Trajectory 

Fuel  Use  (Lb) 

■aHglliBufliiMi 

Steepest  Descent 

7,420 

• 

Opt imum 

Parametric  Optimum 

7,515 

1 . 3 

Conventional 

8,347 

12.5 

(1/3  Level  Flight) 

Thrott le-vs . -range  and  thrott le-vs . -a  1 1 i tude  pro- 
files for  the  three  200-nm  trajectory  types  are  plotted  in 
Figs.  24a  and  24b,  respectively.  The  steepest -descen t 
throttle  profile  starts  to  throttle-down  from  maximum  at 
the  25-nm  range  point.  Thus,  this  reduction  occurs  during 
the  c 1 imb  segment  of  the  trajectory.  The  initial  conven- 
tional profile  throttle  reduction  occurs,  however,  at  th< 
Initiation  of  the  cru l so  flight  mode.  By  reducing  thrott !• 
setting  during  climb  the  steepest-descent  profile  burns 
less  fuel  than  does  the  parametric  trajectory. 

Both  the  parametric  and  steepest-descent  optima! 
trajectories  reduce  throttle  setting  to  idle  at  a range  o! 
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MUM 


a ) •)  throttle  vs.  range 


b)  THROTTLE  VS.  ALTITUDE 


Figure  24  Short-Range-Trajectory 
Throttle  Profiles 

approximately  100  nm.  The  conventional  trajectory  does  not 
employ  idle  throttle  until  about  135  nm  have  been  traversed. 

By  following  a near  maximum  L/D  descent  schedule,  the  optimal 
trajectories  have  a greater  glide  range,  which  contributes  to 
their  lower  fuel  usage.  Although  the  steepest-descent  profile 
provides  a 1.3  percent  fuel  savings  over  the  parametric  tra- 
jectory, the  on/off  nature  of  the  parametric  throttle  set- 
ting may  be  simpler  for  a pilot  to  fly.  On  the  other  hand, 
the  reduced-throttle  climb  of  the  steepest-descent  optimum 
may  be  advantageous  from  the  standpoint  of  minimizing  engine 
maintenance  costs.  In  either  case,  the  throttle  profile 
should  be  accompanied  by  the  proper  altitude-Mach  number 
profile,  shown  in  the  following  figure. 
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MACH  NO. 

•I  CONVENTIONAL  TRAJECTORY 


SWITCH  FtOM  MAX 
TO  MIN  THAOTTIE  • 


MACH  NO. 

bl  PARAMETRIC-OPTIMUM  TRAJECTORY 


SWITCH  TO 


MACH  NO. 


el  STEEPEST-DESCENT-OPTIMUM  TRAJECTORY 


Figure  25  200-nra  Altitude  vs. 

Mach  Number  Profiles 
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The  altitude-vs. -Mach  number  profile  of  the  con- 
ventional trajectory,  Fig.  25a,  is  not  as  smooth  as  those  of 
the  optimal  profiles.  The  "bumps”  in  the  profile  are  due  to 
the  numerous  changes  in  IAS  and  Mach  number  which  define  the 
handbook  climb  and  descent  schedules.  The  parametric  pro- 
file, Fig.  25b,  follows  the  constant  279-KIAS  climb  and  the 
maximum  L/D  descent  segments.  The  fully  optimal  trajectory, 
Fig.  25c,  contains  a constant-Mach  number  climb  and  a ne*r- 
maximum  L/D  descent  similar  to  that  of  the  parametric  pro- 
file. Although  the  two  optimal  climb  profiles  appear  to  be 
quite  different,  the  resultant  effects  on  fuel  reduction  are 
not  that  dissimilar. 

The  effectiveness  of  parametric  optimization  in 
determining  fuel-minimum  flight  profiles  also  is  evident 
from  the  500-nm,  260,000-lb  vehicle,  trajeclory  analyses. 
Fuel  consumption  values  for  this  path  length,  as  deter- 
mined by  parametric  and  variational  optimization  methods, 
are  shown  in  Table  7.  In  this  case,  the  parametric  optimum 
provides  for  fuel  consumption  within  1.1  percent  of  the 
true  minimum.  This  parametrically  defined  flight  path 
consists  of  segments  (constant- IAS  climb,  constant-Mach 
number  climb,  constant-altitude  cruise,  and  constant- 
pitch  angle  descent)  which  can  be  easily  flown  or  monitored 
by  the  pilot  of  the  aircraft.  Thus,  the  parametric 
optimal  path  embodies  both  fuel-optimality  and  ease  of 
flying  characteristics. 


TABLE  7 


FUEL  CONSUMPTION  FOR  500-NM  TRAJECTORIES 


Trajectory 

Fuel  Use  (Lb) 

Percent  Over 
Minimum 

Parametric  Optimum 

15,418 

1.1 

Steepest-Descent  Optimum 

15,248 

- 
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The  effects  of  increasing  range  on  optimal  flight 
paths  are  demonstrated  in  Fig.  26.  Here,  altitude-vs . - 
Mach  number  profiles  for  the  200-nm  and  the  500-nm  fully 
optimal  (steepest-descent ) trajectories  are  compared. 
Although  500  nm  is  classified  as  a short  trajectory  range, 
this  profile  demonstrates  the  same  characteristics  as  longer 
range, fully  optimal  profiles.  The  climb  schedule  consists 
of  a nearly  constant- IAS  path  followed  by  a constant-Mach 
number  phase.  Again,  the  descent  profile  is  flown  with 
the  vehicle  near  its  maximum  value  of  lift-to-drag  ratio. 

The  200-nm  trajectory  is  characterized  by  a constant-Mach 
number  climb  schedule,  i.e.,  the  constant-IAS  speed  seg- 
ment of  climb  disappears  as  the  trajectory  range  decreases. 


a-MMt 


Figure  26  Short-Range  Altitude-vs.-  Mach  Number 

Profiles  (Steepest-descent  trajectories) 
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Elapsed  flight  time,  in  minutes,  and  values  of  fuel 
consumption  for  the  two  short-range  trajectories  are  indi- 
cated at  various  points  on  the  profiles.  The  two  profiles 
intersect  during  the  climb  modes  at  an  altitude  of  about 
20,000  ft  and  a Mach  number  of  approximately  0.53.  This 
intersection  denotes  the  point  at  which  the  aircraft's  spe- 
cific total  energy  (kinetic  plus  potential  energy  divided  by 
weight)  is  the  same  for  both  flight  paths.  At  this  equal 
energy  point,  the  200-nm  vehicle  has  burned  more  fuel  and 
has  flown  further  than  its  longer-range  counterpart.  This 
occurs  because  the  aircraft  which  is  flying  the  short-range 
trajectory  is  maximizing  the  range  travelled  during  the 
climb.  On  the  other  hand,  the  500-nm  vehicle  maximizes 
energy  gained  during  the  climb  and  maximizes  specific  range 
during  the  cruise.  In  this  manner,  total  fuel  use  for  the 
mission  is  minimized. 

As  trajectory  range  becomes  greater  than  500  nm, 
the  cruising  flight  mode  takes  on  increasing  significance. 
Fuel  consumption  and  optimal  profile  shapes  for  longer 
range  trajectories  are  discussed  in  the  following  sections. 


5.3  MEDIUM-RANGE  FLIGHT  PATHS 

Fuel-minimal  trajectory  profiles  for  medium-range 
(1,000  to  2,000  nm)  flight  paths  are  presented  in  this  sec- 
tion. As  in  the  case  of  the  short-range  trajectories, 
optimal  trajectory  shaping  provides  significant  reductions 
in  aircraft  fuel  consumption. 

As  trajectory  range  increases,  the  duration  of  the 
cruising  segment  of  flight  increases  more  quickly  than 
that  of  the  climb  and  descent  modes.  The  predominance 


85 


of  the  cruise  for  these  longer  path  lengths  is  demonstrated  in 
Fig.  27.  As  path  length  varies  between  500  to  4,000  nm,  the 
percentage  of  the  trajectory  distance  which  is  flown  in  the 
cruise  mode  increases  from  about  28  to  90  percent.  Because  the 
cruise  mode  is  efficient  in  terms  of  the  nautical  miles  flown 
per  pound  of  fuel  consumed  (specific  range),  the  longer  range 
trajectories  result  in  larger  values  of  total  specific  range. 


INITIAL  CONDITIONS 
V.  • Ml  Kf  At 
hs  * 3 JO 00  h 
W0  - MO.OOO  Ik 

TERMINAL  CONDITIONS 

h«  - 1.000  it 

V,  - 104-107  KIAS 


Figure  27  Cruise  Duration  as  a Percentage 
of  Path  Length 

The  parametric-optimum  and  the  steepest-descent- 
optimum  profiles  are  similrr  for  the  medium-range  trajec- 
tories. Plots  of  altitude  versus  range  for  these  two  opti- 
mal profiles  are  presented  in  Fig.  28.  The  profiles  are 
indicative  of  a 1,000-nm  trajectory  for  an  Initial  vehicle 
weight  of  260,000  lb.  The  major  difference  between  the 
trajectories  occurs  during  the  transition  from  the  climb 
to  the  cruise  phase  of  flight.  The  fully  optimal  (steepest- 
descent)  profile  performs  a more  gradual  transition  than 
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Figure  28  Fuel-Minimal  Medium-Range 
Trajectory  Profiles 


transition  than  that  of  the  parametric  optimum.  Also, 
the  fully  optimal  profile  does  not  climb  quite  as  much  dur- 
ing cruise. 

Conventional  trajectory  profiles  include  cruise 
at  constant  altitudes  and  Mach  numbers  due  to  air  traffic  con- 
trol regulations  (Refs.  30  and  32).  Because  an  optimal 
cruise  increases  altitude  as  vehicle  weight  decreases,  the 
conventional  constant-altitude  cruise  is  clearly  sub-optimal. 
The  Mach  numbers  flown  on  conventional  trajectories  also 
differ  from  the  C-141A  optimal  value  of  approximately  0.725. 
Comparisons  of  fuel  consumption  for  the  two  optimal  1.000-nm 
profiles  and  four  conventional  trajectories  are  presented 
in  Table  8.  In  addition  to  flying  at  constant  altitude 
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during  cruise,  the  conventional  flight  paths  also  include 
the  handbook  climb  and  descent  schedules  described  in 
Sections  4.3.1  and  4.4.1. 

TABLE  8 

FUEL  CONSUMPTION  FOR  1,000-NM 
TRAJECTORIES 


Trajectory  Type 

1,000-nm  Range 
260, 000- lb  Initial  Weight 

Fuel  Use 
(Lb) 

Percent  Over 
Opt imum 

Steepest -Descent  Optimum 

26,932 

- 

Parametric  Optimum 

27,274 

1.3 

Conventional 

Altitude  (h) 

Mach  Number  (M) 

31,000  ft 

0.740 

30,032 

11.5 

31,000  ft 

0.767 

31,199 

15.8 

35,000  ft 

0.740 

28,452 

5.G 

35,000  ft 

0.767 

29,233 

8.5 

Constant-Altitude  Optimum 

h * 39,114  ft 

27,293 

1.3 

Also  included  in  the  above  table  is  the  fuel  con- 
sumption required  for  a constant-altitude  optimum"  profile. 
This  profile  follows  the  same  climb  and  descent  schedule  as 
the  parametric  optimum,  and  M ■ 0.725,  but  the  cruise  mode  is 
flown  at  a constant  altitude.  The  altitude  is  chosen  to 
be  the  average  value  of  the  Initial  and  final  cruise  altitudes 
experienced  by  the  parametrically  defined  optimum  cruise. 

The  parametric -optimum  flight  path  provides  for  fuel 
consumption  within  1.3  percent  of  the  fully  optimal  trajectory 
This  performance  is  significantly  superior  to  the  results  ob- 
tained with  the  handbook  profiles  (5.6  to  15.8  percent  above 
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optimum).  The  ’’constant-altitude-optimum"  fuel  consumption 
is  basically  the  same  as  that  of  the  parametric  optimum. 

Since  constant  altitude  cruise  is  easily  specified  and  is 
the  method  currently  employed,  it  is  reasonable  to  suggest 
this  profile  as  a convenient  means  of  obtaining  substantial 
fuel  savings. 

Table  9 presents  fuel  consumption  values  for  tra- 
jectories with  path  lengths  of  2,000  nm.  Here,  the  param- 
etric-optimum profile,  which  has  demonstrated  fuel  consumption 
values  extremely  close  to  the  fully  optimal  profiles,  is 
chosen  as  the  point  for  comparison.  As  with  the  1,000-nm  case, 
the  parametric- and  constant-altitude-optimum  trajectories 
yield  fuel  consumption  values  significantly  lower  than  these 
of  the  conventional  profiles. 


TABLE  9 

FUEL  CONSUMPTION  FOR  2.000-NM 
TRAJECTORIES 


Trajectory  Type 
2,000-nm  Range 
260,000-lb  Initial  Weight 

Fuel  Use 
(Lb) 

Percent  Over 
Opt imum 

Parametric  Optimum 

49,291 

— 

Conventional 

Altitude  (h) 

Mach  Number  (M) 

31,000  ft 

0.740 

56,184 

14.0 

31,000  ft 

0.767 

58,918 

19.5 

35,000  ft 

0.740 

52,223 

5.9 

35,000  ft 

0.767 

54,153 

9.9 

Constant-Altitude  Optimum 

h « 40,147  ft 

49,539 

0.5 
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The  fuel  consumption  data  presented  in  Tables  8 and 
9 demonstrate  the  utility  of  trajectory  optimization  in 
reducing  fuel  consumption  for  varying  path  length  missions. 

In  all  cases,  however,  the  initial  vehicle  weight  is  260,000 
lb.  The  effects  of  differences  in  initial  weight  on  fuel 
consumption  are  shown  in  Fig.  29.  In  this  figure',  fuel  con- 
sumption is  plotted  for  a 2,000-nm  trajectory  flown  by  air- 
craft with  initial  weights  of  200,000,  260,000  and  320,000  lb 
For  each  vehicle  weight,  the  fuel  consumption  represents  the 
amount  of  fuel  used  in  flying  the  parametrically-defined 
optimal  profiles.  The  cruise  altitudes  for  each  of  the  pro- 
files also  are  indicated  on  the  figure.  Over  the  range  of 
initial  vehicle  weights  examined,  the  increase  in  fuel  for 
each  pound  of  initial  weight  is  linear,  with  a slope  of  0.19 
lb-fuel/ lb-initial  weight.  The  penalty  paid  for  carrying 
excess  fuel  on  this  2,000-nm  trajectory  is  equal,  therefore, 
to  19  percent  of  the  extra  fuel.  This  points  up  the  desira- 
bility of  not  loading  the  vehicle  with  more  fuel  than  is 
required  to  carry  out  the  mission  with  a reasonable  degree 
of  safety. 


Figure  29  Effects  of  Initial  Aircraft  Weight  on  Total  1’uel 
Consumption  ( Parametric -opt imum  profiles) 
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5.4 


LONG-RANGE  FLIGHT  PATHS 


Potential  fuel  reduction  for  a long-range  mission 
is  even  greater  than  that  of  the  short-  and  medium-range 
flight  paths.  In  this  study,  trajectory  ranges  of  more 
than  2.000  nm  are  categorized  as  *’ long-range . " 

Many  of  the  results  previously  presented  here  are 
determined  for  an  intermediate  value  of  initial  C-141A  air- 
craft weight  (260,000  lb).  The  aircraft  weight  and  maximum 
fuel-loading  data  for  a C-141A  aircraft  indicate  that  at  a 
vehicle  weight  of  316,000  lb,  the  maximum  amount  of  fuel 
which  can  be  loaded  for  flight  is  148,000  lb.  Thus,  for  an 
initial  weight  of  260,000  lb,  about  92,000  lb  of  fuel  is 
available  to  fly  a mission. 

The  parametric-optimum  4.000-nm  trajectory  con- 
sumes 87,285  lb  of  fuel  for  an  initial  260,000-lh  total 
aircraft  weight.  This  value  is  close  to  the  limit  of  fuel 
which  can  be  loaded  for  this  initial  weight.  Thus,  a 4,000-nm 
range  is  approximately  the  maximum  range  of  a C-141A  with 
this  loading.  Fuel  consumption  values  for  this  parametric- 
optimum,  four  handbook  profiles,  and  the  constant-altitude 
optimum  are  presented  in  Table  10.  Note  that  all  the  hand- 
book profiles  require  more  than  the  92,000  lb  of  fuel  avail- 
able. The  relative  percentage  fuel  savings  indicated  in  the 
table  are,  however,  indicative  of  the  potential  fuel  reduc- 
tions which  can  be  achieved  with  trajectories  of  this  range. 
Also,  the  constant-altitude-optimum  profile  again  provides 
performance  which  is  nearly  equal  to  that  of  the  parametric- 
optimum  profile. 

The  maximum  range  capability  of  the  modeled  C-141A 
aircraft  (see  Appendix  B),  is  5,600  nm.  This  long-range 
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TABLE  10 


FUEL  CONSUMPTION  FOR  4,000-NM 
TRAJECTORIES 


Trajectory  Type 
4,000-nm  Range 
260. 000-lb  Initial  Weight 

Fuel  Use 
(Lb) 

Percent  Ov»r 
Optimum 

Parametric  Opt  imam 

87,285 

— 

Conventional 

Altitude  (h) 

Mach  Number  (M) 

31,000  ft 

0.740 

105,850 

21.3 

31,000  ft 

0.767 

111,605 

27.9 

35,000  ft 

0.740 

96,219 

10.2 

35,000  ft 

0.767 

100,511 

15.2 

Constant-Altitude  Optimum 

h - 42,215  ft 

88,636 

1.5 

profile,  Fig.  30,  has  been  generated  by  performing  a param- 
etric optimization  for  a vehicle  with  an  initial  weight  of 

320.000  lb.  The  aircraft  burns  138,328  lb  of  fuel,  only 

10.000  lb  less  than  the  maximum  which  can  be  loaded,  on  this 
maximum  range  mission. 

The  important  characteristics  of  long-range  fuel- 
optimal  trajectory  profiles  are  highlighted  by  the  5,600-nm 
flight  path.  About  9 percent  (500  nm)  of  the  trajectory 
is  flown  during  the  climb  and  descent  segments.  Thus, 
by  optimizing  fuel  consumption  during  cruise,  significant 
reductions  in  fuel  use  are  achieved.  The  increasing  al- 
titude characteristic  of  the  optimal,  constant-W/6  climb 
also  is  clearly  demonstrated  on  this  trajectory.  The 
initial  cruise  altitude  is  approximately  34,500  ft,  and 
the  cruise  phase  terminates  at  an  altitude  of  about 
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n-»04*o 


Figure  30  5,600-nm  Optimal  Trajectory  Profile 

44,500  ft.  While  this  10,000-ft  increase  in  altitude  during 
cruise  is  larger  than  that  which  would  be  experienced  on 
shorter  length  profiles,  the  climbing-cruise  characteristic 
of  the  optimal  profiles  is  the  same. 


5.5  FUEL/TIME/COST  TRADEOFFS 

The  trajectory  shaping  efforts  of  this  study  are 
directed  toward  minimizing  aircra't  fuel  consumption.  By 
lowering  the  quantity  of  fuel  used  on  any  given  flight, 
the  dollar  costs  associated  with  fuel  are  reduced.  Many 
aircraft  mission  costs  can  be  related  to  the  time  duration 
of  the  flight.  Personnel  costs,  frequency  of  aircraft 
overhaul,  and  equipment  life  expectancy  are  among  those 
factors  which  contribute  to  time-related  costs.  It  is 
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demonstrated  in  this  section  that  trajectories  which  are 
fuel-optimal  may  not  be  optimal  in  terms  of  minimizing 
total  costs;  however,  as  fuel  costs  increase,  the  minimum 
fuel  trajectories  do  provide  the  lowest  total  costs. 

A major  factor  influencing  the  interest  in  fuel- 
use  reduction  is  the  current  dramatic  price  increase  for 
petroleum  products.  USAF  total  aircraft  fuel  consumption 
and  fuel  (JP4)  costs  per  gallon  for  the  years  1972  through 
1976,  as  obtained  from  Ref.  25,  are  shown  in  Fig.  31. 

It  is  sean  here  that  while  fuel  consumption  is  expected  to 
remain  relatively  constant,  fuel  costs  are  expected  to  be 
about  four  times  higher  in  1976  than  they  were  in  1973.  As 
the  rise  in  fuel  costs  is  more  rapid  than  the  increase  in 
time-related  costs,  fuel  reduction  becomes  an  increasingly 
important  objective  in  reducing  total  mission  expenses. 


Figure  31  Yearly  Fuel  Consumption 
and  Fuel  Costs 
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The  tradeoff  between  fuel-based  costs  and  time- 
based  costs  is  exemplified  by  the  comparison  of  the  200-nm 
parametric-optimum  and  conventional  trajectories.  The  fuel 
consumed  and  the  elapsed  time  for  these  trajectories  are  pre- 
sented in  Table  11.  Also  shown  is  the  total  cost  (in  dollars) 
for  each  trajectory,  assuming  a fuel  cost  of  45£/gallon, 

6.55  lb  per  gallon,  and  a time-based  cost  of  $400/hr.  The 
minimum-fuel  parametric-optimum  trajectory  uses  932  lb  less 
fuel  than  the  conventional  trajectory,  but  its  flight  time 
is  0.11  hr  longer.  This  results  in  a nearly  equal  value  of 
total  cost  for  the  two  trajectories. 


TABLE  11 

TOTAL  TRAJECTORY  COSTS  — 
200-NM  FLIGHT  PATH 


Trajectory  Type 

Fuel  Use  (Lb) 

Time  (Hr) 

Total  Cost  ($) 

Parametric  Optimum 

7,515 

0.664 

782 

Conventional  (1/3-Level) 

8,347 

0.554 

795 

The  effects  of  increased  fuel  prices  on  the  cost- 
optimality  of  minimum-fuel  trajectories  are  demonstrated  in 
Table  12.  Here,  total  costs  for  the  200-nm  parametric- 
optimum  and  conventional  trajectories  are  shown  for  a vari- 
ety of  fuel  and  time-based  costs.  With  an  increase  in  fuel 
prices,  the  optimal  trajectory  provides  minimum  total  cost 
as  well  as  minimum  fuel  performance. 


TABLE  12 

TRAJECTORY  COSTS  VS.  FUEL/TIME  PRICE  COMBINATIONS  — 

200-NM  PATH  LENGTH 


Trajectory  Type 

Cost  I 
* ($) 

Cost  II 
($) 

Cost  III 
($) 

Cost  IV 
($) 

Parametric  Optimum 

782 

1,126 

1,413 

1,545 

Convent ional 

795 

1,178 

1,496 

1,606 

• COST  I 45(7  gal  Ion,  $400/hr  • COST  II  75^/Rallon,  $40P/hr 

• COST  III  $1 . 00/fial Ion , $400/hr  • COST  IV  $1 . 00/gal Ion , $S00/hr 
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Trajectory  fuel  consumption,  time  duration  and 
total  costs  for  1,000-,  2,000-,  and  4,000-nm  missions  are 
presented  in  Table  13.  As  before,  the  total  costs  are  based 
on  a rate  of  45^/gallon  for  fuel  and  $400/hr  for  time.  Total 


trajectory  expense  for  other  values  of  fuel  and  time  costs 
can  be  computed  from  the  fuel-consumption  and  time  data  con- 


tained ir  the  table.  For  each  trajectory  range,  parametric- 


optimum.  constant-altitude  optimum,  and  four  handbook  tra- 
jectory results  are  shown.  The  handbook  trajectories  are 
denoted  as  HB1  through  HB4,  and  the  constant  altitudes  (h) 
and  constant  Mach  numbers  (M)  which  describe  their  cruise 
conditions  are  as  follows: 


• HB1 

• HB2 

• HB3 

• HB4 


h » 31,000  ft,  M 
h - 35,000  ft,  M 
h - 31,000  ft,  M 
h = 35,000  ft,  M 


0.74 

0.74 

0.767 

0.767 


TABLE  13 


TRAJECTORY  COST  AS  A FUNCTION  OF  PATH 
LENGTH  AND  FLIGHT  PROFILE 


Trajectory 

Range 

Trajectory  Type 

Fuel 

(Lb) 

Time 

(Hr) 

Total 

Cost 

($) 

Excess 

Cost 

(Percent) 

t 

Parametric  Optimum 

27 , 274 

2.671 

2,943 

- 

Constant-Altitude  Optimum 

27,293 

2.671 

2,943 

0.0 

HD1 

30,032 

2.400 

3,023 

2.7 

HB2 

28,452 

2.434 

2,927 

-0.5 

t 

HB3 

31,199 

2.336 

3,078 

4.6 

HB4 

29,233 

2.358 

2,951 

0.3 

•t 

Parametric  Optimum 

49,291 

5.104 

5,428 

— 

Constant-Altitude  Optimum 

49,539 

5. 104 

5,445 

0.3 

HB1 

56,184 

4.701 

5,741 

5.8 

HB2 

52,223 

4.790 

5,504 

1.4 

t 

HB3 

58,918 

4.554 

5,876 

8.3 

HB4 

54,153 

4.635 

5,574 

2.7 

Parametric  Optimum 

87,285 

9.970 

9,987 

_ 

Constant-Altitude  Optimum 

88,636 

9.978 

10,080 

0.9 

HB1 

105,850 

9.286 

10,986 

10.0 

HB2 

96,219 

9.453 

10,391 

4.0 

♦ 

HB3 

111,605 

8.988 

11,262 

12.8 

HB4 

100,511 

9.147 

10,564 

5.8 

*Lowest  Total  Cost 
tMinimum  Fuel 
tUlnimum  Tine 
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For  all  trajectory  ranges,  the  parametric -optimum 
profile  provides  the  lowest  value  of  fuel  consumption.  In 
the  2,000-  and  4,000-ntn  profiles,  the  parametric  optimum 
also  results  in  the  lowest  trajectory  costs.  As  fuel  prices 
continue  to  rise,  the  savings  available  from  these  trajec- 
tories will  become  even  more  pronounced. 


In  the  case  of  the  1,000-nm  trajectory,  the  lowest 
cost  is  obtained  with  a conventional  profile,  HB2.  Although 
this  profile  consumes  1,178  lb  more  fuel  than  does  the  param- 
etric optimum,  the  flight  time  is  reduced  by  0.237  hr.  At  a 
rate  of  $400/hr,  this  is  enough  to  offset  the  increased  fuel 
costs.  Trajectory  costs  for  other  fuel  price  combinations 
are  presented  in  Table  14.  As  fuel  prices  increase,  the 
minimum-fuel  trajectories  also  are  the  minimum-total-cost 
profiles. 


TABLE  14 


TRAJECTORY  COSTS  FOR  A VARIETY  OF 
FUEL/TIME  PRICES 


Trajectory 

Range 

Trajectory  Type 

Cost  I 
($) 

Cost  II 
($) 

Cost  III 
($) 

Cost  IV 
($) 

Parametric  Optimum 

2,943 

■RS 

5,233 

5,767 

Constant  Altitude  Optimum 

2,943 

5,234 

5,770 

1 , 000  nm 

HB1 

3,023 

xWf'i  ^ 

5,545 

6,025 

HB2 

2,927 

-stfrwl 

5,317 

5,804 

HB3 

3,078 

5,698 

6,165 

HB4 

2,951 

IB 

5,406 

5,878 

Parametric  Optimum 

5,428 

7,686 

9,567 

10.587 

Constant-Altitude  Optimum 

5,445 

7,714 

9,605 

10,625 

2,000  nm 

HB1 

5,741 

8,314 

10,459 

11,399 

HB2 

5,504 

7,896 

9,889 

10,847 

HB3 

5,876 

8,574 

10,823 

11,727 

HB4 

5,574 

8,055 

10,122 

11,049 

Parametric  Optimum 

9,987 

13,985 

17,317 

19,307 

Constant-Altitude  Optimum 

10,080 

14,139 

17,523 

19,519 

4 , 000  nm 

HB1 

10,986 

15,834 

19,874 

21,732 

HB2 

10,391 

14,798 

18,471 

20,362 

HB3 

11,202 

16,374 

20,634 

22,432 

HB4 

10,564 

15,168 

19,004 

20,833 

• COST  I 45£/gallon,  $400/hr  • COST  II  75<7gallon,  $400/hr 

• COST  III  $1. 00/gal Ion,  $400/hr  • COST  IV  $1 .00/gallon , $600/hr 
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Tables  13  and  14  demonstrate  fuel-use  trends  for 


vehicles  whose  init ial  total  weights  are  equal.  The  results 
of  comparing  trajectories  with  equal  f inal  weights,  reflect- 
ing decreased  "tankering"  of  excess  fuel , is  even  more  dra- 
matic. Assuming  that  the  weight  ratio,  u,  defined  by 


w 


Wn-AW,  ^ 
0 fuel 


(5.5-1) 


is  unchanging  for  each  mission  length  and  trajectory  type, 
the  fuel  savings  obtained  by  flying  parametric-optimum  tra- 
jectories instead  of  minimum-time  conventional  trajectories 
are  compared,  with  the  following  results: 

• 1,000-nm  trajectory: 

Fuel  saved  = 4,419  lb  (14  percent) 

• 2,000-nm  trajectory: 

Fuel  saved  = 11,865  lb  (20  percent) 


• 4,000-nm  trajectory: 

Fuel  saved  * 36,666  lb  (33  percent) 

In  this  case,  the  fuel  saved  is  fuel  which  is  not  loaded  on 
the  aircraft,  so  the  initial  weights  of  the  aircraft  which 
fly  the  parametric-optimum  trajectories  are  223,334  to 
255,571  lb.  The  minimum-time  aircraft  begins  at  260,000  lb 
in  each  case. 


Another  interesting  comparison  is  payload  weight 
improvement  for  fixed  initial  weight  (260,000  lb).  Compar- 
ing the  parametric-optimum  profile  to  the  minimum-fuel  con- 
ventional profile  indicates  the  following  percentage  increases 
in  available  payload  weight,  assuming  10  percent  fuel  reserve: 

• 1,000-nm  trajectory:  1.5  percent 

• 2,000-nm  trajectory:  5.2  percent 

• 4,000-nm  trajectory:  69.4  percent 
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The  latter  figure  is  seen  to  have  significant  implications 
for  the  strategic  airlift  mission,  because  it  says  that  the 
same  payload  could  be  moved  with  41  percent  fewer  flights. 
Since  the  fuel  used  on  each  flight  also  is  9 percent  less, 
equal  total  payload  weight  would  be  moved  with  46  percent 
less  fuel. 

On  an  actual  strategic  mission,  the  aircraft  would 
operate  closer  to  its  maximum  gross  weight  on  takeoff,  carry- 
ing more  payload  and  fuel  in  both  the  optimum  and  conventional 
cases,  and  the  above  results  do  not  account  for  volumetric 
limitations  on  specific  payloads.  Consequently,  the  percen- 
tage payload  gain  could  be  reduced.  The  gain  would  be  fur- 
ther reduced  by  the  fact  that  the  average  altitude  of  the 
optimal  trajectory  would  be  lower  for  higher  weights  and, 
therefore,  nearer  the  HB2  flight  path.  Nevertheless,  the 
260,000-lb  result  indicates  the  proper  trends  for  the  bene- 
fits obtained  from  fuel-optimal  flight. 

Rather  than  determining  trajectory  profiles  to 
minimize  fuel  consumption,  the  analysis  tools  employed  in 
this  study  could  be  used  to  determine  cost -optimal  profiles. 
Determination  of  a cost-efficient  cruise  phase  might  proceed 
in  a fashion  similar  to  that  of  Section  4.2.  However, 
instead  of  maximizing  specific-range  (nautical  miles 
travelled  per  pound  of  fuel),  a specific-cost  factor 
(nautical  miles  per  dollar)  would  be  substituted.  The 
dollar  costs  associated  with  each  nautical  mile  of  path 
length  are  the  price  of  the  fuel  consumed  plus  the  costs 
associated  with  the  time  duration  for  flying  each  nautical 
mile. 


Fully  optimal  minimum-cost  trajectories  can  be 
determined  by  variational  means,  as  described  in  Section  3.6. 
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The  cost  function  of  Eq.  (3.6-1)  need  simply  be  modified  to 
include  a terminal  penalty  on  time  as  well  as  on  fuel  con- 
sumption. Thus,  the  present  fuel-optimal  program  develop- 
ment (see  Appendix  A)  also  is  capable  of  handling  the 
minimum-cost  trajectory  problem. 


5 . G CHAPTER  SUMMARY 

This  chapter  has  demonstrated  the  utility  of  flying 
optimized  (fuel-minimal)  trajectory  profiles  in  order  to 
reduce  total  fuel  consumption.  The  optimized  flight  seg- 
ments described  in  Chapter  4 serve  as  the  starting  point  for 
the  current  complete  trajectory  optimization.  By  combining 
these  piecewise  optimal  segments,  parametric-optimum  pro- 
files are  constructed.  These  profiles  are  found  to  provide 
excellent  fuel  economy  when  compared  with  the  fully  optimal 
(steepest-descent)  flight  paths. 

The  potential  for  fuel  use  reduction  by  means  of 
trajectory  shaping  has  been  shown  to  exist  for  short-,  medium-, 
and  long-range  trajectories.  The  optimal  short-range  tra- 
jectories demonstrate  a maximum  fuel  savings  of  about  12  per- 
cent. This  percentage  increases  (with  increasing  path  length) 
to  approximately  28  percent  for  the  long-range  flight  paths. 

Constant-altitude-optimal  flight  paths  also  have 
demonstrated  fuel  efficiency  close  to  that  of  the  fully 
optimal  trajectories.  These  constant-altitude-optimal  pro- 
files are  characterized  by  parametrically  defined  climb  and 
descent  segments,  while  the  cruise  phase  is  flown  at  a con- 
stant altitude.  The  altitude  is  chosrtn  to  be  the  midpoint 
value  defined  by  a fuel-minimal,  constant  W/6  flight  path. 
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The  final  point  emphasized  in  this  chapter  is  that 
fuel-minimal  profiles  are  not  necessarily  cost -minimal  pro- 
files. As  cost 8 associated  with  flight  time  are  factored 
into  the  problem,  it  is  found  that  some  conventional  tra- 
jectories are  competitive,  in  terms  of  total  trajectory 
cost,  with  the  fuel-optimal  flight  paths.  As  fuel  costs 
increase  at  a dramatic  rate  (they  have  quadrupled  in  four 
years),  the  fuel-optimal  characteristics  of  mission  pro- 
files become  increasingly  significant.  Simple  methods  are 
presented  which  can  be  employed  to  minimize  total  (fuel 
plus  time)  costs.  These  methods  are  similar  to  the  techniques 
employed  for  the  fuel-minimization  study  and  are,  therefore, 
expected  to  provide  good  results  with  a minimum  of  effort. 

Having  generated  the  fuel-optimal  trajectory  pro- 
files for  a variety  of  flight  segments  and  mission  con- 
straints, the  problem  remaining  is  to  develop  a means  of  fly- 
ing these  profiles.  An  automated  method  for  profile  follow- 
ing which  can  be  implemented  on  actual  aircraft  is  presented 
in  the  following  chapter. 
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VI 


SYNTHESIS  OF  FUEL-OPTIMAL  GUIDANCE  ALGORITHMS 


6.1  OVERVIEW 

The  earlier  chapters  have  established  the  ground- 
work for  fuel-optimal  guidance  by  discussing  use  factors, 
presenting  analytical  tools,  and  developing  a wide  range 
of  fuel-optimal  flight  paths.  This  chapter  presents  the 
means  of  using  these  results  to  define  throttle/energy 
management  algorithms  which  can  be  implemented  in  an 
individual  aircraft.  As  introduced  in  Chapter  1,  the  basic 
approach  is  to 

• Define  fuel-optimal  flight  paths 

• Compute  linear ^optimal  guidance  gain 
schedules 

• Design  a real-time  system  for 
generating  the  reference  flight  paths 
and  following  this  path  in  the  presence 
of  disturbances. 

Section  6.2  describes  the  separation  of  nonlinear 
and  linear  dynamic  models  which  allows  this  approach  to 
throttle/energy  management  to  be  used,  and  it  presents 
the  linear  perturbation  guidance  law  which  is  required 
for  flight  path  following.  The  fuel-optimal  flight  path 
results  of  earlier  chapters  are  generalized  for  on-board 
implementation  in  Section  6.3,  and  corresponding  func- 
tional requirements  on  sub-system  components  are  discussed 
in  Section  6.4.  The  chapter  is  summarized  in  Section  6.5. 
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6.2 


LINEAR-OPTIMAL  GUIDANCE  LAVS 


In  this  section,  the  formulation  of  a guidance  sys- 
tem designed  to  follow  a reference  trajectory  is  presented. 
Although  the  concepts  are  applicable  to  all  trajectories, 
the  current  efforts  are  directed  toward  flying  the  fuel- 
optimal  profiles  previously  described.  The  resulting  algo- 
rithm provides  a guidance  technique  which  minimizes  aircraft 
fuel  consumption  when  used  in  conjunction  with  the  fuel- 
optimal  flight  path  generator  described  in  the  next  section. 


Linear-optimal  guidance  algorithms  (Refs.  38  and  39) 
compute  guidance  gains  to  yield  desired  state  deviation  and 
control  perturbation  response.  This  is  accomplished  by 
determining  guidance  comnands  to  minimize  the  following 
quadratic  cost  function: 


J - AxT(Rf)Pf4x(Rf) 


rRt  T 

v / { Ax  AAx 


♦ Au  BAuJdR 


(6.2-1) 


where  the  state  and  control  vector  deviations  are  defined 
by  Ax(R)  and  Au(RX  respectively,  and  the  matrices, A and  B, 
are  used  to  penalize  in-flight  deviations  in  state  and  con- 
trol. The  matrix,  Pf,  is  used  to  penalize  the  terminal  tra- 
jectory deviations,  and  the  independent  variable,  R,  is 
trajectory  path  length.  By  specifying  large  values  for  the 
elements  of  A,  the  cost  function  places  emphasis  on  main- 
taining small  state  deviations.  Large  values  of  B,  rela- 
tive to  A,  mean  that  the  control  perturbations  will  remain 
small,  even  when  faced  with  large  state  deviations.  The 
choice  of  values  used  in  the  current  study  is  discussed  in 
Section  6.2.2. 
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If  the  dynamic  relationship  between  the  state  and 
control  is  linear,  i.e. , 

Ax'(R)  - F(R)Ax(R)  ♦ G(R)Au(R)  (6.2-2) 

then  the  control  history  which  minimizes  Eq.  (6.2-1)  can  be 
specified  in  a feedback  fashion  as 

Au(R)  - -K(R)Ax(R)  (6.2-3) 

(The  notation  ()'  represents  differentiation  with  respect 
to  path  length). 

In  the  above  equations,  the  matrices  F(R)  and  G(R) 
describe  the  effects  of  state  and  control  deviations  on  the 
state  perturbation  derivative,  Ax'(R),  and  the  elements  of 
K(R)  are  the  desired  perturbation  feedback  gains. 

The  control  perturbations,  Au(R),  of  Eq.  (6.2-3)  are 
linear  functions  of  the  state  deviations,  Ax(R),  and  they  pro- 
vide the  minimum  (opt imal)  cost  for  Eq.  (6.2-1).  Thus,  the 
descriptor  "linear-optimal"  truly  applies  to  the  guidance 
algorithm,  as  long  as  perturbations  from  the  reference  flight 
path  are  described  adequately  by  Eq.  (6.2-2).  The  method  of 
determining  these  guidance  gains  is  presented  in  the  follow- 
ing section. 

6.2.1  Design  Fundamentals 

The  aircraft  equations  of  motion  derived  in 
Section  3.2  are  nonlinear  functions  of  the  state  and  control 
variables.  The  trajectory  perturbations  can  be  expressed 
in  the  linear  form  of  Eq.  (6.2-2)  by  performing  a Taylor 
series  expansion  of  the  exact  equations  and  retaining  only 
first-order  terms. 
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The  perturbed  trajectory  equations  can  be  derived  by 
substitution  In  Eq.  (3.6-2): 

(Xjj+Ax)'  - f ICXjj+AxJ.CUp+Au)!  (6.2-4) 

where  x^  and  u^  are  the  state  and  control  vectors  which  de- 
scribe the  desired  flight  path,  and  Ax  and  Au  represent  tra- 
jectory perturbations.  The  first-order  Taylor  series  expan- 
sion of  Eq.  (6.2-4)  results  In 


*6  + 4*d  * 1(VV  + £x 


Ax  ♦ f Au 
D D 


(6.2-5) 


where  fx  and  fu  are  the  motion  equation  partial  derivative 
matrices  with  respect  to  the  state  and  control  variables,  re- 
spectively, evaluated  along  the  desired  trajectory. 


By  definition,  the  desired  trajectory  satisfies 
the  nonlinear  equations  of  motion,  l.e., 

5d  ” (6.2-6) 


Subtracting  Eq.  (6.2-6)  from  Eq.  (6.2-5),  It  can 
be  seen  that  the  dynamics  of  small  variations  from  the 
nominal  flight  path  ar*  described  approximately  by 

Ax’  - FAx  ♦ GAu  (6.2-7) 

where  the  partial  derivative  matrices,  and  _fu,  have  been 
denoted  simply  as  F and  G.  Equation  (672-7)  has  the  required 
linear  form  of  Eq.  (6.2-2). 

The  determination  of  a perturbation  control  his- 
tory to  minimize  Eq.  (6.2-1)  (Refs.  38  and  39)  is  accomplished 
by  integrating  the  following  matrix  Riccati  equation  from 

R#  to  R : 
z o 
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p* 


(6.2-8) 


- -PP  - PTP  ♦ PGB_1GTP  - A 
subject  to  the  boundary  condition 

P(Rf)  - Pf  (6.2-9) 

Solving  this  Matrix  equation  is  analogous  to 
determining  the  costate  variables  in  the  variational  opti- 
mization problem  described  in  Section  3.6.  The  optimal 
perturbation  guidance  gains  then  are  defined  in  terms  of  the 
Riccati  solution  matrix,  P,  the  control  partial  derivative 
matrix,  G,  and  the  control  penalty  matrix,  B,  as  follows: 

K(R)  - (B(R)]"1  GT(R)P(R)  (6.2-10) 

The  solution  of  Eq.  (6.2-8)  by  means  of  numerical 
integration  techniques  is  a computationally  expensive  task. 

The  matrix  Riccati  equation  used  to  compute  guidance  gains 
accounts  for  time  variations  in  aircraft  dynamics  as  fuel  is 
burned  and  as  altitude  and  velocity  are  varied.  These 
dynamic  variations  are  slow  compared  to  the  natural  fre- 
quencies of  the  system.  Small  time  steps  are  required  in  the 
Riccati  differential  equation  solution  as  a consequence  of 
these  natural  frequencies,  and  this  results  in  an  excessive 
number  of  computations  to  analyze  the  gains  for  a single 
mission. 


In  order  to  surmount  these  numerical  difficulties, 
use  is  made  of  the  fact  that  the  partial  derivative 
matrices,  F and  G,  are  slowly  varying  along  the  aircraft 
trajectories.  Rather  than  integrating  Eq.  (6.2-8)  back- 
wards along  the  complete  trajectory  profile,  only  steady- 
state  gain  values  (for  constant  values  of  F and  G)  at  discrete 
trajectory  points  (e.g.,  every  50  nm)  are  computed.  Because 
of  the  slowly  varying  nature  of  these  gains,  the  resulting 
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steady-state  gains  are  essentially  identical  to  those 
obtained  from  the  complete  numerical  integration  of  the 
Riccati  equation. 


The  computation  of  steady-state  gains  is  signifi- 
cantly simpler  than  the  determination  of  the  time-varying 
values.  Transition  matrix  solutions  (the  "Kalman-Englar 
Method")  are  used  to  replace  the  numerical  integration 
methods  (Ref.  54).  Because  the  partial  derivative  matrices 
(F  and  G)  are  assumed  to  be  constant  at  each  steady-state 
gain  computation  point,  the  necessary  transition  matrix 
need  be  calculated  only  once  at  each  of  these  points.  The 
Kalman-Englar  method  employed  in  this  study  reduced  the  com 
putational  effort  required  for  direct  numerical  integration 
gain  evaluations  by  a factor  of  about  one  hundred. 


6.2.2  Weighting  Matrices  for  System  Design 


The  perturbation  guidance  technique  presented  in 
this  chapter  employs  guidance  gains  which  yield  desired 
state  deviation  and  control  perturbation  response.  These 
gains  are  computed  to  minimize  the  following  scalar  cost 
f unct ion : 


J 


(AxTAAx  ♦ AuTBAu)dR 


(6.2-11) 


This  quadratic  cost  function  is  essentially  the  same  as  that 
of  Eq.  (6.2-1).  The  present  version  does  not  include  a ter- 
minal penalty  (the  term  outside  the  integral  in  Eq.  (6.2-1)), 
as  only  the  steady-state  gain  values  are  of  interest. 
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In  order  to  insure  that  the  cost  function  has  a 
minimum,  the  penalty  weighting  matrices.  A and  B,  must  be 
positive  semi-definite  and  positive  definite, respectively 
(Ref.  54).  This  means  that  for  any  state  or  control  devia- 
tions, the  following  relationships  must  hold: 

AxTAAx  * 0 (6.2-12) 

AuTB Au  > 0 (6.2-13) 

A simple  technique  for  insuring  positive  definite- 
ness is  to  choose  the  A and  B matrices  to  be  of  diagonal 
form  with  positive  or  zero  elements  along  the  diagonal. 

For  a diagonal  matrix,  Eq.  (6.2-12)  can  be  evaluated  as 

AxTAAx  - a^  AxJ  + a22  Ax2  + ...  anQ  Axjj  (6.2-14) 

where  the  quantities,  a^,  are  the  diagonal  elements  of  the  A 
matrix,  and  the  terms,  Ax^  are  the  elements  of  the  state  per- 
turbation vector.  Because  of  its  quadratic  form,  Eq.  (6.2-14) 
is  guaranteed  to  be  greater  than  or  equal  to  zero  as  long 
as  the  coefficients,  a^,  are  greater  than  or  equal  to  zero. 

If  the  values  of  the  elements  of  the  A matrix  are 
chosen  to  be  large  compared  to  those  of  B,  the  state 
deviations  are  penalized  more  heavily  than  are  the  control 
perturbations.  This,  in  turn,  leads  to  large  guidance  gain 
values,  whereby  large  values  of  control  would  be  used  to 
null  even  small  deviations  from  the  reference-trajectory 
states.  Conversely,  if  the  B matrix  elements  are  chosen  to 
be  large  compared  to  those  of  A,  the  control  variables  are 
more  heavily  penalized,  and  only  small  control  perturbations 
will  be  employed  to  bring  the  vehicle  state  close  to  the 
reference  values. 
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The  choice  of  values  for  A and  B In  the  current 
fuel  optimization  etudy  were  derived  from  Information  con- 
cerning deelred  trajectory  state  and  control  response 
obtained  from  U8A7  personnel  (Ref.  2S).  This  information 
was  provided  In  terms  of  the  maximum  trajectory  deviations 
that  one  could  allow  In  an  actual  flight  and  the  typical 
control  perturbation  which  would  be  exercised  to  null  them. 
The  current  analyses  considered  velocity,  AV,  altitude,  Ah, 
and  flight  path  angle,  Ay,  perturbations  as  the  state 
variables  to  be  controlled.  Ay  corresponds  to  an  altitude 
rate  perturbation,  since  Ay  * A&/V.  Pitch  angle,  A6,  and 
throttle  setting,  AT,  are  the  available  perturbation  control 
variables.  The  maximum  deelred  values  of  these  quantities 
are  presented  in  Table  IS 

TABLE  IS 

DESIRED  PERTURBATION  GUIDANCE  CHARACTERISTICS 


Maximum  State  Deviations 

Maximum  Control  Perturbations 

Velocity  (AV):  10  kt  (17fps) 
Flight  Path  Angle  (Ay):  0.17  deg 
Altitude  (^h):  100  ft 

Pitch  Angle  (A6):  1.5  deg 

Throttle  Setting  (AT): 

10  percent 

Specifications  of  A and  B penalty  matrices  to 
reflect  the  maximum  trajectory  deviations  of  Table  6.2-1  is 
accomplished  by  a straightforward  procedure.  The  diagonal 
elements  of  the  penalty  matrices  are  set  equal  to  the  square 
of  the  Inverse  of  the  maximum  allowable  perturbations,  l.e., 

a. .  - 1/ ( Ax±  )2  (6.2-15) 

max 

b. .  - l/(Au.  )2  (6.2-16) 

max 
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where  the  variables  AxifflftX  and  Auimax  are  the  specified 
maximum  allowable  state  and  control  deviations,  respectively. 
This  method  assures  that  the  penalties  associated  with  each 
maximum  trajectory  deviation  contribute  an  equal  amount 
to  the  total  cost  function  of  Eq.  (6.2-11).  For  example,  a 
10-kt  velocity  deviation  causes  the  same  increase  in  the 
cost  function  as  a 100-ft  altitude  deviation  or  a 10-percent 
throttle  perturbation. 


Using  the  values  presented  in  Table  6.2-1  and 
Eqs.  (6.2-15)  and  (6.2-16),  the  A and  B weighting  matrices 
for  the  current  study  are  given  as 


A 


3.46xl0‘3 

0 

0 

0 


0 

l.llxlO5 

0 

0 


0 

0 

l.xlO-4 

0 


0 

0 


•J 


(6.2-17) 


B 


C 


1.445x10' 
0 


° ) 
1.0x10  2J 


(6.2-18) 


In  the  above  equations,  Ay  and  A6  penalties  are  scaled  in 
radians  and  the  AV  penalty  is  scaled  in  fps.  The  fourth 
diagonal  element  of  matrix  A represents  the  penalty  on 
vehicle  mass  deviations.  By  setting  this  element  ^44)  to 
zero,  no  direct  penalty  is  imposed  on  mass  variations.* 
Mass  values  should  be  well  controlled,  however,  as  long  as 
the  actual  trajectory  closely  follows  the  reference  fuel- 


*Note  that  positive  or  negative  mass  variations  are  penalized 
equally  in  Eq.  (6.2-11),  i.e.,  both  cause  J to  increase; 
hence,  the  quadratic  cost  function  cannot  minimize  fuel  use* 
directly — it  can  only  minimize  its  deviation  from  a nominal 
value. 
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optimal  profile.  Dynamic  coupling  in  the  state  equations 
given  by  Eq.  (6.2-6)  results  in  small  gains  on  mass  in  the 
linear  perturbation  guidance  law  even  though  there  is  no 
direct  penalty  function  applied  to  the  mass  deviations. 

6.2.3  Linear-Optimal  Guidance  Gains 

The  preceding  sections  of  this  chapter  have  de- 
scribed the  procedure  by  which  linear-optimal  guidance  gains 
are  computed  and  the  numerical  values  of  the  weighting 
matrices  which  affect  perturbed  trajectory  response.  In 
this  section,  the  actual  guidance  gains  and  simulated  tra- 
jectory responses  are  presented  for  a typical  C-141A  flight 
profile. 


For  the  current  study,  the  vector  perturbation  con- 
trols, Eq.  (6.2-3)  can  be  written  in  component  form  as: 

A0  - k^Vjj-V)  + k12(YD-Y)  + k13(hD-h)  + k^mp-m) 

- eD  - e (6.2-19) 

AT  “ k21^VD~V^  + k22^YD~Y^  + k23^hD_h^  + k24^mD-m* 

* T T (6.2-20) 

D 

Here,  the  perturbed  controls,  pitch  angle,  0D~0 , and 
throttle,  Td-T,  are  related  to  the  state  deviations  through 
the  eight  guidance  gains  (k^,  The  subscript 

D on  the  trajectory  states  represents  the  values  of  the  state 
variables  along  the  desired  flight  path.  The  gain  values, 
kAj , represent  the  amount  by  which  the  controls  are  perturbed 
for  a given  value  of  state  deviation.  The  gain  k^ , for 
example,  represents  the  ratio  between  pitch  command  and 
velocity  deviation,  A0/AV,  while  the  gain  kgg  is  the  ratio 
of  throttle  perturbation  to  altitude,  AT/ Ah.  The  remaining 
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gains  represent  similar  control-to-state  perturbation  ratios 
as  defined  by  Eqs.  (6.2-10)  and  (6.2-20). 

Actual  gain  values,  computed  at  50-nm  Intervals 
along  a 525-nm  trajectory,  are  plotted  in  Fig.  32.  For  the 
purposes  of  simulating  perturbed  trajectories,  the  gain  values 
aro  linearly  interpolated  between  computed  points.  For  this 
reason,  the  plots  of  Fig.  32  present  the  gain  histories  in  the 
form  of  straight-line  segments. 

The  reference  trajectory  along  which  these  gain 
values  are  evaluated  consists  of  a climb  phase  (Oto  200  nm), 
a 200-nm  cruise  segment  (h  “ 39,000  ft),  and  a descent  seg- 
ment (400  to  525  nm).  The  six  major  gains  — A0/AV,  A6/Ay, 
A0/Ah,  AT/AV,  AT/Ay,  AT/ Ah  — are  seen  to  be  relatively  con- 
stant during  cruise,  with  major  gain  changes  occurring  dur- 
ing climb  and  descent  (Fig.  6.2-1).  The  two  gains  which 
relate  control  perturbations  to  deviations  of  vehicle  mass, 
A6/Am  and  AT/ Am,  obtained  average  values  of  only  0.017  deg/ 
1,000  lb  and  0.32  percent /I ,000  lb  along  this  reference  tra- 
jectory. Plots  of  these  gains  as  functions  of  range  are  not 
Included  in  the  figure. 

The  guidance  gain  plots  give  an  indication  of  the 
control  perturbations  which  are  employed  by  the  guidance 
algorithm  in  reducing  state  deviations.  For  example,  dur- 
ing the  cruise  segment,  the  plots  show  that  a 100-ft  altitude 
deviation  will  cause  a pitch-angle  perturbation  of  -1.5  deg 
and  a throttle  perturbation  of  approximately  -1.5  percent  to  be 
commanded.  The  true  merits  of  the  guidance  gain  values  can 
only  be  determined,  however,  by  observing  the  perturbed  tra- 
jectory response  provided  by  the  complete  set  of  gains. 
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Perturbation  state  and  control  profiles  for  each 
of  the  major  flight  segments  (climb,  cruise,  und  descent) 
are  presented  in  Figs.  33  through  35.  In  each  flight 
segment,  the  effects  of  initial  deviations  in  velocity 
(AV  ■ 17  fps),  flight  path  angle  (Ay  ■ 0.20  deg),  and  altitude 
(Ah  - 100  ft)  are  examined  individually.  The  initial  statt 
perturbations  are  set  to  the  maximum  valuos  which  define 
the  desired  guidance  system  response  (Table  15).  In  addi- 
tion to  plots  of  the  perturbed  state  profile,  the  control 
perturbations,  A0  and  AT,  commanded  by  the  linear-optimal 
guidance  algorithm  to  null  the  trajectory  deviations  also 
are  depicted.  An  initial  perturbation  in  a single  state 
variable  causes  the  other  state  variables  also  to  deviate 
somewhat  from  the  reference  trajectory.  These  secondary 
deviations  arise  through  the  coupling  of  the  trajectory 
states  defined  by  the  vehicle's  equations  of  motion.  Also, 
the  control  perturbations  employed  to  reduce  the  initial 
trajectory  deviation  simultaneously  affect  the  initially 
unperturbed  state  variable  profiles.  The  present  analysis 
demonstrates  that  the  induced  deviations  in  the  originally 
unperturbed  states  are  small  compared  to  the  specified 
trajectory  response  characteristics.  Therefore,  only  pro- 
files of  the  initially  perturbed  state  variables  ate  plotted. 
The  maximum  values  of  the  induced  deviations  are  noted  for 
each  perturbed  trajectory  simulated. 

Examination  of  the  perturbed  climb  profiles  (Fig. 

33)  shows  that  the  initial  trajectory  deviations  in  velocity 
and  flight  path  angle  are  eliminated  in  a range  of  less  than 
12  nm.  In  particular,  flight  path  angle  deviation  is  elimi- 
nated in  about  1 nm.  The  initial  altitude  perturbation  is 
substantially  reduced  in  the  figure  but  has  not  yet  reached 
a zero  perturbation  level  in  the  12-nm  range. 
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Figure  33  Climb  Phase  Perturbation  Response 
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Figure  34  Cruise  Phase  Perturbation  Response 
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Figure  35 


Descent  Phase  Perturbation  Response 
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The  perturbed  profiles  for  the  cruise  and  descent 
flight  segments,  Figs.  34  and  35,  respectively,  demonstrate 
characteristics  similar  to  those  of  the  climb  phase.  As 
before,  initial  perturbations  in  flight  path  angle  are  nulled 
almost  instantaneously.  The  initial  deviations  in  velocity 
and  altitude  take  longer  to  remove  but  are  significantly 
reduced  over  the  15-nm  profiles.  In  all  flight  modes,  not 
only  are  the  initial  perturbations  reduced  effectively,  but 
the  induced  secondary  perturbations  also  are  small. 

The  perturbed  flight  profiles  demonstrate  the 
excellent  trajectory-following  capabilities  of  the  linear- 
optimal  guidance  system.  Implementing  this  system  on  an 
actual  aircraft  would  necessitate  some  means  of  gain  deter- 
mination other  than  the  table  lookup  as  a function  of 
range  employed  in  the  current  simulations.  Because  each 
actual  flight  would  be  different,  i.e.,  path  length  and 
trajectory  profiles  would  vary,  an  adaptive  gain  determin- 
ation, based  upon  actual  flight  conditions,  would  be  needed. 
The  following  section  describes  an  adaptive  gain  scheduling 
technique. 


6.2.4  Correlation  of  Gains  With  Flight  Condition 

This  section  presents  a three-step  procedure 
for  scheduling  the  linear-optimal  guidance  gains.  The  gains 
are  scheduled  by  finding  functional  relationships  between 

t 

certain  aircraft  flight  conditions  and  the  control  gains  at 
these  flight  conditions.  This  results  in  an  adaptive  method 
of  gain  scheduling  which  could  be  utilized  in  an  on-board 
guidance  system. 
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Previous  methods  for  scheduling  control  gains  have 
been  successful  and  indicate  that  gain  scheduling  is  a sound 
approach  (Refs.  55  and  56).  The  methodologies  typically  are 
based  on  single-input/single  output  concepts  (e.g.,  main- 
taining constant  loop  gain);  therefore,  they  provide  inade- 
quate insight  for  scheduling  multivariable  systems. 


The  method  employed  here  (Ref.  57)  is  a logical 
extension  of  previous  work  to  multivariable  systems.  It 
involves  three  steps  and  places  minimum  reliance  on  past 
experience  and  intuition.  The  three  steps  are: 

• The  determination  of  a mean/standard 
deviation  table  which  aids  in  finding 
the  gains  that  should  be  kept  constant 
instead  of  being  scheduled. 

• The  determination  of  a correlation 
coefficient  table  which  aids  in 
specifying  the  strongly  correlated 
flight  conditions  that  will  be  used  in 
gain  scheduling  algorithms. 

• The  determination  of  relationships  (or 
curve  fits)  between  the  chosen  flight 
conditions  and  gains. 


The  first  step  in  the  procedure  is  to  determine  the 
mean  and  standard  deviations  of  the  optimal  guidance  gains 
over  a wide  range  of  flight  conditions.  The  mean  value  of  a 
gain,  k,  is  given  by 


(6.2-21) 


where  kA  is  the  value  of  the  gain  observed  at  flight  point  i , 
and  N is  the  number  of  points  at  which  the  gain  has  been 
computed.  The  gain's  standard  deviation,  o(k),  is  computed  as 
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o(k) 


(6.2-22) 


If  both  the  mean  and  standard  deviation  of  a 
gain  are  small,  then  the  gain  does  not  affect  the  control 
perturbations  and  need  not  be  included  in  the  closed-loop 
guidance  algorithm.  If  the  mean  value  is  significant,  but 
the  standard  deviation  is  small  with  respect  to  the  mean, 
then  the  gain  value  is  essentially  constant  over  the 
flight  conditions  examined.  Rather  than  being  scheduled 
as  a function  of  flight  conditions,  the  t;ain  should  b‘ 
considered  constant. 

The  means  and  standard  deviations  of  the  guidance 
gains  Inscribed  in  Section  6.2.3  are  presented  in  Table  16. 


TABLE  16 

MEAN  VALUE  AND  STANDARD  DEVIATION  OF 
PERTURBATION  GUIDANCE  GAINS 


Gain 

Mean  Value 

Standard  Deviation 
(Percent  of  Mean) 

1. 

A0/AV 

-7 . 04xl0“3  (deg/fps) 

108.8 

2. 

A0/Ay 

-8.61  (deg/deg) 

2.5 

3. 

A0/Ah 

-1.41x10-2  (deg/ft) 

11  . 1 

4. 

At) /Am 

1 . 72xl0“5  (deg/lb) 

19.9 

5. 

AT/AV 

-8  43x10“* (percent/ fps) 

10.7 

6. 

AT/Ay 

-5. 27xl0“1( percent /deg) 

35.8 

7. 

AT/Ah 

-1. 22xl0-2(percent/ft ) 

21.8 

8. 

AT  /Am 

3. 13x10-4 (percent /lb ) 

43.8 
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These  values  represent  gain  statistics  for  a twenty-point 
sample.  The  mean  gain  values  for  A0/Am  and  AT/ Am  are  ex- 
tremely small.  A 10,000-lb  weight  deviation  causes  a change 
in  pitch  command  of  only  0.17  deg  and  a change  in  throttle 
setting  of  only  3 percent.  These  gains  could  be  omitted, 
therefore,  when  implementing  an  on-board  system. 

Gains  2,  3,  and  5 ( A0/Ay , A0/Ah,  and  AT/AV  respec- 
tively), have  standard  deviations  equal  to  less  than  12  percent 
of  their  mean  value.  Because  the  variation  is  low,  it  is 
reasonable  to  schedule  these  gains  as  constants  in  the  guid- 
ance algorithm. 

The  remaining  gains  (1,  6,  and  7)  have  significant 
mean  values  and  do  vary  as  functions  of  flight  conditions. 

The  second  step  in  scheduling  these  gains  is  to  determine 
flight  conditions  with  which  they  are  correlated.*  The 
search  for  gain/flight  variable  dependencies  begins  by  deter- 
mining correlation  coefficients  between  the  gains  and  all 
available  flight  variables.  One  method  of  determining  the 
correlation  coefficient  between  a set  of  gains  (dependent 
variables)  and  flight  variables  (independent  variables)  is 
given  by  the  following  (Ref.  58): 


P(k,d) 


] cov [k , d] I 
o(k)  o(d) 


N 

J](ki-k)(di-d) 

i-1 


^i>i-s>2  J 


(6.2-23) 


♦Two  variables  are  correlated  if,  by  observing  the  value  of 
one  variable,  a conclusion  can  be  drawn  about  the  value  of 
the  other  variable. 
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In  Eq.  (6.2-23),  N is  the  number  of  flight  conditions  for 
which  the  gains  are  known,  k^  is  the  value  of  the  gain 
observed  at  flight  point  i,  and  d^  is  the  measured  value  of 
the  flight  variable  at  flight  point  i.  The  variable,  k, 
is  the  mean  value  of  the  gain,  and  d is  the  mean  value  of 
the  flight  variable.  The  closer  the  value  of  the  correla- 
tion coefficient,  p,  is  to  one,  the  better  the  correlation 
between  the  gain  and  the  flight  variable. 

The  independent  flight  variables  considered  for  gain 
correlation  in  the  current  study  are  true  airspeed  (TAS), 
atmospheric  pressure  ratio  (6),  weight  divided  by  pressure 
ratio  (W/6),  weight  divided  by  dynamic  pressure  (W/q),  and 
angle  of  attack  (a).  While  by  no  means  an  exhaustive  set, 
these  flight  variables  represent  major  dynamic  variations 
with  vehicle  flight  conditions. 

In  addition  to  determining  the  correlation  coeffi- 
cients for  the  flight  variables  stated,  correlations  between 
the  guidance  gains  and  functions  of  the  flight  variables  such 
as  the  square,  the  inverse,  the  inverse  square,  etc..,  of  each 
flight  variable  can  be  obtained.  Also,  correlation  coeffi- 
cients relating  the  gains  and  a combination  of  the  flight 
variables  can  be  determined  assuming  the  following  functional 
dependency : 

k = bQ  + b1d1  + b2d2  + . . . +bndn  (6.2-24) 

A 

where  k is  the  computed  guidance  gain,  d^  to  dfi  are  the  flight 
variables,  and  the  coefficients,  b^  provide  the  "best"  gain 
estimate  as  a function  of  the  varying  flight  variables.  The 
estimate  provided  by  the  bA  is  "best,"  in  that  it  provides  a 
least-squares  fit  to  the  data  points.  That  is,  the  k^  are 
computed  to  minimize  the  following  cost  function: 
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(6.2-25) 


j ■ t,  <w2 

i=l 

The  techniques  for  determining  the  functional  coefficients, 
b^  as  well  as  the  correlation  coefficients  for  the  multiple 
flight  condition  relationships  are  detailed  in  Ref.  57. 

A summary  of  the  correlation  coefficients  for  the 
three  guidance  gains  (A0/AV,  AT/ Ay,  AT/Ah)  of  interest  here 
is  presented  in  Table  17.  Correlation  coefficients  are 
shown  for  those  flight  variables  which  are  highly  corre- 
lated with  these  gains.  The  first  gain  A0/AV,  correlates 
strongly  with  both  W/q  and  a.  Because  a may  be  a difficult 
quantity  to  measure  in  flight,  it  would  not  bo  unreasonable 
to  schedule  A^/AV  as  a function  of  W/q  only.  The  sixth  guic  - 
ance  gain,  AT/Ay  also  correlates  most  strongly  with  W/q. 

The  combination  of  W/q  and  a does  not,  in  fact,  yield  a 
larger  correlation  for  this  gain. 

TABLE  17 

GUIDANCE  GAIN  CORRELATION  COEFFICIENTS 


Gain 

Flight 

Conditions 

W/q 

a 

(1/TAS)2 

W/q, a 

1.  A0/AV 

0.83 

0.89 

0.81 

0.89 

2.  AT /Ay 

0.93 

0.90 

0.83 

0.93 

3.  AT/Ah 

0.26 

0.51 

0.45 

0.80 

The  gain  relating  throttle  commands  to  altitude 
deviations,  AT/Ah,  has  its  strongest  single-flight-variable 
correlations  with  a and  with  (1/TAS)2.  The  correlation  of 
AT/Ah  with  a linear  combination  of  these  two  variables  did 
not  show  any  improvement;  however,  the  correlation  of  AT/Ah 
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with  the  combined  flight  variables,  W/q  and  a,  correlated 
significantly,  with  a coefficient  equal  to  0.80. 

Employing  W/q  to  schedule  both  A0/AV  and  AT/Ah 
and  the  combination  of  W/q  and  a to  schedule  AT/Ah  results 
in  the  following  gain  scheduling  equations: 

A6/A  * 2.994X10"2  - 2.344xl0"5  W/q  (6.2-26) 

AT/Ay  = 0.5245  - 6.707xl0"4  W/q  (6.2-27) 

AT/Ah  = 0. 5257xl0”2  - 0.1755xl0-4  W/q  + 0.4386X10-2  a 

(6.2-28) 

In  the  above  equations,  the  coefficients  of  the  flight  vari- 
ables are  those  obtained  from  the  polynominal  regression 
analysis  described  in  Ref.  57.  The  gain  A0/AV  has  units  of 
deg/fps,  while  the  gains  AT/Ay  and  AT/Ah  have  units  of  per- 
cent/deg and  percent/ft,  respectively  . The  flight  variable 
W/q  is  expressed  in  units  of  ft^,  and  a is  measured  in  degrees. 
The  use  of  these  coefficients  provides  for  the  correlations 
between  the  guidance  gains  and  flight  variables  presentel 
in  Table  6.2-2. 

This  section  has  described  a means  by  which  the 
linear-optimal  guidance  gains  can  be  scheduled  in  an  adap- 
tive fashion.  This  enables  the  use  of  these  gains  in  an 
on-board  guidance  algorithm.  The  results  presented  provide 
reasonably  good  correlations  between  the  gains  and  flight 
conditions.  Better  correlations  could  be  obtained  by  em- 
ploying more  independent  flight  variables  and  using  higher- 
order  polynomial  representations.  The  purposes  of  this  sec- 
tion are  served,  however,  by  the  demonstration  of  the  appli- 
cation of  the  correlation-coefficient  gain-scheduling  method 
to  the  fuel-optimal  perturbation  guidance  problem. 
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6.3 


STRUCfURE  OF  THE  THROTTLE /ENERGY  MANAGEMENT  ALGORITHM 


The  previous  section  has  shown  how  the  fuel-mini- 
mizing guidance  algorithm  can  be  separated  into  two  parts: 
the  generation  of  a desired  flight  path,  Xp,  and  its 
associated  nominal  guidance  variable,  Upi  and  the  correction 
of  flight  path  errors  using  a linear  guidance  law 

Au  - -KAx  « -K(Xp-x)  (6.3-1) 

The  total  guidance  command,  u^,  is  the  sum  of  the  nominal 
value  and  the  linear  perturbation: 

U£  - uD  + Au  (6. 3-2 ) 

It  can  be  seen  in  Eq.  (6.3-2)  that  the  guidance  command  has 
an  open- loop  component,  Up,  which  corresponds  to  the  trim 
values  of  pitch  angle,  0,  and  throttle  setting,  T,  plus 
the  closed-loop  component,  Au,  derived  from  the  feedback 
of  the  difference  between  the  desired  and  actual  state  to 
the  pitch  and  throttle  channels.  The  result  is  a combined 
autopilot/autothrottle  command  which  provides  coordinated 
control  of  altitude  and  velocity. 

The  guidance  commands,  0C  and  Tc,  should  be  con- 
sidered as  set  points  for  the  inner-loop  control  law.  Thus, 
it  is  assumed  that  the  throttle/energy  management  logic  is 
followed  by  pitch  attitude  and  throttle  setting  regulators, 
as  shown  in  Fig.  36.  Typically,  the  pitch  channel  would 
use  attitude  and  rate  feedback,  while  the  throttle  channel 
could  use  Engine  Pressure  Ratio,  EPR,  and  low-speed  com- 
pressor RPM,  N^. 

Figure  36  illustrates  the  basic  features  of  the 
throttle/energy  management  structure.  The  desired  flight 
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path  Is  specified  by  the  flight  path  generator  logic,  which 
is  discussed  below.  The  linear  perturbation  guidance  law 
and  gain  scheduling  were  discussed  in  the  previous  section, 
and  the  trim  logic  is  presented  in  this  section.  Require- 
ments for  crew  display  appear  in  a later  section.  The 
advantages  of  this  guidance  structure  are  that 

• The  system  generates  near-optimal 
flight  paths; 

• These  flight  paths  are  followed  using 
reasonable  commands;  and 

• The  throttle/energy  management  is  directly 
applicable  to  all  large  transport  and 
bomber  aircraft , although  the  present 
results  are  numerically  specific  to  the 
C-141A. 


It  is  important  that  this  guidance  structure  be 
implemented  in  a practical  way,  if  it  is  to  be  used  on  an 
operational  basis.  This  means  that  fuel-optimal  flight 
paths  must  be  stored  in  an  efficient  manner  while  providing 
sufficient  adaptability  to  changing  flight  conditions  and 
flexibility  with  respect  to  mission  constraints.  Unless 
reasonable  simplifications  are  made  in  the  storage  and 
computation  of  the  Xp  and  Up  profiles,  the  computer  burden 
of  the  throttle/energy  management  algorithm  could  be 
excessive.  This  section  illustrates  that  partitioning  each 
mission  into  flight  path  segments  leads  to  efficiency  in 
specifying  Xp,  while  using  trim  integrators  to  compute 
Up  is  consistent  with  managing  fuel  use  with  a simple  sys- 
tem. 


6.3.1  Trim  Computation 

The  ad  hoc  guidance  logic  used  t:o  generate  param- 
etric-optimum and  reference  flight  paths  (Chapter  3)  pro- 
vides the  design  algorithm  for  the  trim  logic  (for  those 
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flight  conditions  in  which  the  trim  settings  are  not  defined 
explicitly).  Integrating  these  rate-command  guidance  laws 
(Eqs. (3.4-1)  and  (3.4-2))  leads  to 


e 


c 


(Vjj-V)dt  + kQ  (VD-V) 
2 


(6.3-3) 


(6.3-4) 


The  linear-optimal  guidance  laws  derived  in  the  previous 
section  (Eqs. (6. 2-19)  and  (6.2-20))  provide  disturbance  rejec- 
tion by  using  proportional  feedback.  Dropping  the  propor- 
tional term  in  Eq.  (6.3-3),  the  ad  hoc  guidance  laws  provide 
integral  compensation  which  assures  zero  steady-state  error 
in  following  the  aircraft's  velocity  and  specific  total 
energy.  The  trim  values,  0D  and  T^,  corresponding  to  the 
desired  flight  path  are  computed  by 


0, 


/<v 


V)dt 


(6.3-5) 


\/(v 


E)dt 


(6.3-6) 


where  k@D  and  kTD  are  chosen  to  provide  trim  adjustment  times 
which  do  not  conflict  with  the  proportional  perturbation 
guidance.  (For  example,  the  gain  kgp  is  chosen  to  be  a small 
percentage  of  the  value  of  A9/AV  given  in  Eq.  (6.2-26).  Since 
both  velocity  and  altitude  contribute  to  E (Eq.  (3.4-3)),  k-rD 
is  small  compared  to  both  AT/AV  and  AT/Ah  in  Table  16  and 
Eq.  (6.2-28). 


The  trim  logic  should  be  relatively  insensitive  to 
the  oscillatory  phugoid  mode  for  several  reasons.  The  gains 
in  Eqs. (6.3-5)  and  (6.3-6)  are  small,  and  periodicities  in  the 
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integrands  will  be  averaged  out.  Furthermore,  the  phugoid 
mode  is  a constant-energy  oscillation,  and  it  should  be 
largely  unobservable  in  the  integrand  of  the  throttle  equa- 
tion. Finally,  the  perturbation  guidance  law  augments  damping 
of  the  phugoid  mode 

6.3.2  Fuel-Minimizing  Flight  Path  Segments 

Integrated  flight  profiles  must  be  computed  to 
determine  the  fuel-optimal  means  of  flying  from  one  point  to 
another,  but  once  these  profiles  are  understood,  distinct 
flight  path  segments  can  be  identified  and  used  as  the  basis 
for  real-time  implementation.  The  segments  are  useful  because 
not  every  parameter  is  significant  to  fuel  use  in  every  seg- 
ment; hence,  efficiency  is  gained  by  using  the  minimum  number 
of  appropriate  parameters  to  define  each  segment. 

It  has  been  shown  that  there  are  two  basic  flight 
path  segments  for  short-range  trajectories  (climb  and  descent) 
and  three  segments  for  longer  trajectories  (climb,  cruise, 
and  descent).  Figure  37  illustrates  that  8 segment  types 
can  be  defined  for  a throttle/energy  management  system. 

These  are 

• Takeoff-Climb  Transition 

• Climb 

• Climb-Cruise  Transition 

• Cruise 

• Step  Climb 

• Cruise-Descent  Transition 

• Descent 

• Descent-Approach  Transition 

The  following  paragraphs  describe  the  major  features 
of  each  segment.  In  every  case,  V represents  true  airspeed 
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Figure  37  Flight  Path  Segments  for  the  Throttle/ 
Energy  Management  System 

(TAS),  and  h is  the  pressure  altitude.  It  is  assumed  that 
estimates  (or  measurements)  of  weight,  range  to  go,  wind 
speed,  and  sea  level  pressure  are  available  as  each  segment 
progresses.  In  many  cases,  indicated  airspeed  (IAS)  or 
Mach  number  is  a more  useful  variable  than  TAS;  hence,  V 
typically  is  a function  of  h.  Numerical  definitions  which 
support  these  descriptions  are  given  in  the  next  section. 

Takeof f-Climb  Transition  - This  segment  is  initiated 
when  the  landing  gear  and  flaps  are  fully  retracted  and  when 
the  throttle/energy  management  system  is  enabled  by  the 
crew.  It  lasts  for  approximately  one  minute.  Desired  values 
of  the  state  and  control  variables  initially  are  as  follows: 

• Vjj  * V,  i.e.,  the  current  value 
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• hjj  ■ Zero  to  current  value 

• hp  ■ Constant  to  current  value 

• 0n  ■ Value  which  provides  horizontal 

flight  (trim  off) 

• TD  ■ 1.0  (trim  off) 

Velocity-error  feedback  and  perturbation  throttle  commands  are 
suppressed;  hence,  the  throttle  command  is  entirely  "open- 
loop",  and,  perturbation  pitch  angle  regulates  altitude  as 
velocity  is  increasing.  During  the  initial  phase  of  the 
transition,  0p  decreases  as  velocity  increases. 


During  the  climb  phase,  0 is  used  as  a velocity  con- 
trol,! and  h is  treated  as  an  independent  variable;  therefore, 
a gradual  change  in  the  function  of  0 must  be  accomplished 
as  the  climb  speed  is  approached.  The  change  process  is 
called  "fade  out/fade  in",  and  it  is  obtained  by  interpolation 
of  VD,  hp,  and  hp.  Assume  that  the  fade-out/fade-in  time  is 
tp  sec  and  that  the  vehicle  is  accelerating  at  a rate  of  V. 

The  change  begins  at  time,  t , when  the  measured  velocity  is 

c 


v(tc) 


UCLIMB  c c 


(6.3-7) 


At  the  beginning  of  the  change,  Vp  is  V(tc),  hp  is  zero, 
and  hp  is  the  altitude  at  the  beginning  of  the  transition. 
During  the  transition,  0p  is  fixed  at  @D(tc)>  and 

VD(t)  - V(tc)  ♦ V(tc)  At/tF  (6.3-8) 


hp(t)  - h(t ) At/tj,  (6.3-9) 

hD(t ) - W ♦ (h(t)-hD(tc)]At/tF  (6.3-10) 


where 


At  - t 


(6.3-11) 
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In  other  words,  the  desired  velocity  gradually  takes  its 
climb  value,  and  ftp  and  hp  gradually  become  the  current 
values  (nullifying  the  errors  in  the  perturbation  guidance 
equations).  The  transition  is  complete  when  At  » ty  . 

Climb  - This  segment  begins  when  the  IAS  has 
reached  its  iuitial  climb  value,  and  it  lasts  for  15  to  40 
min  (depending  on  vehicle  weight),  the  initial  cruising 
altitude  (for  long-range  flight  paths),  and  the  transition  to 
descent  (for  short-range  flight  paths).  Desired  values  of 
the  state  and  control  variables  are  the  following: 

• Vp  = Value  which  yields  constant  IAS 

or  M,  as  defined  below 

• ho  “ Current  value  of  h 

• hD  * Independent  variable 

• 0D  “ k0Dj(VD-V)dt  (trlm  on) 

• Tn  ■ 1.0,  or  piloted-selected  reduced 

L setting  (trim  off) 

VD  is  a function  of  vehicle  weight,  throttle  setting, 
altitude,  air  temperature,  mission  range,  and  air  traffic  con- 
trol constraint.  At  most,  there  are  three  sub-segments  of 
the  velocity  profile: 

• Vn  provides  250  KIAS  below  10,000-ft 

altitude. 

• Vn  provides  optimal  IAS  until 

M - 0.721. 

• Vjj  provides  M ■ 0.721  until  cruise 

condition  is  reached. 

The  first  segment  is  absent  if  ATC  approval  for  higher 
speed  is  obtained,  or  if  the  optimal  IAS  is  less  than  250  KIAS 
(as  is  the  case  for  light  vehicles).  The  third  segment  is 
absent  on  short-range  missions  which  do  not  reach  cruising 
conditions,  as  shown  in  Chapter  4,  the  optimal  IAS  during 
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climb  increases  as  vehicle  weight  increases,  and  the  climb 
IAS  is  increased  for  short-range  missions  without  cruise  seg- 
ments; the  desired  value  is  computed  taking  these  factors  into 
account.  It  also  can  be  shown,  using  energy-state  methods, 
that  climb  IAS  is  approximately  proportional  to  the  square 
root  of  engine  thrust;  therefore,  a 10-percent  thrust  reduc- 
tion should  be  accompanied  by  a 5-percent  speed  reduction. 
Climb  thrust  reduction  increases  fuel  use  but  may  be  selected 
for  reasons  of  reduced  engine  wear  or  in  an  engine-out 
condition.  In  addition,  the  crew  should  have  the  option  of 
specifying  desired  values  of  IAS  and  II  for  comparative  test- 
ing. 


Initialization  logic  in  the  throttle/energy 
management  system  must  determine  whether  or  not  each  mission 
contains  a cruise  segment.  If  cruise  is  absent,  the  descent- 
initiation  prediction  procedure  described  in  Section  6.3.3  must 
be  executed  during  the  climb  segment.  If  cruise  is  present, 
the  decision  to  switch  to  climb-cruise  transition  is  made  in 
the  climb  segment.  The  switch  is  made  when  the  time  to 
climb  to  cruise  altitude  (either  computed  from  W/6  or  speci- 
fied by  the  crew)  equals  the  desired  transition  time,  assuming 
that  the  climb  proceeds  at  hp.  If  cruise  is  absent,  the 
climb-cruise  transition  is  immediately  followed  by  cruise- 
descent  transition,  and  it  is  selected  on  the  basis  of  range 
to  go. 


Altitude  and  altitude-rate  errors  are  suppressed 
in  the  linear  A8  and  AT  controllers  because  altitude  is  used 
as  the  independent  variable  for  climb;  therefore,  the  present 
altitude  is  the  desired  value,  and  the  error  is  zero.  Throttle 
commands  are  limited  to  1.0  or  less.  When  climbing  at  nor- 
mal-rated thrust  (Tp  ■ 1.0),  negative  AT  commands  are  executed, 
but  positive  AT  commands  are  not.  For  Tp  < 1.0,  perturbation 
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commands  of  both  signs  are  issued  up  to  the  point  that 
Tc  * 1.0. 


Climb-Cruise  Transition  - This  segment  begins  when 
the  climb-segment  prediction  indicates  that  cruise  is  imminent. 
The  dominant  features  of  this  segment  are  the  adjustment  of 
velocity  to  cruise  value,  the  reduction  of  climb  rate,  and 
the  achievement  of  the  desired  initial  altitude  for  cruising 
flight.  The  principal  functions  of  this  segment  are  to 
"fade  out"  climb  desired  values  and  to  "fade  in"  cruise 
desired  values  using  transition  interpolation;  therefore, 
the  cruise  option  selected  by  the  crew  (described  below) 
must  be  accessible  during  this  segment  and  should  be  used  to 
compute  initial  values  of  the  desired  cruising  flight  vari- 
ables. Given  this  information  plus  the  range  to  the  final 
destination,  the  desirability  of  step  climbs  during  cruise 
should  be  evaluated.  The  transition  segment  is  complete 
when  velocity  and  altitude  reach  their  respective  cruise 
values. 


Assuming  a transition  time  of  ty^,  the  desired 
velocity  changes  from  the  climb  value,  at  time,  tk, 

to  the  cruise  value,  VD£R: 


CL 


♦ (Vr 


-V 


CR 


n )At/tp 
°CL  Fk 


(6.3-12) 


where 


At  * t - tjj 


(6.3-13) 


while  the  altitude-rate  and  altitude  errors  are  established 
by  fading  in  the  desired  cruise  values,  h^R  and  hp^,R: 


hD(t)  - h(t)  + (hjj  -h(t)JAt/tF 

CR  k 


(6.3-14) 
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(6.3-15) 


hD(t)  - h(t)  ♦ (hp  -h(t )] At/tp 

CR  k 

The  transition  should  begin  early  enough  to  prevent 
an  altitude  overshoot  with  the  established  climb  rate  and 
closed-loop  guidance  response.  Section  6.2  indicates  that 
Ay  response  (which  is  equivalent  to  Ah  response)  of  the 
perturbation  controller  is  rapid;  therefore,  t^  can  be  chosen 
on  the  basis  of  load  factors  required  to  reduce  h from  the 
climb  value  to  the  cruise  value. 

On  a short-range  flight,  the  climb-cruise  transition 
is  followed  immediately  by  the  cruise-descent  transition.  In 
this  case,  the  primary  function  of  the  climb-cruise  transition 
is  to  arrest  the  climb  rate  and  to  begin  throttle  reduction. 

A nominal  reduction  in  desired  velocity  should  be  specified 
for  vDcr*  along  with  zero  value  for  Hdcr  and  current  altitude 
for  h0cR. 


Cruise  - Three  options  are  implemented  in  the  cruise 
segment  logic: 

• Breguet  Climbing  Cruise 

• Constant-Altitude  Cruise  With  Optimal 

(Varying)  Velocity 

• Constant-Altitude  Cruise  with  Constant 

Velocity 

In  the  first  option,  h^  increases  to  maintain  W/6  at  its 
optimal  value,  and  V^  is  chosen  to  maintain  optimal  M. 

Optimal  W/6  is  a function  of  c.g.  location,  drag  coefficient, 
and  air  temperature,  while  optimal  M is  a function  of  air 
temperature  and  wind  speed.  Mach  number  decreases  in  con- 
stant-altitude cruise  for  the  second  option,  providing  a 
relative-optimum  condition  as  fuel  is  burned  and  W/6  varies. 
The  third  option  allows  a conventional  cruise  segment  to  be 
flown  automatically  for  comparative  purposes. 
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Subject  to  these  options,  the  desired  values  of 
the  state  and  control  variables  are  as  follows: 


Vp  * Value  which  yields  desired  IAS 
or  M,  as  defined  above 

hD  = °,  or  function  of  Vjj  and  h[>  in 
climbing  cruise  (Section  6.3.3) 

hp  = Value  which  provides  optimal  W/6, 
or  constant,  as  defined  above 


(Vp-V)dt  (Trim  on) 


(Trim  on ) 


The  full  perturbation  guidance  algorithms 
(Eq.  (6.2-19)  and  (6.2-20))  are  implemented  for  altitude 
and  velocity  control  during  the  cruise  phase. 

A particular  cruise  segment  can  end  in  one  of  two 
ways:  a step  climb  (followed  by  another  cruise  segment)  or 

a descent.  Background  computations  to  predict  when  either 
segment  change  could  occur  should  be  conducted  during  the 
cruise  segment.  The  step  climb  prediction  is  computed  only 
during  constant-altitude  cruise,  and  it  is  assumed  that  the 
climb  must  occur  in  a 4,000-ft  increment.  The  current  air- 
craft weight,  c.g.  location,  air  temperature  (adjusted  to 
account  for  standard  lapse  rate),  and  any  known  variation 
in  the  particular  aircraft's  drag  coefficient,  are  used  to 
compute  the  optimal  W/6  and  the  corresponding  pressure 
altitude.  A step  climb  is  indicated  when  the  optimal  alti- 
tude is  2,000  ft  above  the  current  hp.  Prediction  of  the 
descent  initiation  is  carried  out  by  estimating  the  range 
increment  during  descent  plus  the  range  increment  during 
cruise-descent  transition  and  subtracting  the  result  from 
the  range  to  go  (Section  6.3.3). 
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Although  each  prediction  is  individually  useful,  the 
two  can  be  used  together  to  decide  whether  or  not  a step 
climb  is  desirable.  Following  the  4,000-ft  climb,  the  air- 
craft will  be  2,000  ft  above  its  optimal  altitude.  The  step 
climb  produces  measurable  fuel  savings  if  the  aircraft  spends 
enough  time  at  the  higher  altitude  for  the  fuel  burnoff  to 
make  the  new  altitude  optimal.  Otherwise,  the  fuel  consumed 
in  gaining  potential  energy  is  not  offset  by  the  improved 
specific  range.  If  the  range  to  be  flown  to  descent  initia- 
tion (computed  for  the  higher  altitude)  would  not  require 
sufficient  fuel  burnoff,  then  the  step  climb  should  not  be 
made. 


Step  Climb  - The  procedures  followed  during  the  step 
climb  are  similar  to  those  of  the  earlier  climb  segment. 
Altitude  again  is  the  independent  variable,  full  throttle  is 
used,  and  pitch  trim  is  engaged.  Referring  to  Fig.  10,  it 
can  be  seen  that  the  optimal  Mach  number  on  the  higher  con- 
stant-altitude cruise  segment  will  be  greater  than  it  was  on 
the  lower  cruise  segment  (since  the  relative-optimal  M is 
defined  by  the  zero-slope  tangeucies  on  the  constant-specific 
range  contours,  Fig.  10);  therefore,  VD  should  increase 
linearly  during  the  step  climb  to  match  the  relative-optimal 
cruise  points  at  the  beginning  and  end  of  the  climb. 

Cruise-Descent  Transition  - This  important  flight 
Dath  segment  is  a near-constant-altitude,  minimum-thrott le 
deceleration  from  the  cruise  Mach  number  to  the  IAS  required 
for  (L/D)  descent.  The  primary  objective  of  this  seg- 

max 

ment  is  to  provide  a smooth  transition  which  does  not  excite 
the  phugoid  mode.  (The  latter  would  lead  to  velocity- 
altitude  oscillations,  which  would  be  perceived  as  load- 
factor  oscillations  by  passengers  and  the  crew. ) Unlike  the 
takeoff-climb  and  climb-cruise  transitions,  the  distance 
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travelled  during  the  transition  is  important  for  flight 
path  management . The  descent  initiation  point,  i.e.,  the 
start  of  the  cruise-descent  transition,  must  be  predicted 
for  maximum  fuel  savings,  so  the  transition  range  increment 
should  be  designed  to  be  easy  to  compute.  This  infers  that 
the  transition  velocity  profile,  which  is  integrated  to 
obtain  wind-relative  distance  travelled,  should  be  well- 
defined  and  followed  closely. 

Desired  state  and  control  variables  for  the  cruise- 
descent  transition  are  the  following: 

• vd  “ Cruise  VD  to  descent  VD 

• hp  ■ Zero  to  descent  value 

• hD  * Cruise  altitude  to  descent  value 

• eD  “ keD  $(VD-V)dt  (Trim  on) 

• TD  - J(ED-E)dt  (Trim  on) 


This  segment  can  be  sub-divided  into  the  throttle- 
down  and  deceleration  phases,  both  of  which  involve  fade-out/ 
fade-in  logic.  In  the  first  phase,  the  total  throttle 
command  is  reduced  linearly  from  its  cruise  setting,  T^,  to 
the  descent  setting  (normally  idle),  T^,  over  a period  of 
ty  sec.  Assuming  the  transition  begins  at  t . , 

L ^ 

T (t)  - T (t  ) - [T  (t  )-TM] At  /t„  (6.3-16) 

C Ct  Ct“tIt 

where 

AtT  - t - tt  (6.3-17) 

TD  is  ignored  during  the  throttle  reduction  and  reinitialized 
to  zero  at  the  end  of  the  phase.  During  the  throttle-down 
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there  is  a nominal  decrease  in  Vp  to  the  value  , hp  is 
zero,  and  hp  is  the  cruise  altitude.  During  the  deceleration 
the  desired  velocity  decreases  from  to  and  hp  and 

hD  are  faded  to  the  initial  descent  values  over  a period  of 

tp.  sec  using  Eq.  (6.3-9)  and  (6.3-10). 

1 d 

Descent  - Guidance  during  the  descent  is  analogous 
to  guidance  during  .he  climb,  but  there  is  an  additional 
function  to  be  conducted:  range  control.  If  the  descent 

has  been  initiated  at  the  proper  point  , then  the  vehicle 
will  arrive  at  the  terminal  altitude  and  range  with  no 
error;  however,  it  is  more  likely  that  disturbances  would 
cause  the  vehicle  to  miss  the  terminal  condition  in  an  open- 
loop  descent.  It  is  necessary,  therefore,  to  minimize  posi- 
tion errors  as  well  as  fuel  use  during  this  trajectory  seg- 
ment. By  reducing  flight  path  errors  early  in  the  descent, 
range  control  allows  operational  constraints  to  be  satisfied 
with  a minimal  fuel  penalty. 


The  basic  descent  profile  is  a maximum  lift-to-drag 
ratio  flight  path  with  throttle  set  at  idle.  (Reference  47 
indicates  that  air  conditioning  and  pressurization  can  be 
maintained  at  idle  throttle,  so  there  is  no  need  to  consider 
powered  descent.)  This  is  a maximum-range  glide  path;  in 


the  absence  of  wind,  the  corresponding  constant  flight  path 
angle,  , is 


opt  ’ 


Y «.  * -cot-*(L/D) 
’opt  v ' 'max 


(6.3-18) 


which  is  independent  of  aircraft  weight.  Consequently,  there 
is  a straight-line  relationship  between  altitude  and  range 
on  the  optimal  descent  path,  and  it,  too,  is  independent  of 
weight.  Calling  the  range  to  go  (to  reach  the  destination), 
Rg0,  and  the  terminal  altitude,  h Gp,  the  desired  altitude  is 
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‘b  “ hG0  " RG0  tan  Yopt 


(6.3-19) 


The  associated  altitude  rate  is  a function  of  weight, 
as  is  the  time  to  go,  t^,  because  the  dynamic  pressure  (and, 
therefore,  the  indicated  airspeed)  required  to  maintain  "1-g" 
flight  varies.  The  desired  altitude  rate  is 

(6.3-20) 


(6.3-21) 

and  CLopt  is  the  lift  coefficient  at  (L/D)max-  Note  that 
the  specific  total  energy  during  descent  is  fully  defined  by 
Eq.  (6.3-19)  and  (6.3-21). 

The  descent  flight  path  must  be  adjusted  for  wind 
effects.  The  air-relative  variables  are  unchanged,  but  the 
effective  flight  path  angle,  ye,  is  made  shallower  by  tail 
winds  (vw>0)  and  steeper  by  head  winds  (Vw<0).  Since 
tan  y£  = -h /ft, 

Ye  * -tan_1[hD/(VDcosYopt+Vw)J  (6.3-22) 

and,  assuming  a constant  wind  profile, 

»D  * hGO  ' RG0tanYE  (6.3-23) 

The  desired  state  and  control  variables  are  then  defined  by 

• VD  » Eq.  (6.3-21) 

• hD  » Eq.  (6.3-20) 
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Eq.  (6.3-23) 


• "d  * 


* TD  ' kT-J(ED  -E)dt  (Trim  on) 

There  are  two  choices  for  descending  through  a 
varying  wind  profile:  either  update  Yjj  as  the  variations 

occur  or  use  an  average  value  of  to  compute  a fixed 
descent  path.  The  first  approach  is  more  nearly  fuel- 
optimal,  but  the  second  produces  a more  predictable  path 
and  is,  therefore,  more  compatible  with  air  traffic  con- 
trol. 


j(VD-V)dt  (Trim  on) 


The  full  perturbation  guidance  algorithms  are 
employed  during  descent.  For  the  most  part,  it  can  be 
assumed  that  pitch  control  will  be  used  when  the  vehicle 
is  flying  "long"  (i.e. , above  or  faster  than  the  nominal 
flight  path)  and  that  throttle  control  will  be  used  when 
the  vehicle  is  flying  "short"  (the  opposite).  In  the 
first  case,  E > Ep,  which  normally  would  decrease  TD:  how- 
ever, Tp  is  nominally  zero  and  cannot  become  smaller, 

so  6n  will  increase  to  lower  the  specific  kinetic  energy 
2 u 

(V  /2g).  In  the  second  case,  E < E^,  and  the  vehicle  does 
not  possess  sufficient  energy  to  glide  to  its  destination. 
Tq  will  increase  to  provide  the  necessary  energy. 

Descent-Approach  Transition  - Once  the  vehicle  has 
descended  from  altitude,  the  sink  rate  must  be  arrested  to 
restore  the  aircraft  to  level  flight . The  descent-approach 
transition  can  be  regarded  as  establishing  a low-altitude 
cruise  condition  by  increasing  throttle  setting  and  pitch 
angle.  The  desired  state  and  control  variables  at  the  end 
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of  the  transition  are: 


• vd  “ Pilot-selected  value 

• " 0 

• hp  * Pilot-selected  value 

• 0D  - kQ  J(VD-V)dt  (Trim  on) 

• Td  « ^ (Ejj-E)dt  (Trim  on) 

The  fade-out/fade-in  logic  used  for  the  climb- 
cruise  transition  can  be  used  here  as  well.  The  guidance 
algorithm  enters  a cruise  segment  at  the  completion  of  the 
transition. 

6.3.3  Computational  Details 

Additional  aspects  of  the  throttle/energy  manage- 
ment algorithm  are  introduced  here.  These  include  generation 
of  the  nominal  flight  paths,  prediction  of  the  descent 
initiation  point,  computation  of  descent  speed,  and  support 
calculations  for  crew  display  and  throttle  control  inter- 
face. 


Nominal  Profiles  - In  order  to  provide  efficient 
storage  of  data,  nominal  flight  profiles  are  stored  in 
natural  dimensions,  and  auxiliary  calculations  are  imple- 
mented where  necessary.  Rather  than  storing  true  airspeed, 
which  is  only  indirectly  related  to  energy  management, 
nominal  values  of  indicated  airspeed  and  Mach  number  are 
stored,  as  appropriate.  Similarly,  W/6  is  stored  in  place  of 
altitude.  Nonzero  values  of  altitude  rate  are  used  only  in 
the  descent  and  Breguet-cruise  guidance  algorithms.  If 
future  testing  indicates  the  need  for  improved  damping  during 
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climb,  the  hp  required  to  form  the  feedback  error  can  be 
computed  from  aircraft  characteristics,  Vp,  and  hp.  During 
a constant-IAS  segment, 


F-CDSPvg/2 


(6.3-24) 


where  F is  engine  thrust,  and  air  density  (p)  and  scale  height 
(6)  are  functions  of  hp.  During  a constant-M  segment,  only 
potential  energy  changes  (neglecting  the  effects  of  sound- 
speed  variation  on  U),  leading  to 

hD*  Y ( F-CpSpVp/2)  (6.3-25) 


If  the  aircraft  is  climbing  at  constant  M to  keep  W/6  constP.it, 


W 

(W/6)  || 


(6.3-26) 


Chapter  4 has  indicated  that  IAS,  M,  and  W/6  are 
smoothly  varying  functions  of  flight  conditions  and  vehicle 
characteristics ; therefore,  these  quantities  are  most  easily 
stored  in  functional  form.  The  algorithms  used  to  develop 
gain  schedules  in  Section  6.2  can  be  employed  to  schedule 
desired  values  of  IAS,  M,  and  W/6  as  functions  of  the  major 
variables. 

Descent  Initiation  Prediction  - The  point  at  which 
the  cruise  segment  ends  and  the  descent  begins  must  be  com- 
puted precisely  if  maximum  fuel  savings  are  to  be  realized. 
Given  the  cruise  conditions  during  medium-  to  long-range 
flight  and  the  climb  conditions  during  short-range  flight, 


143 


the  descent  should  be  initiated  when  Rqq  equals  the  sum  of 
four  incremental  ranges: 

• Climb-cruise  transition  range  (short- 
range  flight  path  only) 

• Throttle-down  range 

• Deceleration  range 

• Idle  thrust  descent  range 

On  a short-range  mission,  the  descent  is  initiated 

during  the  climb,  and  the  climb-cruise  transition  segment 
is  used  to  arrest  the  climb  rate.  The  range  travelled  dur- 
ing this  segment  is  approximately 

Ro  - IV,  *<v  +V  >/21t  (6.3-27) 

Q LL  Ln  K 

where  is  the  prevailing  wind  velocity. 
w0 

The  throttle-down  velocity  and  range  can  be  computed 
using  an  analytical  expression  for  the  deceleration  which 
occurs  during  a linear  decrease  in  throttle  setting. 

Assuming  CD  and  altitude  remain  constant  during  the  throttle 
reduction , 

V - -ktAtTV2  (6.3-28) 

where  AtT  is  defined  by  Eq.  (6.3-17),  and 

kt  - CDpS/2mtFt  (6.3-29) 


Calling  the  initial  and  final  velocities,  VpCR  and 
Vot,  Eq.  (6.3-28)  is  integrated  to  obtain 


- V /(1+V  lu.t;  / 2) 
Dt  °CR  °CR^  Ft 


CR 


(6.3-30) 


and  the  corresponding  ranges,  Rq  and  R^, 
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(6.3-31) 


Next,  the  vehicle  must  decelerate  from  Vd^  to  the 
desired  initial  velocity  for  descent,  Vd^.  The  deceleration 
model  is 

V - -*dV2  (6.3-32) 

where 


kd  - CDpS/2m 


(6.3-33) 


Equation  (6.3-32)  is  integrated  once,  implicitly  defining 
the  deceleration  time,  ty^,  and  again  to  obtain  the  range 
at  the  end  of  deceleration,  Rg: 


(6.3-34) 

(6.3-35) 


The  range  increment  during  deceleration,  RWt  due  to  the 

prevailing  wind,  Vw  , is 

t 

RWt  - vwt(tFt+tFd>  (6.3-36) 


The  range  travelled  during  descent  is  obtained  from  the 
developments  of  Section  6.3.2.  The  wind-free  range  increment 
during  descent  is 


R3  " <h-hGO)(L/D>max  (6-3-37) 

assuming  that  the  current  altitude,  h,  and  the  altitude 
following  deceleration  are  the  same.  The  wind  effect  during 
descent  is  a function  of  the  descent  time  as  well  as  the  wind 
speed.  The  descent  (glide)  time,  tg,  can  be  approximated  by 
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t « _2_  fG"8h/2  _ -8hG0/2 

g kg8  6 > 


where 

k 


' 2Wcosy 


g 


921  sinYont 
S °Pt 


Jopt 


GO 


(6.3-38) 


(6.3-39) 


and  W is  the  average  value  of  aircraft  weight  during  the 
descent.  Variations  in  atmospheric  temperature  and  pressure 
can  be  reflected  in  the  values  of  scale  height,  8,  and  air 
density  at  the  destination  altitude,  p^q.  With  an  average 
wind  of  Vwg  during  descent,  the  wind-induced  range  increment 
is 


Vw  tg 
g 


(6.3-40) 


The  total  range  increment  from  descent 
to  destination  is 

R ■ R1  + R2  + "3  + % + Rwg 


initiation 

(6.3-41) 


and  the  descent  is  initiated  when  R > Rqq.  The  descent  range 
and  range  to  go  should  be  evaluated  at  frequent  intervals 
so  that  the  segment  switch  occurs  at  the  proper  time. 


Crew-Display  and  Throttle  Control  Interface  - Although 
throttle  settings  are  computed  (as  a function  of  the  available 
normal  rated  thrust),  the  principal  indicators  of  engine 
thrust  on  the  C-141A  are  engine  pressure  ratio  (EPR)  and 
low-speed  compressor  RPM  (N^).  For  the  crew  to  exercise 
manual  throttle/energy  management  or  for  an  automatic  control 
closure  to  operate  most  effectively,  the  proper  values  of  EPR  and 
should  be  computed.  (The  computations  are  optional  for  the 
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automatic  system,  as  the  throttle  trim  integrator  provides 
the  appropriate  steady-state  throttle  setting  with  or  without 
knowledge  of  EPR  and  N^.)  The  necessary  transformations  are 
presented  graphically  in  Ref.  32,  and  supplemental  infor- 
mation is  contained  in  Ref.  47.  The  two  transformations  can 
be  stored  as  functions  of  air  temperature  and  altitude, 
taking  anti-ice  and  other  engine  bleed  air  requirements  into 
account . 


6.4  SYSTEM  SYNTHESIS 

The  throttle/energy  management  structure  has  been 
defined  for  a fully  automatic  system,  as  shown  in  Fig.  6.3-1, 
but  the  figure  also  suggests  alternatives  for  manual  par- 
ticipation. There  are  four  levels  at  which  fuel-minimizing 
throttle/energy  management  can  be  implemented: 

• Handbook/Calculator  Flight  Path 
Computation  (Level  1) 

• Automatic  Flight  Path  Display 
(Level  2) 

• Automatic  Flight  Path  and  Command 
Displuy  (Level  3) 

• Automatic  Throttle/Energy  Management 
(Level  4) 

To  implement  Level  1,  the  functions  conducted  in 
the  Flight  Path  Generator  are  assembled  in  handbook  format, 
and  supplemental  computations  are  defined  for  a programm- 
able electronic  calculator.  The  crew  is  required  to  make 
pre-flight  calculations,  to  read  aircraft  instruments  in 
flight,  to  compute  desired  flight  variables,  and  to  follow 
Vp,  hp,  and  hp  with  standard  piloting  procedures.  As  in 
the  case  of  current  EPR  setting  (Ref.  47),  computations 
should  be  made  at5,000-ft  intervals  during  the  climb  until 
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the  climb  rate  has  decreased  below  1,000  ft/min,  when  more  fre- 
quent computations  are  warranted.  A similar  schedule  should 
be  followed  during  descent.  During  cruise,  computations 
should  be  made  at  30-min  intervals  (Ref.  35),  with  decreased 
intervals  as  the  descent  initiation  point  is  approached. 

The  need  for  data  lookup  and  computation  by  the 
crew  is  eliminated  in  Level  2 implementation,  as  the  throttle/ 
energy  management  computer  generates  Vp,  and  hp  from 
stored  data  and  crew  inputs.  The  computer  operates  without 
direct  input  of  flight  data  from  aircraft  sensors,  so  it  is 
necessary  for  the  crew  to  enter  (manually)  all  mission-specific 
data.  (Aircraft-specific  data  is  stored  within  the  com- 
puter. ) Computer  outputs  occur  on  a continuous  basis,  and 
they  are  followed  with  standard  piloting  procedures.  The 
Level  2 system  provides  a major  reduction  in  crew  workload, 
and  it  is  independent  of  all  aircraft  systems.  Conceptually, 
the  Level  2 system  can  be  contained  in  a standard  mini- 
computer which  is  supported  by  a standard  keyboard-input/ 
alphanumeric-output  terminal  device. 

The  Level  3 system  incorporates  feedback  information 
from  aircraft  sensors.  The  crew  workload  is  reduced  further, 
and  the  throttle/energy  management  system  displays  Vp,  hp, 
hp  ©c , and  Tc  (or  EPRC  and  N^c).  Crew  interaction  with  the 
Level  3 system  consists  of  mission  definition,  mode  switching, 
and  standard  piloting  procedures  (aided  by  suggested  settings 
for  pitch  angle  and  throttle,  which  are  computed  by  the  per- 
turbation guidance  law  and  the  trim  integration).  This 
implementation  automatically  performs  all  of  the  fuel- 
minimizing throttle/energy  management  functions,  but  it  leaves 
the  piloting  to  the  pilot.  The  Level  3 logic  is  contained  in 
a flight  computer,  which  is  supported  by  real-time  peripheral 
input  devices  as  well  as  the  crew  terminal. 
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The  Level  4 system  is  fully  coupled  to  autopilot/ 
autothrottle  actuators  (via  inner-loop  regulators),  although 
the  system  monitors  provide  a Level  3 capability  as  an  option. 
As  in  the  previous  case,  the  crew  must  provide  mission  and 
mode  information,  but  the  pilot  relinquishes  inner-loop  con- 
trol to  the  automatic  system  under  normal  flying  conditions. 

He  retains  full  authority  to  change  flight  level  and  Mach 
number  during  cruise,  to  step  climb,  to  command  alternate 
climb  and  descent  speeds,  to  modify  data  in  flight  (e.g. , to 
change  the  destination  location,  terminal  altitude,  or  air- 
speed, or  to  respecify  aircraft  weight  and  c.g.  location 
following  an  air  drop),  or  to  assume  total  control  in  an 
emergency. 

Although  not  considered  in  the  current  study,  the 
throttle/energy  management  algorithm  can  be  expanded  readily 
to  a 4-D  guidance  law.  Lateral  guidance  requires  the  addition 
of  heading  and  cross-track  information,  leading  to  roll-angle 
commands.  Time-slot  control  is  affected  by  altering  VD> 
Fuel-minimizing  guidance  is  totally  compatible  with  air  traffic 
control  and  with  the  stringent  requirements  of  time-based 
navigation. 

The  remainder  of  this  section  is  devoted  to  the 
requirements  for  implementation  of  an  experimental  Level  4 
throttle/energy  management  system,  including  crew  interface, 
sensors  and  actuators,  digital  flight  computer,  and  software 
development. 

6.4.1  Crew  Interface  Requirements 

At  a minimum,  the  Level  4 system  requires  an  alpha- 
numeric display  for  crew  monitoring  of  system  performance  and 
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a keyboard  for  mission  data  input,  mode  switching,  and 
experimental  procedures.  Manual  piloting  in  the  Level  3 
mode  would  be  aided  by  a graphical  or  flight  director  dis- 
play , and  various  control  functions  can  be  implemented  with 
dedicated  switches.  The  following  complement  would  be  use- 
ful in  an  experimental  system,  with  P indicating  location  at 
the  pilot's  station  and  E representing  the  flight  engineer's 
station : 

• CRT  display,  or  flight  director/LED 
display  (P) 

• Switch  panel  (P) 

• CRT  display  (E) 

• Full  keyboard  (E) 

An  example  of  a display/keyboard  panel  for  the 
pilot's  station  using  a cathode-ray  tube  (CRT)  and  calculator 
button-type  input  is  illustrated  in  Fig.  38.  Actual  and 
desired  values  of  engine  parameters  and  pitch  angle  are  dis- 
played on  the  CRT  screen  with  alphanumeric  presentation  of 
velocity,  altitude,  and  mode  information.  The  pilot  can 
enter  switching  and  numerical  information  through  the  key- 
board, and  special  system  status  information  (e.g. , failure 
indication,  request  for  response,  and  pending  mode  switching) 
can  be  displayed  on  annunciator  lights.  Pilot  inputs  are  dis 
played  on  the  CRT  and  are  not  incorporated  until  the  pilot 
verifies  their  correctness  and  pushes  "ENTER". 

6.4.2  Sensor  and  Actuator  Requirements 

The  throttle/energy  management  algorithm  requires 
data  for  determining  the  proper  nominal  flight  path,  for 
perturbation  guidance  and  trim,  and  for  gain  scheduling. 
Manually  entered  data  includes 

• Mission  range  and  terminal  conditions 
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Figure  38  Display/Keyboard  Panel  for  Throttle/ 
Energy  Management  System 

• Air  traffic  control  constraints 

• Initial  weight  and  c.g.  location 

• Pressure  altitude  calibration 

• Wind  data  (unless  derived  on-board 

from  navigation  data) 

• Experimental  parameters 

Data  measurements  required  for  primary  guidance  include 

• Range  to  go 

• Total  temperature 

• Static  pressure 

• Total  pressure 

or  corresponding  air  data  computer  outputs.  Measurements 
which  are  not  required  but  which  would  improve  performance 
of  the  primary  guidance  logic  include 

• Fuel  quantity 

• Fuel  flow  rate 

• EPR,  N^,  or  throttle  setting 
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• Attitude  reference 

• Inertial  velocity 

• Vertical  acceleration 

Section  6.2  demonstrated  the  use  of  angle  of  attack  for 

gain  scheduling,  but  this  requirement  can  be  eliminated  in 
the  flight  system. 

The  information  listed  above  must  undergo  analog-to- 
digital  conversion  at  a 10-  to  12-bit  quantization  level, 
which  corresponds  roughly  to  resolutions  of  0.1  to  0.025 
percent  of  full  scale.  While  this  is  more  than  adequate 
for  most  measurements,  difficulties  could  be  encountered  with 
the  derivation  of  altitude  from  quantized  static  pressure 
measurements,  particularly  at  the  higher  altitudes.  For 
example,  with  12-bit  pressure  quantization,  and  a maximum 
value  equivalent  to  sea  level  pressure,  a single  bit  is 
equivalent  to  an  increment  of  25  ft  at  an  altitude  of 
40,000  ft.  The  problem  is  solved  by  double-word  quantization 
in  the  present  example,  the  altitude  resolution  could  be 
improved  to  better  than  0.01  ft  by  using  a double  word  for 
static  pressure. 

There  are  no  explicit  actuator  requirements  for 
throttle/energy  management  system,  because  its  outputs  are 
assumed  tc  command  an  inner-loop  controller.  Ultimately,  the 
Level  4 system's  commands  would  reach  hydraulic  actuators 
for  control  surface  and  throttle  displacement,  e.g.,  those 
used  in  the  C-141A  Fly-By-Wire  Program  (Ref.  48). 

6.4.3  Computer  Requirements 

It  is  assumed  that  the  throttle/energy  management 
logic  would  be  implemented  in  a digital  flight  computer 
for  a number  of  reasons.  The  logic  itself  involves  a number 
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of  multiplications  and  nonlinear  functions,  and  it  requires 
precision,  room  for  growth,  flexibility,  low  cost,  and 
moderate  reliability.  Previous  experience  with  analog 
fuel  management  computers  indicates  that  because  they 
lack  these  traits,  their  value  is  greatly  diminished  in 
practice.  Furthermore,  the  throttle/energy  management  sys- 
tem is  compatible  with  navigation  system  requirements,  and 
the  logic  could  be  implemented  in  a standard  navigation 
computer  (given  sufficient  memory),  which  almost  certainly 
would  be  digital. 

For  maintainability,  hardware  and  software  system 
support,  and  reliability,  it  is  desirable  to  use  off-the- 
shelf  equipment.  The  energy  management  function  for  a 
large  jet  transport  does  not  pose  stringent  requirements  on 
size,  weight,  or  computer  speed,  but  if  the  system  is  to 
be  considered  for  fleet-wide  retrofit,  logistical  consid- 
erations are  paramount. 

For  an  experimental  Level  4 system,  it  is  desirable 
that  the  flight  computer  have  hardware  floating-point  logic 
and  mandatory  that  it  be  supported  by  a higher-order  alge- 
braic programming  language  (i.e.,  that  a program  compiler 
exists).  Countless  recent  surveys  illustrate  the  cost  savings 
which  a compiler  affords  in  software  development.  Since  the 
higher-order  language  allows  the  program  designer  to  speci- 
fy operations  as  engineering  equations,  the  full  benefit  of 
the  compiler  is  realized  only  if  the  scaling  of  these  opera- 
tions is  automatic.  This  is  most  readily  accomplished  by 
working  in  floating-point  format.  Although  solution  speed 
is  a secondary  consideration  for  fuel -minimizing  guidance, 
the  one-  to  two-order  of  magnitude  improvement  in  execution  time 
which  the  floating-point  hardware  provides  minimizes  guidance 
sequencing  problems  and  improves  the  possibility  of  expanding 
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the  logic  into  a full  4-D  navigation  system.  It  also  allows 
the  experimental  system  to  execute  off-line  monitoring  compu- 
tations which  can  aid  the  testing  program  by  providing 
"quick-look"  analysis  of  system  performance. 

It  is  estimated  thet  the  throttle/energy  management 
logic  can  be  implemented  with  a memory  of  less  than  8,192 
("8K")  16-bit  words.  Experience  shows  that  the  memory  pro- 
vided should  be  at  least  double  the  requirement,  not  only  to 
allow  for  contingencies  but  to  allow  software  development 
to  proceed  without  the  pressures  of  squeezing  coding  into 
the  last  few  percent  of  available  memory.  The  computer  also 
should  be  supported  by  mass  storage  for  program  loading, 
data  input,  and  logging  of  system  outputs.  Again,  because 
speed  is  seconds-  ' and  operational  convenience  is  important, 
a cartridge  or  cassette  tape  drive  is  suggested  for  this 
function. 


6.4.4  Software  Development 

The  throttle/energy  management  structure  pre- 
sented in  Section  6.3  can  be  used  as  a preliminary  specifi- 
cation for  software  development,  and  the  preparation  of 
functional  flow  charts  and  coding  from  this  section  should 
be  direct.  The  figures  of  this  section  illustrate  the 
functions  which  remain  to  be  accomplished  in  developing 
computer  code  for  throttle/energy  management. 

Figure  39,  adapted  from  Ref.  59,  divides  the  soft- 
ware development  into  four  phases: 

• Analysis  and  design 

• Implementation 

• Evaluation  and  procedures 

• Delivery  (for  flight  test) 
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Figure  39  Phases  of  Software  Development 

The  first  phase  translates  engineering  requirements  into 
program  specifications  and  structure.  The  second  phase 
produces  coding  for  the  flight  computer.  The  third  phase 
involves  acceptance  testing  and  the  definition  of  crew 
procedures,  and  the  fourth  phase  provides  system  support. 

The  activities,  documentation,  and  milestones  which 
characterize  this  development  plan  are  summarized  in  Fig. 

40.  It  is  important  to  identify  these  elements  at  an 
early  stage  of  development  and  to  observe  certain 
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formalities  in  design  and  documentation  (e.g.,  anomaly 
reports,  change  notices,  and  test  reports)  even  for  small 
programs.  Efficient  management  of  available  resources  and 
assurance  that  milestones  are  accomplished  requires  visi- 
bility and  control  of  the  development  program. 

Development  of  fuel-minimizing  throttle/energy 
management  currently  is  in  the  first  (Analysis  and  Design) 
phase.  There  are  elements  of  this  phase  which  remain  to 
be  conducted,  including  the  preparation  of  flow  charts  and 
preliminary  design  review.  The  structure  of  the  Level  4 sys- 
tem is  defined,  the  preliminary  software  specifications  are 
ready  for  flowcharting  and  coding,  and  the  expected  perfor- 
mance of  the  system  has  been  estimated. 


6.5  CHAPTER  SUMMARY 

This  chapter  has  presented  the  basic  information 
which  is  necessary  to  design  and  develop  a fuel-minimizing 
guidance  system  for  the  C-141A  aircraft.  The  throttle/energy 
management  algorithms  consist  of  a flight  path  generator 
which  defines  a fuel-minimizing  reference  flight  path,  a 
perturbation  guidance  and  trim  law  which  automatically  con- 
trols the  aircraft  to  follow  the  reference  path,  and  gain 
scheduling  functions  which  assure  acceptable  guidance  response 
at  all  flight  conditions.  The  system  is  synthesized  for 
two-dimensional  (vertical-plane)  flight  paths  (Fig.  41)  and 
can  be  extended  to  include  lateral  guidance  and  time-based 
navigation. 

Four  levels  of  system  implementation  are  discussed. 
Handbook/Calculator  Flight  Path  Computation  (Level  1)  requires 
no  special  equipment , but  it  imposes  a large  workload  on  the 
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crew,  and  it  is  likely  to  save  less  fuel  than  the  other  levels. 
Automatic  Flight  Path  Display  (Level  2)  requires  a minicomputer 
for  data  storage  and  computation,  and  it  eliminates  the 
requirements  for  many  manual  functions.  The  computer  is 
independent  of  aircraft  systems,  relying  on  the  crew  for 
all  inputs,  and  it  could  be  moved  form  one  aircraft  to  another 
with  relative  ease.  Automatic  Flight  Path  and  Command  Display 
(Level  3)  uses  direct  inputs  of  aircraft  data  to  compute  a 
continuous  output  of  state  and  control  variables,  minimizing 
data  input  by  the  crew  and  allowing  the  pilot  to  follow  the 
fuel-optimal  flight  path  as  closely  as  he  deems  necessary. 
Automatic  Throttle/Energy  Management  (Level  4)  closes  the 
guidance  loop  through  autopilot /autothrottle  actuators, 
relying  on  the  crew  to  perform  executive  functions  only 
(mission  specification  and  choice  of  optional  flight  param- 
eters). 


Concepts  for  testing  throttle/energy  management  on 
operational  aircraft  are  presented  in  the  following  chapter. 
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VII  . 


FLIGHT  TEST  CONCEPTS 


7.1  OVERVIEW 

The  confirmation  of  fuel-saving  flight  policies  by 
testing  in  the  actual  aircraft  is  made  difficult  by  uncer- 
tainties in  the  environment  and  in  the  aircraft  itself.  To 
achieve  the  goal  of  reducing  fuel  use  in  practice  and  on  a 
fleet-wide  basis,  it  is  necessary  to  consider  operational 
constraints  and  preferences  during  test  planning.  The  cost 
of  testing,  as  well  as  the  fuel  consumed  to  provide  conclusive 
results,  also  should  be  considered. 

The  primary  objectives  for  flight  testing  are  the 
following : 

• Demonstrate  fuel-optimal  flight  paths 

• Evaluate  throttle/energy  management 
technique  at  four  levels  of  implementation. 

• Collect  data  regarding  air  crew  and  air 
traffic  controller  acceptance  of  fuel- 
saving flight  policies. 

• Estimate  the  probable  impact  of  fuel- 
minimizing throttle/energy  management 
on  Air  Force  fuel  consumption  and 
operating  costs. 

• Determine  the  level  of  throttle/energy 
management  which  is  most  likely  to  pro- 
vide the  greatest  benefits  to  Air 
Force  flight  operations. 

This  chapter  presents  a plan  for  evaluating  throttle/ 
energy  management  in  the  face  of  the  uncertainties  and  limi- 
tations which  are  mentioned  above.  There  are  two  basic 
tenets  of  the  suggested  flight  test  program: 
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• Fuel-saving  flight  policies  should 

be  demonstrated  on  line  aircraft  (using 
Level  1-2  implementation)  and  on  an 
aircraft  equipped  with  Level  3-4  equip- 
ment . 

• Evaluation  of  test  results  will  be 
aided  materially  by  computerized  data 
analysis  techniques,  described  below. 

The  capital  equipment  required  for  Level  1-2  implementation 
is  low;  therefore,  test  results  can  be  obtained  during 
normal  operations  and  at  low  cost.  Nevertheless,  fuel- 
minimizing guidance  requires  Level  3-4  instrumentation, 
which  should  be  installed  on  at  least  one  test  aircraft. 
Variations  in  atmospheric  conditions,  aircraft  weight, 
etc. , will  make  a direct  comparison  of  fuel  used  on  two 
flights  (for  example,  an  optimal  flight  path  and  a conven- 
tional flight  path)  difficult;  therefore,  a three-way  tra- 
jectory analysis  procedure  is  developed  to  eliminate  uncer- 
tainties in  data  evaluation. 

Experiment  descriptions,  equipment  and  data 
requirements,  and  data  evaluation  are  described  in 
Sections  7.2,  7.3,  and  7.4,  respectively.  The  chapter  is 
summarized  in  Section  7.5. 

7.2  EXPERIMENT  DESCRIPTIONS 

A series  of  tests  is  suggested  in  this  section. 

The-  tests  fall  into  two  categories:  those  which  use  minimal 

equipment  (Levels  1 and  2)  and  are  intended  to  identify 
operational  use  factors,  and  those  which  use  more  complete 
throttle/energy  management  equipment  (Levels  3 and  4)  and 
are  devoted  to  system  development . The  first  category  of 
tests  can  be  conducted  on  line  aircraft  performing  scheduled 
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exercise  and  airlift  operations.  The  second  category  of 
tests  requires  the  installation  of  aircraft-computer  inter- 
faces and  an  input-display  device  at  the  pilot's  and  engineer’s 
stations.  This  system  ultimately  is  tested  in  line  operations, 
although  it  must  be  available  for  flights  devoted  exclusively 
to  system  tests,  as  well. 

7.2.1  Tests  of  Operational  Use  Factors 

Six  test  routes  should  be  chosen  for  these  tests,  two 
each  in  the  short-,  medium-,  and  long-range  categories  defined 
in  Chapter  5.  Where  possible,  the  routes  should  be  flown 
regularly  by  MAC  personnel  in  both  directions,  allowing 
head-tail  wind  and  flight  level  effects  to  be  evaluated. 

They  should  be  in  areas  of  relatively  light  commercial  air 
traffic,  so  that  clearances  for  special  flight  paths  can  be 
obtained,  as  required.  Conducting  tests  at  night  will  increase 
the  flexibility  in  choosing  flight  paths  by  reason  of  lighter 
traffic.  The  routes  should  be  in  areas  of  good  Air  Weather 
Service  coverage,  in  order  that  meteorological  data  required 
for  data  analysis  is  accurate  and  available. 

It  is  preferable  that  a limited  number  of  aircraft 
(4  or  less)  be  used  in  the  testing  program  and  that  up  to 
six  flight  crew  members  (in  each  rating)  conduct  the  tests. 
Characteristics  of  each  aircraft  (by  tail  number)  and  its 
engines  should  be  documented,  including  time  since  major 
overhaul,  minor  performance-related  maintenance  actions, 
and  known  deviations  from  standard  configuration  (e.g. , 
actual  thrust  and  specific  fuel  consumption  measured  at 
last  engine  overhaul,  throttle  and  control  surface  rigging 
anomalies,  etc.).  Each  aircraft  need  not  fly  each  route 
segment,  although  it  is  desirable  for  each  segment  to  be 
flown  by  more  than  one  aircraft.  Flight  crew  members  should 
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provide  a representative  range  of  experience  levels  and 
should  be  able  to  adapt  to  non-standard  procedures. 

Two  test  engineers  accompany  each  flight.  One 
engineer  is  responsible  for  executing  the  Level  1-2  throttle/ 
energy  management  technique  and  transmitting  flight  path 
information  to  the  crew.  The  other  engineer  is  responsible 
for  logging  data,  verifying  crew  procedures,  and  recording 
the  crew's  comments  on  fuel-saving  procedures.  Both  engi- 
neers should  develop  a summary  of  general  use  factors  after 
each  flight. 

Ten  to  twenty  flights  should  be  conducted  in  each 
of  the  three  major  route  lengths,  with  the  flights  being 
divided  approximately  equally  between  fuel-saving  and  con- 
ventional procedures.  The  former  are  flown  to  evaluate 
potential  fuel  savings,  while  the  latter  provide  reference 
data  plus  sensitivities  to  major  variables,  e.g.,  cruise 
altitude  and  speed.  Each  chosen  route  segment  should  be 
flown  at  least  twice:  once  optimally  and  once  con*  entiona  ly. 

Accumulated  flying  time  (by  route  length)  is  approximately 

• Short  Range:  10  to  25  hours 

• Medium  Range : 25  to  75  hours 

• Long  Range : 75  to  200  hours 

for  Level  1-2  testing.  Tests  should  extend  over  a period 

of  six  to  nine  months  to  allow  time  for  data  evaluation  and 
incorporation  of  early  results  in  final  testing.  It  is 
recommended  that  each  test  flight  be  conducted  as  a unified 
experimental  profile  containing  climb,  cruise  (if  required), 
and  descent  segments  that  (together)  form  an  optimal  or  con- 
ventional "hop,''  although  individual  flight  segments  can  be 
analyzed. 


Examples  of  flight  types  are  given  in  Table  18, 
which  shows  15  profile  variations  (10  conventional,  5 optimal). 
The  optimal  profiles  are  flown  using  Level  1 and  2 throttle/ 
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energy  management;  hence,  the  flight  types  are  equally 
divided.  In  practice,  it  may  be  necessary  to  repeat  tests, 
and  all  15  flight  types  need  not  be  conducted  in  every  range 
category. 

TABLE  18 

FLIGHTS  FOR  OPERATIONAL  USE  TESTING 


Trajectory  Type 

(eight 

Climb 

Crulae 

Deeceat 

Cemmeate 

roeventloeal 

Light 

Staadard 

Staadard 

Staadard 

Crulae  Mach  Humber  - 0.74 

N 

Medium 

N 

N 

N 

H 

H 

Heavy 

H 

•e 

N 

M 

Optimal 

Light 

Optimal 

Optimal 

Optimal 

Levele  1 aad  2 

Medium 

H 

H 

N 

Heavy 

ee 

N 

M 

Coavaatioaal 

Medium 

Staadard 

Staadard 

Staadard 

Crulae  Mach  Humber  • 0.767 

N 

H 

IAAtSO  KIAH 

SS 

»• 

H 

•9 

Staadard 

N 

N 

N 

H 

M 

M 

Idle 

Crulae  Altitude  t 4000  ft 

H 

H 

M 

H 

Peaetratloa 

M 

H 

St 

H 

Staadard 

Head  vlad 

H 

H 

H 

H 

H 

Tall  wind 

Optimal 

" 

Optimal 

Optimal 

Optimal 

Tall  wind;  Levele  1 aad  2 

M 

H 

M 

H 

H 

Tall  vlad:Le«ele  1 aad  2 

Conventional  profiles  should  be  planned  and  executed 
using  standard  procedures.  It  may  be  desirable  for  these 
reference  paths  to  be  flown  by  pilots  who  are  not  familiar 
with  the  throttle/energy  management  program,  in  order  to 
get  an  unbiased  measure  of  fuel  consumption  on  truly  con- 
ventional profiles  (or  these  flights  can  be  conducted 
before  the  optimal  flights).  Repetition  of  optimal 
flight  paths  may  improve  fuel  conservation  by  giving  the 
pilots  the  opportunity  to  learn  differences  in  procedures. 

7.2.2  Tests  for  System  Development 

Flight  tests  for  the  development  of  a Level  3-4 
system  are  similar  to  the  Level  1-2  tests,  but  a greater 
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emphasis  is  placed  on  calibration  of  the  aircraft.  Con- 
sequently a single  aircraft  is  used  for  system  development 
Details  of  fuel  consumption  along  the  flight  path,  as  well 
as  end-to-end  consumption,  are  important.  In  addition, 
system  checkout , i.e.,  the  verification  of  proper  operation 
and  the  validation  of  basic  concepts,  is  a major  testing 
goal.  In  this  respect,  Level  3-4  testing  is  directed  at 
familiar  goals  of  instrument  and  control  system  development 

The  major  objectives  of  system  checkout  are 
evaluation  of 

• Fuel-minimizing  flight  path  generation 

• Perturbation  guidance 

• Mode  switching 

• Command  step  response 

• Disturbance  rejection 

• Human  factors 

• Sensor-computer  interfaces 

• Installation  and  maintenance  requirements 

Flight  time  devoted  to  system  checkout  can  be 
minimized  by  judicious  use  of  the  flight  computer  to  adapt 
to  self-testing  procedures,  e.g.  , injecting  artificial 
disturbances,  biasing  altitude  inputs  to  simulate  cruise, 
and  accelerated  sequencing  through  flight  path  segments. 
Prior  to  actual  fuel-use  testing,  checkout  functions  should 
be  completed  with  the  equivalent  of  single  short-range  and 
medium-range  flight  paths  — about  5 to  10  flight  hours. 

Fuel -use  testing  of  the  Level  3-4  system  can  be 
divided  into  three  phases: 

• Concept  Validation  (Phase  1) 

• MAC  Operational  Testing  (Phase  2) 

• Integrated  Testing  with  4-D 
Navigation  (Phase  3) 


165 


Concept  Validation  extends  for  approximately  60  flight  hours, 
allowing  the  throttle/energy  management  techn.  »ue  to  be 
demonstrated  and  evaluated  over  a wide  variety  of  mission 
ranges,  air  traffic  control  constraints,  vehicle  weights,  and 
atmospheric  conditions.  To  the  extent  possible,  testing  in 
this  phase  is  combined  with  other  test  programs  or  MAC 
exercise  flights.  During  MAC  Operational  Testing,  the  sys- 
tem is  flown  by  Airlift  Command  pilots  during  up  to  200  hours 
of  scheduled  flying.  Operational  use  and  maintenance  factors 
are  evaluated  during  this  phase,  and  the  system  is  demonstrated 
for  appropriate  personnel.  Whereas  Phase  1 is  supported  by 
detailed  data  analysis,  numerical  analysis  is  deemphasized 
in  Phase  2,  when  operational  factors  are  most  important. 

For  Integrated  Testing  with  4-D  Navigation,  the  system  is 
returned  to  research  status. 


7.3  EQUIPMENT  AND  DATA  REQUIREMENTS 

In  addition  to  the  system  synthesis  requirements 
discussed  in  Chapter  6,  flight  test  equipment  and  reference 
data  are  required.  For  all  tests,  synoptic  meteorological 
data  should  be  obtained  from  the  Air  Weather  Service. 

Standard  forms  used  by  the  Airlift  and  Logistics  Commands  to 
document  each  flight  stage  that  is  part  of  the  test  program 
should  be  obtained.  This  information  is  valuable  for 
statistical  analysis  of  test  results  and  for  calibration  of 
the  computer  analysis  described  in  the  next  section. 

No  additional  instrumentation  is  required  for  Level 
1-2  testing  (other  than  the  calculator  or  computer  used  for 
throttle/energy  management).  Test  engineers  record  necessary 
position  and  velocity  from  standard  air  data  and  navigation 
instrumentation  for  later  use  in  approximate  reconstruction  of 
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the  aircraft's  flight  path.  Standard  flight  recorder  infor- 
mation could  be  useful  in  this  regard. 


The  Level  3-4  aircraft  should  be  equipped  with 
precision  navigation  and  fuel-use  instrumentation  during 
Phase  1 testing  to  allow  accurate  reconstruction  of  the 
flight-path,  thrust,  and  fuel  consumption  profiles.  Engine 
condition  monitoring  equipment  could  improve  calibration  of 
the  power  plants. 


7.4  DATA  EVALUATION 

Because  no  two  flight  profiles  can  be  evaluated 
under  identical  vehicle  and  atmospheric  conditions  during 
flight  testing,  a direct  comparison  can  only  be  made  compu- 
tationally.  For  example,  if  a given  aircraft  flies  from  A 
to  B optimally  and  then  returns  to  A on  a conventions  path, 
the  winds  are  different  and  vehicle  weight  is  decreased. 

Two  aircraft  leaving  A at  the  same  time  to  fly  to  B are 
subject  to  performance  differences  on  the  order  of  5 to  10 
percent  due  to  differences  in  engine  wear,  control  rigging, 
payload  weight  and  distribution,  and  many  other  items  (Ref. 
60).  As  the  fuel  savings  which  the  testing  is  supposed  to 
demonstrate  are  of  the  same  order,  a single  comparison  may 
not  be  conclusive.  The  test  could  be  repeated  with  the 
roles  of  the  two  aircraft  reversed,  but  the  desired  result 
is  obtained  only  by  expending  twice  as  much  fuel  and  flight 
time  and  still  is  subject  to  atmospheric  uncertainties. 

Substantially  more  information  can  be  obtained  from 
each  flight  if  the  data  is  used  in  a three-way  numerical 
comparison,  consisting  of  computed  profiles  of  the 


167 


• Optimal  flight  path 

• Throttle/energy  management  flight  path 

• Conventional  flight  path 

There  are  two  steps  in  the  numerical  comparison : 

• Model  calibration  to  flight-test  data 

• Three-way  comparison  of  flight  paths 

The  process  is  illustrated  in  Fig.  42,  assuming  that  the 
flight  test  result  was  obtained  with  the  Level  4 throttle/ 
energy  management  system  engaged. 

MNH 


•I  MODEL  CALIBRATION  b)  THREE-WAY  COMPARISON 

Figure  42  Computer  Analysis  of  Flight 
Test  Data 


The  basic  approach  uses  a flight  path  simulation 
and  optimization  program  such  as  ENERGY,  described  in 
Appendix  A.  The  flight  test  data  includes  sampled  time 
histories  of  the  state  and  control  variables.  If  the  air- 
craft mathematical  model  (Appendix  B)  was  precise,  if 
atmospheric  conditions  were  known,  and  if  the  measurements 
contained  no  errors,  the  flight  path  could  be  reconstructed 
exactly  by  integrating  the  equations  of  motion, 

x(t)  - f I(x(t),u(t)J  (7.4-1) 
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using  the  measured  control  history  as  an  input;  however, 
errors  and  uncertainties  preclude  an  exact  solution.  The 
best  that  can  be  expected  is  that  the  reconstructed  flight 
path  minimizes  an  error  criterion  in  fitting  the  data.  For 
example,  let  the  measured  state  be  3^  and  the  computed  state 
be  x^,.  Forming  the  error  vector 


Ax  - 


in, 

—TO 


— C 


(7.4-2) 


a set  of  parameters,  pfc,  is  chosen  to  minimize  a quadratic 
cost  function  of  the  fitting  error: 


J 


T 

Axf  Pf  Axf 


Ax^(t)  AAx(t)dt 


(7.4-3) 


as  in  Section  3.5.  In  other  words,  the  data  fitting  problem 
is  formulated  as  a parametric  optimization,  and  the  param- 
eters are  calibration  coefficients  of  the  aircraft's  mathe- 
matical model.  The  determination  of  these  coefficients 
follows  the  discussion  of  Section  3.5  and  is  illustrated  in 
Fig.  A. 3-1. 


It  is  important  to  identify  a minimal  set  of  cali- 
bration coefficients,  and  in  the  present  case,  it  appears 
that  the  minimal  number  is  three.  Referring  to  Eq.  (3.2-1) 
to  (3.2-7),  it  can  be  seen  that  the  aircraft  model  enters  four 
equations  ($,Y.£,m)  but  the  thrust  and  aerodynamic  terms 
are  identical  in  two  equations  (y,£).  For  the  purposes  of 
comparing  fuel  used  on  various  trajectory  types,  it  can  be 
assumed  that  the  form  of  the  aircraft  model  (Appendix  B)  is 
known  and  that  the  major  variations  from  one  flight  to  the 
next  are  contained  in  the  relative  scaling  of  the  V,  y,  and 
m equations.  Three  parameters  then  are  necessary  to  describe 
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• Thrust -drag  scaling 

• Li ft -weight  scaling 

• Fuel  flow  rate  scaling 

and  these  parameters  account  for  dynamic  variations  due  to 
engine  parameters,  weight,  aerodynamics,  and  atmospheric 
effects.  The  disadvantage  of  this  calibration  procedure 
is  that  these  parameters  do  not  identify  the  source  of 
variations.  The  advantage  is  that  they  provide  a simple 
calibration  of  the  aircraft  model  for  comparing  fuel-optimal 
and  conventional  flight  operations. 

The  best  definitions  tor  the  scaling  parameters 
are  subject  to  confirmation;  however,  the  following  set 
appears  reasonable:  the  coefficients  C^,  C^,  and  m are 

replaced  by 


°D 

-P1CD 

(7.4-4) 

<£ 

' P2CL 

(7.4-5) 

* 

0) 

■ P3“ 

in  the  equations  of  motion,  and  the  p^  are  chosen  to  minimize 
Eq.  (7.4-3),  calibrating  the  computer  model  (as  shown  in 
Fig.  42a)  and  completing  the  first  step. 

Program  ENERGY  is  used  to  compute  the  remaining  two 
flight  paths  in  the  three-way  comparison  using  the  cali- 
brated model.  If  the  flight  test  used  throttle/energy 
management,  an  optimal  trajectory  and  a conventional  flight 
path  are  computed  (Fig.  42b);  if  the  test  path  was  con- 
ventional, then  the  optimal  and  throttle/energy  management 
paths  are  computed.  The  flight  test  program  allows  the 
validity  of  this  comparison  to  be  assessed,  as  both 
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conventional  and  throttle/energy  management  results  can  be 
used  in  model  calibration.  A single  flight  is  thus  capable 
of  producing  results  which  are  equivalent  to  four  test 
flights  using  the  two-aircraft  approach  mentioned  at  the 
beginning  of  this  section.  This  data  analysis  technique 
can  be  used  with  either  Level  1-2  or  Level  3-4  flight  tests 
and  with  any  number  of  subject  aircraft,  because  the  com- 
puter model  is  recalibrated  for  each  comparison. 


7.5  CHAPTER  SUMMARY 

Objectives  for  flight  testing,  descriptions  of 
experiments,  data  requirements,  and  a new  technique  for 
flight  path  evaluation  have  been  presented  here.  It  is 
shown  that  tests  should  demonstrate  not  only  the  technical 
characteristics  of  throttle/energy  management  techniques  but 
the  operational  use  factors  as  well.  Real  fuel  savings  will 
be  achieved  only  if  fuel-saving  systems  are  acceptable  to 
air  crews.  Testing  of  essentially  manual  (Level  1-2)  tech- 
niques in  a MAC  exercise-airlift  environment  requires  about 
100  to  300  flight  hours.  Testing  of  automatic  (Level  3-4) 
throttle/energy  management  requires  an  equivalent  amount  of 
time.  Evaluation  of  the  true  potential  for  saving  fuel  in 
flight  operations  is  assisted  by  using  a self-calibrating, 
three-way  computer  comparison  method  which  minimizes  flight 
testing  requirements. 
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VIII. 


CONCLUSION 


The  objectives  of  this  investigation  have  been  to 

a Determine  fuel  savings  which  can  be 

realized  in  the  flight  operations  of 
the  C-141A  aircraft. 

• Synthesize  a throttle/energy  management 
technique  which  is  capable  of  guiding 
the  C-141A  aircraft  along  fuel-optimal 
flight  paths. 

• Establish  methods  for  demonstrating 
throttle/energy  management  techniques  by 
flight  testing. 

• Identify  fuel  conservation  principles 
which  could  be  applied  to  other  jet 
transport /bomber  aircraft. 


In  the  course  of  conducting  this  work,  the  fuel- 
optimal  performance  of  the  C-141A  was  defined  for  a wide 
variety  of  vehicle  and  atmospheric  conditions;  for  climb, 
cruise,  and  descent  segsmnts;  and  for  fully  integrated 
flight  oaths  with  ranges  of  200  to  5,600  nm.  These  numerical 
results  were  obtained  using  a newly  developed  computer  pro- 
gram (ENERGY)  which  generates 

e Conventional  flight  profiles 
e Optimal  flight  profiles 

e Flight  profiles  with  closed-loop 

guidance 

a Linear-optimal  gain  values  for  an 

automatic  guidance  (combined  autopilot/ 
autothrottle)  system 

This  program  incorporates  a detailed,  nonlinear  model  of  the 
C-141A  aircraft,  as  well  as  functions  for  duplicating  con- 
straints due  to  air  traffic  control  and  vehicle  limitations. 
It  uses  new  techniques  for  efficient  computation  of  flight 
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paths  and  rapid  convergence  of  optimal  solutions.  The  pro- 
gram can  be  applied  to  other  aircraft  by  changing  aerodynamic, 
thrust,  and  inertial  data.  For  Jet  transport /bomber /util ity 
aircraft,  e.g. , C-5A,  B-52,  KC-135,  YC-14,  YC-15,  B-l,  T-39, 
and  wide-bodied  or  conventional  commercial  jet  transports, 
data  conversion  would  be  a simple  and  direct  process.  Fur- 
thermore, the  program' 8 flexibility  could  be  of  substantial 
benefit  in  comparative  flight  testing  of  fuel  conservation 
techniques,  reducing  flight-hour  requirements  by  up  to  75 
percent . 


The  throttle/energy  management  technique  developed 
here  has  four  major  elements: 

• Fuel -optimal  flight  path  generator 
e Perturbation  guidance  and  trim  law 

• Guidance  gain  scheduling 

• Crew  input  and  display 

The  technique  can  be  implemented  at  four  levels  of  system 
complexity,  ranging  from  essentially  manual  to  fully  auto- 
matic operation: 

• Handbook/Calculator  Flight  Path 
Computation 

• Automatic  Flight  Path  Display 

• Automatic  Flight  Path  and  Command 
Display 

• Automatic  Throttle/Energy  Management 

Fuel-minimizing  throttle/energy  management  is  shown  to  be 
compatible  with  conventional  flight  operations  at  the  first 
level  and  with  integrated  4-D  navigation  systems  at  the 
fourth  level. 

The  suggested  flight  test  program  is  divided  into 
two  parts:  the  evaluation  of  operational  use  factors  and 

system  development.  Level  1 and  2 testing  falls  into  the 
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first  category,  and  it  is  conducted  largely  by  MAC  personnel. 
Level  3 and  4 falls  into  the  second  category,  which  is  fur- 
ther subdivided  into 

• Concept  Validation 

• MAC  Operational  Testing 

• Integrated  Testing  with  4-D  Navigation 

Uncertainties  and  measurement  errors,  which  make 
comparative  evaluation  of  fuel  use  difficult  (because  the 
levels  of  fuel  saving  and  uncertainty  are  themselves  com- 
parable), are  largely  overcome  by  a three-way  numerical 
comparison  of  flight  data.  The  computer  model  is  calibrated 
by  the  flight  test  data,  and  the  analysis  program  is  used  to 
generate  three  flight  paths  which  are  dynamically  equivalent 
to  the  measurements: 

e Optimal  Flight  Path 

• Throttle/Energy  Management  Flight  Path 

e Conventional  Flight  Path 

It  is  easy  to  overstate  the  potential  fuel  and 
cost  savings  associated  with  throttle/energy  management 
because  there  is  so  much  variability  in  "conventional" 
flight  paths.  A high-speed,  low-altitude  cruise  is  known  to 
use  more  fuel  than  a Breguet  cruise-climb;  potential  savings 
relative  to  the  first  path  are  impressive,  but  they  may  be 
small  compared  to  the  second  path.  Actual  fuel  savings  depend 
on  actual  flight  policies.  Taking  the  optimal  flight  profile 
as  a base,  this  study  has  demonstrated  that  some  handbook" 
profiles  use  over  30  percent  more  fuel  than  the  optimum,  while 
proper  choice  of  cruise  altitude  and  use  of  idle-thrust  , 
(L/!>)max  descent  can  narrow  the  gap  to  a few  percent.  The  key 
to  minimizing  fuel  use  on  military  Jet  transport  is  to  make 
fuel  saving  the  "natural*’  way  to  fly  the  airplane.  A fuel- 
minimizing throttle/energy  management  system  can  achieve  this 
Koal . 
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APPENDIX  A 


COMPUTER  TOOLS  FOR  TRAJECTORY  ANALYSIS 


A. 1 PROGRAM  DEVELOPMENT 

In  order  to  develop  and  analyze  potential  fuel  reducing 
throttle/energy  management  systems  for  aircraft,  a computer 
program,  ENERGY,  has  been  developed.  Program  ENERGY  provides 
the  capability  to  exercise  a wide  variety  of  trajectory  shap- 
ing and  analysis  techniques  and  to  validate,  by  simulation, 
any  proposed  on-board  guidance  algorithms.  A list  of  the 
major  applications  of  program  ENERGY  include 

e Aircraft  trajectory  generation 

e Parametric  trajectory  optimization 

e Variational  ( steepest-descent ) trajectory 
optimization 

e Linear-optimal  perturbation  guidance 
gain  computation 

e Comparative  evaluation  of  varying  trajec- 
tory shapes  and/or  on  board  gjidance 
algorithms . 

The  modular  structure  of  program  ENERGY  (Fig.  4-1) 
facilitates  its  modification  for  other  types  of  trajectory 
analysis  and  enables  each  module  to  be  used  in  various 
applications.  For  example,  both  variational  optimization  and 
linear-optimal  guidance  gain  computations  necessitate  evaluat- 
ing partial  derivatives  of  the  dynamic  equations  of  motion 
along  the  aircraft’s  trajectory.  By  programming  the  partial 
derivative  computations  as  a separate  program  module,  dupli- 
cation of  effort  is  eliminated  and  program  checkout  and  vali- 
dation is  simplified. 


17fS 


•- 14791 


Figure  A-l  Program  Structure  of  the 
Engineering  Simulation 
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The  advantages  of  ENERGY'S  modular  programming  structure 
are  also  demonstrated  by  the  programs  guidance  algorithm  options 

• Open-Loop  Trajectories  - Guidance  commands 
are  specified  as  a function  of  an  inde- 
pendent variable,  e.g.,  time  or  range. 

• Ad  Hoc  Guidance  - Guidance  commands  are 
generated  based  upon  arbitrary  perturbation 
gains. 

• Optimal  Perturbation  Guidance  - Guidance 
commands  are  generated  using  optimally  de- 
rived perturbation  gains. 

These  guidance  options  are  useful  for  generating  ref- 
erence trajectories,  performing  trajectory  optimization,  and 
comparing  guidance  algorithms.  These  guidance  techniques 
are  by  no  means  an  exhaustive  set  of  possible  algorithms; 
however,  the  modularity  of  program  ENERGY  simplifies  the 
addition  of  any  options  which  might  be  of  interest. 

Those  sections  of  ENERGY  which  define  the  vehicle 
aerodynamics  also  are  programmed  as  self-contained  modules. 

The  model  of  the  C-141A  aircraft,  which  was  the  subject  of 
the  current  study,  is  presented  in  Appendix  B.  Other  air- 
craft can  be  analyzed  by  simply  replacing  the  C-141A  model 
subroutines  with  those  that  correspond  to  the  vehicle  of 
i nterest . 


A. 2 TRAJECTORY  SHAPING  AND  EVALUATION 

Trajectory  optimization  consists  of  determining  a tra- 
jectory and  control  profile  which,  in  some  sense,  is  "best". 
The  approach  used  in  ENERGY  to  compare  different  trajectories 
is  to  assign  a scalar  cost  to  each  trajectory.  Typically,  a 
general  trajectory  cost  function  has  the  form 
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I 


where 


(A. 2-1) 


J ■ scalar  cost  value 

x(R. ) ■ final  values  of  trajectory  state 
variables  (velocity,  flight  path 
angle,  altitude,  mass,  and  time) 

u - control  variables  (pitch-angle, 
throttle,  and  roll) 

RQ,R_  ■ initial  and  final  values  of  inde- 
pendent variable  (range) 

L ■ scalar  integrated  cost  function 

4>  * scalar  function  of  final  state 
variables 

The  optimal  or  "best"  trajectory  is  the  one  which  provides 
the  lowest  cost  value. 


The  inclusion  of  numerous  penalty  functions, 
terminal  conditions,  *(x,(Rj)),  and  in-flight  constraints, 

L(x,  u,  RX  in  program  ENERGY  provides  a wide  degree  of  flexi- 
bility in  shaping  and  evaluating  trajectory  and  control  pro- 
files. A summary  of  the  penalty  functions  which  can  be 
specified  in  ENERGY  is  presented  in  Table  A-l . 


A simple  example  demonstrates  how  these  penalty 
functions  can  be  used  for  trajectory  optimization.  If  it  is 
desired  that  a trajectory  terminate  at  a particular  altitude, 
Hp,  and  not  exceed  a maximum  Mach  number,  MACH|(AX,  while 
minimizing  fuel  use,  Am,  the  appropriate  cost  function  to 
minimize  would  be  given  as 
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TABLE  A-l 


TRAJECTORY  PENALTY  FUNCTIONS 


Terminal  Constraint* 
•<*(*,)) 

State  Constraint* 
I'x.u.R) 

Control  Constraint* 
t(x.u.R) 

Altitude 

Flight  Path  Angle 
Max* 

Velocity 

Tie* 

Maximum  Altitude 
Minimum  Altitude 
Xxcluded  Altitude  Bead 
Maximum  Altitude  Rate 
Maxlmua  Mach  Rumber 
Buffet  Limit* 

Maximum  IAS 
Minimum  IAS 

Maximum  Throttle 
Minimum  Throttle 
Maximum  Pitch  Angle 
Minimum  Pitch  Angle 
Maximum  Angle-of-AttacL' 
Minimum  An.jle-of-Attack 

J * C1(HD-H(Rf))2  + CgUnO+J*1,  L(MACH)  dR  (A. 2-2) 


where 


(C3  (MACHMAX-IIACH)  , if  MACH  > MACH^ 


, if  MACH  < MACH, 


MAX 


(A. 2-3) 


and  C^,  Cg,  are  positive  scalars 


Examination  of  the  cost  function  of  Eq.  (A. 2-2) 
shows  an  increase  in  cost  for  missing  the  desired  terminal 
altitude,  for  usincr  a large  amount  of  fuel,  Am,  or  for 
exceeding  the  maximum  Mach  number  along  the  trajectory. 


The  penalty  functions  listed  in  Table  A-l  show 
that  program  ENERGY  can  describe  the  characteristics  of  a 
desired  trajectory  in  some  detail.  In  keeping  with  the 
flexible  structure  of  the  program,  the  list  of  penalty 
functions  can  be  modified  or  expanded  to  take  into  account 
any  other  desirable  trajectory  characteristics. 
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PROGRAM  APPLICATIONS 


A. 3 


Sections  A.l  and  A. 2 described  the  modular  structure 
of  program  ENERGY  and  the  methods  of  determining  a cost  for 
each  trajectory,  in  order  to  compare  different  trajectory  and 
control  profiles.  The  present  section  briefly  describes  the 
major  program  modules  and  demonstrates  how  they  can  be  executed 
in  appropriate  order  to  perform  complex  trajectory  and  guidance 
algorithm  analysis. 


A list  of  the  major  computational  blocks  and  a des- 
cription of  their  functions  is  given  below: 

• Trajectory  Generation  - Given  an  ini- 
tial vehicle  state  (velocity,  altitude, 
flight  path  angles  and  mass)  and  control 
profile,  the  trajectory  generation  logic 
integrates  forward  from  the  initial  con- 
ditions the  complete  trajectory  profile. 

Either  open-loop,  ad  hoc,  or  linear- 
optimal  perturbation  guidance  algorithms 
can  be  employed. 

• Cost  Computation  - For  each  complete 
trajectory,  a scalar  cost  is  computed 
as  described  in  Section  A. 2. 

• Partial  Derivatives  - The  partial  deriva- 
tives of  the  vehicle  dynamics  and  the 
penalty  functions  with  respect  to  the 
state  and  control  variables  are  evaluated 
along  the  computed  trajectory. 

• Steepest-Descont  Control  Perturbation  - 
Using  the  conditions  and  partial  deriva- 
tives  of  the  nominal  trajectory,  this 
program  section  integrates  the  costate 
variable  equations  backward  from  the 
terminal  conditions  and  computes  the  con- 
trol perturbation  history  to  most  effectively 
reduce  trajectory  costs. 
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• Linear-Optimal  Gain  Computation  - Given  a 
nominal  trajectory  to  follow,  this  routine 
solves  the  optimal  regulator  control  problem 
to  determine  linear  perturbation  guidance 
gains.  The  solution  to  the  optimal  regula- 
tor problem  is  achieved  by  computing  the 
steady-state  values  of  a matrix  Riccati 
equation. 


Flow  diagrams  depicting  three  of  the  major  trajectory 
analysis  techniques  which  can  be  exercised  with  ENERGY  are 
presented  in  Fig.  A-2.  These  diagrams  demonstrate  how  the 
ordering  of  the  program's  computational  blocks  can  be  rearran- 
ged to  perform  significantly  different  forms  of  trajectory 
analysis. 


Many  of  the  computational  blocks  are  used  for  more 
than  one  type  of  trajectory  analysis.  The  trajectory  genera- 
tion block,  for  example,  is  used  for  steepest -descent  optimi- 
zation, for  parametric  optimization,  and  for  optimal  pe-tur- 
bation  guidance  computations.  In  each  application,  however, 
the  guidance  algorithm  is  chosen  to  fit  the  specific  needs 
of  the  problem.  For  steepest-descent  optimization,  tra- 
jectories are  generated  using  open-loop  guidance  commands 
which  are  computed  to  minimize  total  trajectory  cost.  Para- 
metric optimization  trajectories  utilize  closed-loop  guidance 
commands  to  fly  the  profiles  which  are  specified  by  the  param- 
eters to  be  optimized.  Linear-optimal  guidance  gain  analysis 
requires  trajectory  generation  with  either  open-loop  or  closed- 
loop  guidance  algorithms  to  define  the  nominal  path  and  its 
associated  linear  model.  Performance  of  the  linear-optima 
guidance  law  is  evaluated  using  both  linear  and  nonlinear  mathe- 
matical models. 

Program  ENERGY'S  applications  can  be  made  more  complex 
by  chaining  together  the  analysis  tools  depicted  in  Fig.  A-2. 
Rather  than  initiating  the  guidance  gain  computations  with  an 
arbitrary  trajectory,  the  trajectory  defined  as  a result  of 
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STEEPEST-  DESCENT  OPTIMIZATION 
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A-2  Flow  Diagrams  of  Trajectory  Analysis  Tools 


steepest -descent  or  parametric  optimization  could  be  employed. 
The  guidance  gains  then  computed  always  would  attempt  to 
follow  the  optimal  reference  trajectory.  Section  A. 4 shows 
that  ENERGY'S  control  logic  facilitates  this  type  of  chain- 
ing operation. 


A. 4 PROGRAM  CONTROL 

The  control  logic  of  ENERGY  provides  lor  both 
sequential  and  iterative  processing  of  the  major  trajectory 
analysis  algorithms.  Sequential  processing  enables  chaining 
together  specific  program  applications,  such  as  trajectory 
optimization  and  optimal  guidance  gain  computations.  Itera- 
tive program  operation  is  useful  for  examining  the  effects 
of  vehicle  or  environmental  variations  on  system  performance 
(e.g.,  optimal  trajectory  shapes,  nominal  trajectory  costs, 
on-board  guidance  algorithm  operation,  etc.). 

The  two  modes  of  program  operation  — sequential 
and  iterative  — are  controlled  by  the  program's  RUN-CASE- 
SWEEP  control  logic.  This  control  logic  is  basically  a set 
of  nested  repetitive  loops  as  shown  in  Fig.  A-3.  The  inner- 
most loop  is  the  SWEEP,  the  middle  loop  is  the  CASE,  and  the 
outer  loop  is  the  RUN  type. 

The  RUN  control  point  specifies  the  particular 
type  of  trajectory  analysis  to  be  performed,  such  as  para- 
metric trajectory  optimization.  A complete  list  of  analysis 
options  is  given  in  Section  A.l.  The  CASE  control  point 
specifies  the  desired  terminal  conditions  and  inflight  con- 
straints to  be  used  in  determining  trajectory  costs.  The 
available  penalty  functions  are  presented  in  Section  A. 2. 

The  SWEEP  control  point  provides  for  modification  of  the 
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V 


MOIM 


PROGRAM  INITIATION 


Figure  A- 3 RUN-CASE-SWEEP  Control  Logic 
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aircraft  and  environmental  conditions.  The  conditions  which 
can  be  modified  at  the  sweep  control  point  are: 

e Initial  vehicle  state:  Velocity,  altitude, 

flight  path  angle,  and  mass. 

e Aircraft  modifications:  Locations  of  vehicle 

c.g.,  number  of  engines  operating. 

e Environmental  effects:  Wind  profile  vs. 

position  and  altitude,  atmospheric  tempera- 
ture variations. 

The  use  of  sequential  analysis  program  operation 
is  achieved  by  consecutively  executing  two  RUN  types.  That 
is,  if  only  one  case  and  one  sweep  are  specified  for  each 
RUN  type,  control  looping  occurs  only  at  the  RUN  control 
point.  By  this  means,  an  optimal  guidance  gain  computation 
could  be  performed  immediately  following  a trajectory  opti- 
mization step.  The  optimized  trajectory  would,  therefore, 
be  available  as  the  nominal  profile  to  be  used  in  the  gain 
computation  analysis. 

Iterative  program  operation  is  achieved  by  specify- 
ing more  than  one  CASE  and/or  SWEEP  option  for  each  RUN 
type.  This  iterative  processing  of  multiple  CASE  or  SWEEP 
options  enables  comparisons  of  trajectory  costs  as  a function 
of  changing  vehicle  and  atmospheric  conditions. 

An  example  of  program  ENERGY  output  is  illustra- 
ted in  Fig.  A-4 , which  is  a partial  listing  of  the  200-nm 
steepest-descent  optimal  trajectory.  The  columns  of  the 
flight  path  output  are  defined  as  follows: 

ENERGY  Specific  total  energy  times  g, 

ft2/sec2 

K.E.  Specific  kinetic  energy  times  g, 

ft2/sec2 

HEIGHT  Altitude,  ft 


185 


VELOCITY 

True  airspeed,  fps 

GAMMA 

Flight  path  angle,  deg 

RANGE 

Range  to  go,  nm 

XI 

Heading  Angle,  deg 

AGO 

Azimuth  Angle,  deg 

MACH 

Mach  number 

ALPHA 

Angle  of  attack,  deg 

PHI 

Roll  angle,  deg 

L/D 

Lift-to-drag  ratio 

DYNP 

Dynamic  pressure,  psf 

GLOAD 

Load  factor,  Mg'"s 

CL 

Lift  coefficient 

CD 

Drag  coefficient 

TIME 

Elapsed  time,  sec 

IJK 

Step  number 

IHLF 

Number  of  step-size  halvings 

PITCH 

Pitch  angle,  deg 

THROTL 

Throttle  setting,  (0  to  1) 

FLOW 

Fuel  flow  rate,  lb/hr 

NRT 

Normal-rated  thrust,  lb 

WEIGHT 

Vehicle  weight,  lb 

FORCE 

Thrust , lb 

WT/DELTA 

Weight  divided  by  pressure 
ratio,  lbs 

DELTA 

Atmospheric  pressure  ratio 

HDOT 

Climb  rate,  ft/min 

SR 

Specific  range,  nm/lb  of  fuel 

IAS 

Indicated  airspeed,  kt 

The  flight  path  listing  is  followed  by  particulars 
of  the  optimization,  including  terminal  errors,  penalty 
costs  of  the  trajectory  and  dollar  costs  of  the  mission. 
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Figure  A-4  Segments  of  2 ,000-nra  Steepest -Descent  Trajectory  Listing. 

(Throttle-Down  Begins  at  1098.44  Sec)  (Continued). 


APPENDIX  B 


MATHEMATICAL  MODEL  OF  THE  C-141A  AIRCRAFT 


The  investigation  of  aircraft  fuel  minimization 
techniques  has  concentrated  on  the  C-141A  aircraft,  and 
mathematical  models  of  the  C-141A  have  been  included  in 
the  computer  program  ENERGY  which  performed  all  the  trajec- 
tory analysis  related  to  this  task. 

Five  sources  of  vehicle  data  have  been  used  for 
defining  the  vehicle  model  (Refs.  1,  2,  32,  47,  and  48). 
These  sources  provide  what  is  believed  to  be  a vehicle  model 
with  a high  degree  of  fidelity.  The  appropriate  equations 
and  the  data  which  comprise  the  model  are  presented  in  the 
following  sections  of  this  appendix. 


B.l  C-H1A  AERODYNAMIC  MODEL 

The  aerodynamic  data  model  for  the  C-141A  aircraft 
has  been  obtained  from  Ref.  48.  The  model  computes  the  air- 
craft's coefficients  of  lift  and  drag,  C^  and  C^,  as  func- 
tions of  angle  of  attack  (a),  Mach  number  (M),  altitude  (h), 
center  of  gravity  (c.g.),  total  engine  thrust  (F),  and 
dynamic  pressure  (q).  Because  this  study  employs  only 
throttle  setting  and  pitch  angle  as  guidance  commands,  the 
vehicle's  flap  and  spoiler  settings  are  assumed  to  be  zero. 
As  only  two-degree-of-f reeuom  dynamics  are  used  for  optimi- 
zation, the  total  pitching  moment  on  the  vehicle  also  is 
assumed  to  be  zero,  i.e. , the  vehicle  is  trimmed. 
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With  the  above  assumptions,  the  computations  needed 
to  compute  and 
coefficient  is  computed  as 

CLBAS  ■ CLALP*CRCLA* ( a-CRALO*aQ)  (B.l-1) 


are  as  follows.  First,  the  basic  lift 


where 


CLALP  (M)  ■ lift  curve  slope 
CRCLA()(,h)  ■ lift  flexibility  parameter 
aQ(M)  - zero-lift  angle  of  attack 
CRALO(M,h)  - zero-lift  flexibility  parameter. 

The  total  lift  coefficient  is  then  given  by 

CL  « CLBAS  - 0. 2979(CMIT*STAB'  (B.l-2) 


where 


CMIT  *=  pitch  moment  coefficient  due  to 
stabilator  deflection 

STAB  ■ Stabilator  deflection  angle 

The  coefficient  0.2979  is  the  ratio  of  the  aircraft Ts  reference 
chord  to  tail  moment  arm  and  reflects  the  fact  that  the  tail 
provides  lift  as  well  as  pitching  moment. 

The  zero  total  pitching  moment  assumption  allows 
the  quantity  CMIT  • STAB  to  be  expressed  in  terms  of  the  other 
pitching  moments  acting  on  the  vehicle.  Making  these  sub- 
stitutions results  in 

p CLBAS  ( 1+0 . 2979 » CMCL » CRCMC ) +0 . 2979  ( CMO  * CRCMO+CMTH ) 

UL  [1-0. 2979(CCG-0. 241)] 


(B. 1-3) 
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where 


CMCL(M)  ■ pitching  moment  coefficient  due  to  lift 
CRCMC(M.h)  * CMCL  flexibility  parameter 
CMO(M)  * pitching  moment 
CRCMO(M,H)  * CMO  flexbility  parameter 

CMTH  * pitch  moment  due  to  thrust 

The  thrust-induced  pitching  moment,  CMTH,  which  ap- 
pears in  Eq.  (B.l-3)  is  computed  from  the  engine  thrust  and 
dynamic  pressure  as 


CMTH  * F/(17960*q ) (B.l-4) 

Having  computed  the  vehicle's  lift  coefficient,  the  drag 
coefficient  is  determined  as  a function  of  CL  and  M,  i.e., 

Cp  - (CDCCL(M,Cl))  (B. 1-5) 

Either  total  lift  coefficient  (CL)  or  wing  lift 
coefficient  (CLBAS)  could  be  used  to  compute  the  drag 
coefficient  (Eq.  (B.l-5)).  The  first  approach  neglects 
trim  drag,  as  the  down  load  on  the  tail  (required  for  static 
trim)  is  not  considered  explicitly.  Since  CLBAS  is  greater 
than  CL,  trim  drag  is  included  in  the  second  approach, 
has  been  computed  using  CL  for  optimal  and  conventional 
flight  path  computations  and  using  CLBAS  in  the  specific 
range  analyses. 

The  computations  of  vehicle  C^  and  Cp  described  in 
Eqs.  (B.l-1),  (B.l-3)  and  (B.l-5)  necessitate  performing 
table  lookups  to  obtain  appropriate  values  of  the  tabulated 
data.  The  table  lookups  are  performed  using  either  linear 
or  cubic  polynominal  internolation.  The  linear  lookup  rou- 
tines provide  satisfactory  results  for  all  trajectory 
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simulations.  The  cubic  interpolations  are  employed  whenever 
partial  derivatives  of  the  data  are  required,  as  this  infor- 
mation is  available  with  the  cubic  interpolation  routines. 

The  actual  data  values  used  in  the  computation  of 
vehicle  lift  and  drag  are  in  the  form  of  single-  and  two- 
dimensional  data  arrays.  To  complete  the  C-141A  aerodynamic 
model  specification,  these  data  values  are  presented  below. 

The  data  tables  are  stored  as  functions  of  Mach 
number,  altitude,  and  CL.  The  arrays  which  define  the  values 
of  the  Mach  number-only  dependent  variables  are: 

• Lift  curve  slope  (CLALP(M)) 

• Zero  lift  angle-of-attack  (aQ(M)) 

• Pitching  moment  coefficient  due  to  lift  (CMCL(M)) 

• Pitching  moment  (CMO(M)) 

The  data  points  for  these  variables  as  a function  of  Mach 
number  are  shown  in  Table  B-l. 
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TABLE  B-l 


MACH  NUMBER-DEPENDENT  AERODYNAMIC  VARIABLES 


Mach  number 

CLALB 

a0 

CMCL 

CMO 

0.0 

0.0835 

-3.  37 

-0.246 

0.085 

0. 1 

0.0851 

-3.39 

-0.245 

0.084 

0.2 

0.0867 

-3.41 

-0.244 

0.083 

0.3 

0.0882 

-3.43 

-0.241 

0.080 

0.4 

0.0900 

-3.45 

-0.242 

0.080 

0.5 

0.0917 

-3.46 

-0.247 

0.082 

0.6 

0.0937 

-3.48 

-0.254 

0.085 

0.65 

0.0968 

-3.42 

-0.257 

0.089 

0.7 

0.0998 

-3.36 

-0.260 

0.092 

0.75 

0.1076 

-3.27 

-0.255 

0.098 

0.  77 

0.1100 

-3.22 

-0.263 

0.096 

0.8 

0.1055 

-3.25 

-0.283 

0.100 

0.83 

0.0870 

-3.28 

-0.311 

0.080 

0.86 

0.0861 

-2.56 

-0.317 

0.061 

0.9 

0.1030 

-1.30 

-0.293 

0.115 

The  data  points  for  those  variables  which  are  func- 
tions of  two  independent  variables  are  the  following: 


• Lift  flexibility  parameter  (CRCLA(M,h)) 


\h(ft) 

U \ 

20,000 

40,000 

0.0 

1.00 

1.00 

1.00 

0.1 

1.00 

1.00 

1.00 

0.2 

.983 

.992 

.994 

0.3 

.949 

.976 

.987 

0.4 

.900 

.952 

.976 

0.5 

.854 

.926 

. 966 

0.6 

.806 

.995 

.955 

0.65 

.775 

.870 

.944 

0.70 

.744 

.845 

.933 

0.75 

.712 

.816 

.912 

0.77 

.712 

.804 

.893 

0.80 

.712 

.796 

.883 

0.83 

.712 

.810 

.905 

0.86 

.712 

.814 

.994 

0.90 

.712 

.749 

.873 
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CMO  flexibility  parameter  (CRCMOfU.h 


20,000  40,000 


l.'OO 

1.00 

.088 

.928 

.860 

.756 

.670 

.639 

.608 

.550 

.520 

.490 

.440 

.466 

.496 


Coefficient  (CDCCL(M,CL)) 


DDDBI 


.6  0.7 


.0234 
.0353 
.0498 
0595  I . 0610 


.0213  | 

.0253 

.0215 

.0253 

.0213 

.0253 

.0213 

.0253 

.0214 

.0254 

.0216 

.0257 

.0221 

.0262 

.0224 

.0266 

.0228 

.0270 

.0234 

.0276 

.0240 

.0283 

.0267 

.0320 

.0393 

.0476 

.0536 

.0610 

.0637 

.0697 

,0304 
,0308 
,0313 
,0318 
.0323 
0333 
,0347  1.0493 
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B. 2 C-141A  THRUST  MODEL 

The  installed  thrust  model  of  the  C-141A  aircraft 
provides  the  normal  ratud  thrust,  NRT,  of  each  of  the  air- 
craft's engines.  The  model  included  in  program  ENERGY  was 
provided  by  the  fuel  management  studies  of  Ref.  2. 

Altitude,  Mach  number,  and  non-standard  atmospheric 
temperature  effects  on  NRT  are  accounted  for  in  the  current 
engine  model.  A thirteen-term  polynomial,  which  is  a func- 
tion oi  these  parameters,  is  evaluated  to  determine  engine 
thrust.  The  NRT  equation  is  given  by 

NRT  - 5763.499  - 498.8629*  UN  - 5590.819-HN 

- 123. 8276*UN2-136. 1916*  UN*MN 

+ 77.995*  MN2  + 1177. 55* UN-IIN+731 . 959-MN-HN 

+ 1609. 4427* HN2  + 108. 217*UN*MN*HN 

- 875. 1862* UN* HN2  - 805. 6959* MN*HN2 

+ 385.875*  HN3  (B. 2-1) 

where 

NRT  * normal-rated  thrust  (lb) 

UN  ■ normalized  temperature  deviation 
HN  * normalized  altitude 
MN  = normalized  Mach  number 

The  normalized  terms  in  Eq.  (B.2-1)  are  computed  as: 
UN  » ATEMP/40.  (B.2-2) 

HN  - (h-25500.  V29500.  (B.2-3) 
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MN  - (M-0. 55)/0. 25 


(B.2-4) 


where  ATEMP  is  the  tenperature  deviation  (°F)  from  the  1962 
Standard  Atmosphere,  h is  the  aircraft  altitude  in  feet,  and 
M is  the  aircraft's  Mach  number. 

The  polynomial  form  of  the  NRT  computation  eliminates 
the  need  to  execute  numerous  table  Interpolations.  It  also 
simplifies  the  evaluation  of  the  partial  derivatives  of 
engine  thrust  with  respect  to  the  vehicle  state,  which  is 
needed  for  trajectory  and  guidance  gain  optimization.  The 
partial  derivatives  of  NRT  with  respect  to  the  vehicle  state 
are  computed  in  closed  form  by  taking  the  appropriate  partial 
derivatives  of  Eq.  (B.2-1). 


B.  3 ENGINE  FUEL  FLOW  MODEL 

A general  fuel  flow  model  for  the  C-141A  was  obtained 
from  Ref.  1.  The  data  represents  flow  rates  for  a TF33 
Pratt  and  Whitney  engine  installed  on  the  C-141A  aircraft. 

Figure  B-l  presents  the  fuel  flow  data  in  graphical 
form.  In  this  figure,  normalized  flow  rato  is  plotted  as  a 
function  of  normalized  engine  thrust  for  three  Mach  numbers. 

The  atmospheric  normalization  factor,  6,  is  the  ratio  of 
atmospheric  pressure  at  altitude,  p^,  to  pressure  at  sea-level, 
PSL,  i. e. , 


6 “ Ph/PSL  (B. 3-1 ) 

The  normalization  factor,  0,  is  given  by 

0 - TEMPh/TEMPgL  (B.3-2) 
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where 


TEIIPh  - temperature  at  altitude 
TEMPgL  “ sea-level  temperature 


For  the  analysis  of  this  study,  engine  thrust  is  computed 
simply  as  the  product  of  throttle  setting,  T,  and  the 
engines ' NRT , 


F - T • NRT  (B. 3-3 ) 

where  throttle  setting  ranges  between  0 and  1. 


TF33  InttolUd  Eneina  Do  to* 


Figure  B-l  TF33  Installed  Engine  Data* 


*Fuel  flow  rates  modeled  in  program  ENERGY  have  been  increased 
by  11  percent  to  match  C-141A  Category  II  Flight  Performance 
Test  results. 
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The  data  presented  in  Fig.  B-l  is  incorporated  in 
program  ENERGY  as  a two-dimensional  table  lookup.  Fuel  flow 
data  is  stored  as  a function  of  engine  thrust  and  Mach  number. 

As  with  the  vehicle  aerodynamic  data,  linear  table  interpolation 
is  used  for  most  trajectory  analysis.  For  those  functions 
which  require  partial  derivatives  of  fuel  flow  rate  with 
respect  to  the  vehicle  state,  cubic  interpolation  routines 
are  employed  to  accurately  determine  the  required  derivatives. 

In  order  to  validate  the  C-141A  vehicle  model  of  pro- 
gram ENERGY,  the  maximum  range  factor  and  best  Mach  number  were 
computed  for  a wide  range  of  aircraft  altitudes  and  weights. 

The  results  of  this  study  then  were  compared  with  the  data 
available  from  the  C-141A  Category  II  Performance  Tests  (Ref. 
47).  It  was  found  that  the  range  factor  performance  of  the 
simulated  vehicle  was  approximately  11  percent  too  high.  The 
fuel  flow  rates  depicted  in  Fig.  B-l  were  Increased  by  this 
percentage  and  the  range  factor  comparisons  repeated.  The 
close  agreement  between  the  modeled  results  and  those  of  the 
vehicle  tests  is  shown  in  Fig.  B-2.  The  trajectory  analysis 
parformed  in  the  current  study  employed  the  fuel  flow  data 
adjusted  by  the  11  percent  factor  in  order  to  represent  an 
actual  vehicle  as  accurately  as  possible. 


B.  4 VEHICLE  LOADING  CONDITIONS 

Typical  C-141A  vehicle  weights  and  fuel  loadings  are 
documented  in  Ref.  32.  Those  values  which  were  useful  for 
the  current  program  development  are  listed  here: 

e Maximum  ramp  weight:  318,000  lb 

e Maximum  normal  landing  weight:  257,500  lb 

e Maximum  flight  fuel  load:  148,120  lb 
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1-20144 


Figure  B-2  Comparison  of  Vehicle  Model  and 
Actual  Aircraft  Performance 

The  vehicle  reference  area,  needed  for  computations 
of  specific  force  on  the  vehicle  was  obtained  from  Ref.  45 
as  3,228  ft^.  At  a ir  xirnum  vehicle  weight  of  318,000  lb, 
level  "1-8"  flight  would  necessitate  a wing  loading  on  the 
order  of  99  lb/ft^. 
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